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SOAP USER'S GUIDE 
1. INTRODUCTION 

The Subsatellite Orbital Analysis Program (SOAP) Is an all-digital functional 
simulation designed to model the Interaction of several Independent fliers In 
the on-orbit environment. Equations of motion are Integrated for both rota- 
tion and translation of each vehicle. Each has its own pilot, control system, 
sensors, and environmental interactions. The properties of each vehicle and 
the fidelity with which each is modeled are determined by the user, The SOAP 
supplies a library of math models for various vehicles, a structure for setting 
up the relationships among various math models, and a mechanism for general 
input and output. 

The SOAP is written in FORTRAN and is designed for use with a UNIVAC 1 1 00 
series computer. It differs from most FORTRAN simulations in that the order 
In which computations are made is determined not only by logical branching 
within subroutines or within a driver program, but on the basis of simulation 
time, The SOAP maintains an Internal clock and also maintains a schedule of 
events. These events can Include such procedures as execution of control 
system routines, sampling the vehicle dynamic state by sensors, response of 
control jets to a pilot's commands, the recording of vehicle states for 
output, or the pilot checking the status of a fuel gauge. The user need only 
develop subroutines modeling the desired events and specify their whereabouts 
to the SOAP executive program. Events which occur at regular time intervals 
are scheduled by alerting the SOAP executive. Routines which occur only 
irregularly or which occur only when some specified condition is met can be 
scheduled or called by other routines during the execution of the simulation. 

The advantage in using the SOAP executive lies in the fact that the sequence 
of events can be arranged to reflect realistic order and time of occurence. 
There is no need for the times at which events occur to be integer multiples 
of a single time unit. In addition, the SOAP automatically creates a driver 
to call all the subroutines desired and to produce appropriate input and output. 
This facilitates the setup and execution of new simulations. 


From the user's point of view, each vehicle has a separate pilot, control 
system, dynamics, and sensors. Each vehicle and its pilot may or may not 
have knowledge of the activities of any other. Typically, communication 
between the control system, pilot, and dynamics takes place by Issuing com- 
mands in the form of setting locations In COMMON blocks. Communication with 
the user typically takes the form of inputing into locations in COMMON blocks 
at the start of simulation, printing out information during the simulation 
(either at fixed time intervals or at the beginning and end), and of plots 
generated after the simulation. 


Performing a simulation using the SOAP typically Involves four stages of 
program execution: 

a. Preprocessing - The math models to be used, the logical and/or temporal 
sequence of events, and the input/output desired are specified, A pro- 
gram is executed which sets up the driver program and writes several sub- 
routines which %ci as subdrivers. 

b. Mapping or Object Code Collecting - Instructions are provided which de- 
scribe the location of all object code which will be needed by the 
simulation. 

c. Simulation Time Proceeds - Input to the simulation is typically performed 
by specifying values for COMMON block locations at the start of the simu- 
lation. The simulation proceeds until it is terminated due to the attain- 
ment of some internal condition. Printed output can be obtained during 
this stage. 

d. Postprocessing - Data which has been collected during the simulation is 
plotted. 

From the user's point of view, steps a, b, and d are not of particular interest. 
They simply require a series of steps which must be followed in order to pro- 
duce the simulation desired. They are, therefore, documented by a user's 
guide (section 2) and by brief descriptions of individual subroutines. 

Step c, the simulation, is discussed further. 

The SOAP system maintains a clock and a driver program. It takes care of 
input and output automatically, but most other activities are specified by 
the user. These routines, typically, are either scheduled or called by other 
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routines. The SOAP system is not aware of the names of the routines until 
the preprocessing stage. 

Several activities have a special status within a SOAP simulation, One such 
activity is the process of computing the position and attitude of each vehicle. 
This activity requires that forces and torques due to such effects as gravity, 
gravity gradient, aerodynamic drag, and reaction control jets be computed. 

It also requires that the mass and the inertia tensor of each vehicle be 
updated. This series of computations is termed the "Environment" (of the 
Flight Control System and pilot), In most situations the entire dynamics 
computation is performed at once (for all vehicles) , The modeling of each 
vehicle need not, however, be the same. Both the causes of force and torque 
and the way in which the causes are modeled can differ. 

The Environment communicates with the rest of the simulation and the user via 
COMMON blocks. The Environment routines and COMMON blocks have a special status 
within the SOAP, wecause the SOAP executive knows their names. This knowledge 
of the names of the COMMON blocks Involved is required to enable the user to 
input data from the outside. The knowledge of the names of model s allows the 
SOAP executive to create a driver. It should be noted that any of these names 
can be changed, but a system change is required. 

Figure 1-1 displays the organization of the Environment routines. The names 
of the routines displayed are really names of driver routines. Each driver 
does 1 ittle except call the routines which perform the computation for each 
vehicle. The usual execution involves the driver program ACTVEH calling 

ar • 

vehicle models for each vehicle (Space Shuttle Grbiter, Subsatellite, Manned 
Maneuvering Unit, etc. ). These vehicle models typically call all of the force 
and torque computations involved with a single vehicle. 

Two other special categories of model exist: "Passive Vehicles" which repre- 

sents a free-flying vehicle without active control, and Sensors. Each of 
these categories is special to the SOAP executive in that the routine name 
and supporting COMMON block: name is recognized. The Sensor block consists 
of a sensor control driver which, in turn, calls model s for the sensors 
(e.g. , gyros, horizon sensors) associated with each vehicle. 


^ C tcTU 9 r 0rCe and torque coniputers » integrates equations of motion) COWON blocks, ACTVEH, ACTV1C, 
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All routines In the Environment, Sensors, or Passive Vehicles are specified 
to the Preprocessor by listing the version names of those drivers which are 
to be used to the Preprocessor (see section 2). All routines used, both 
drivers and actual models, must have the relevant object code collected by 
mapping Instructions during the map execution (phase b above), The actual 
execution of these routines is triggered by the necessity of updating the 
dynamics (as requested from other parts of the simulation) or of carrying 
out a command (e.g. , turn on jets) Issued by the flight control or pilot. 

The dynamic state can also be updated at regular intervals by scheduling a 
driver program termed CALENV at the desired Intervals. 

The SOAP differs from previous similar simulation programs (notably the Space 
Shuttle Functional Simulator (SSFS)) in that the computations supporting the 
vehicle dynamics can be performed entirely In double precision. 

All portions of the simulation not previously specified are specified by the 
user at the time of Preprocessing (step a) and/or mapping (step b). 

Any operation desired by the user can be written as a subroutine, either 
callable by another routine or occurring on the basis of time or logical 
decision. Any routine which is not solely called by other routines is 
"scheduled." That Is, the user tells the SOAP executive the name of a particu- 
lar simulation subroutine, the simulation time at which it begins to be used, 
the time Interval between executions, and the time when it ceases to be used, 
If any routine Is to be executed before or after a scheduled routine but at . 
the same simulation time, this can also be specified. The preprocessor uses 
this information to create a calling routine (SELECT) with knowledge of the 
routines to be called. The map Instructions must include the whereabouts of 
object code for all routines used, whether scheduled or called out, 

Some models, such as control jet models, accept commands from the pilot or 
control system and may even delay the execution of the action commanded. Such 
a delay would be specified by the environment model accepting the commands. 

The SOAP master clock ensures that the delay is .applied, Table 1-1 lists 
commands which are in use by some of the SOAP models. The SOAP system assigns 
the usage of these command indices at the time of simulation, so there is no 
necessity to avoid the command Indices used by a model when the model is not 
In use, 


TABLE 1-1.- SOAP ENVIRONMENT COMMAND SERIES ASSIGNMENTS 


Command nosJ 

Mode 

500 - 550 

SSO - RCS Model 

551 - 574 

RCS free flier 1, 24- jet system 

575 - 599 

RCS free flier 2, 24-jet system 


^Command numbers TOO - 1100 are available 


Figure 1-2 displays the organization of a simple two-vehicle simulation. 

The two vehicles involved are the Maneuverable Television (MTV) and the 
Space Shuttle Orbiter (SSO). The interaction between the two in the simu- 
lation displayed is solely through the "master clock" which represents the 
entire SOAP executive and scheduling system. The information written near 
the interface arrows describes some of the variables, in COMMON blocks, which 
are shared by the models. Anyone of these variables is available for output 
at regular time intervals by specification to the SOAP system executive at 
the preprocessing stage. Such operations are transparent to the user once 
they are specified. 

The simulation in figure 1-2 shows no sensors or passive vehicles, but these 
would form separate, additional blocks driven by the master clock. Pilots 
for the vehicles would form additional blocks within flight control. Infor- 
mation about each vehicle can be passed to any other through the COMMON blocks. 
Such passage would involve additional lines joining the flight control blocks 
of the vehicles. Vehicles can be added simply by replicating the blocks which 
model the vehicle features and by supplying the appropriate vehicle parameters. 

Output from a SOAP simulation is typically initiated by the user by specifying 
those variables to be printed or plotted and the time interval at which they 
are to be recorded. The specifications are made during the preprocessing 
operation. Up to six dependent variables can be displayed as a function of 
a single "independent" variable within each plot. Three of the variables are 
represented by the position of symbols of the user's choosing, and three by 
the angles at which these symbols are drawn. Simplified plots can be produced 
on a line printer. 
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The SOAP executive provides diagnostic capabilities in the form of a state 
image print and a restart capability. The state, image print reproduces all 
of the COMMON block locations in decimal or Hollerith form. Such a print is 
typically supplied at the start and end of a simulation and can be requested 
(scheduled) at intervals during the simulation. The restart capability allows 
the user to initialize the simulation to replicate some point in the middle 
of another simulation; thus, a problem maneuver could be restarted immediately 
prior to a point of difficulty to avoid repeated simulation of an uninteresting 
time period, Instructions for obtaining a restarted run can be found in 
section 2.4.2. 

Up to this point the SOAP description has been deliberately vague. The 
simulations which can be performed using the SOAP structure are not clearly 
limited. Those simulations which have already been performed and for which 
runstreams and models are already set up form a much more limited set. 

These simulations are documented by the runstreams in section 2 and by the 
models in section 4. It is hoped that these models will give the user a 
headstart in creating a model for any operations required. 


Models of vehicle flight control systems are typically documented by flowcharts, 
scheduling algorithms, and program descriptions. The SSO software has, how- 
ever, a special status. The operation of the SSO software is described in 
detail by Functional System Software Requirements (FSSR) documents. Much 
of the SSO software used by SOAP was developed for the SSFS and was designed 
to produce a functional model of the flight system. These routines, therefore, 
are not as extensively documented as are those routines developed for the 
SOAP alone. 


The SOAP has roots in the SSFS. This document itself owes much of its develop- 
ment to documentation for the SSFS (refs. 1, 2, and 3). The authors would 
like to thank those who developed the SOAP framework and to refer the reader 
to this documentation as an example of a SOAP-type simulator. 
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2, SOAP USER'S GUIDE 


2.1 PROGRAM FILES 

This guide describes construction of runstreams which cause all four SOAP 
operations* i .e. 

Preprocessing 

Mapping 

Simulation Execution 
Postprocessing 

These instructions are geared for the user of the UNIVAC 1110 with Exec 8. 

The SOAP systems and library routines are maintained on tapes and secure files 
in the UNIVAC 1110. Both source and relocatable codes are available, though 
the casual user will interact with these routines only by scheduling and by 
mapping the relocatable elements. Typical selections of models and mapping 
instructions are available in section 2.2. 

2.2 RUNSTREAM 

Table 2-1 displays the basic logical unit assignments for a SOAP run. Table 
2 - I I describes the setup of the files required for a typical simulation. 
Typical lengths are associated with these files to ensure the availability 
of adequate storage. 


TABLE 2-1.- INPUT/OUTPUT FILE REQUIREMENTS FOR SOAP EXECUTION 




(a) Preprocessor 

Usage 

Jm. 

Purpose 

Required 

File 

Output images for subroutine SELECT, 
CALMOD, and DUMMYS 

Required 

Fite 

Output images for the collector directives 

Required 

Reader 

Card input (can be file) 

Required 

Printer 

Printed output (can be file) 

Required 

File 

Scratch input/ output 



(b) Simulation 

Usage 

1m. 

Purpose 

Required 

File 

Scratch input/output 

Required 

Reader 

Card input (can be file) 

Required 

Printer 

Printed output (can be file) 

Optional 

Tape 

State image dump input file; used for 
RESTART runs only 

Optional 

Tape 

State image dump output file; used for 
RESTART or RECYCLE runs 

Optional 

File 

Headlines print retention file; required 
if subroutine HDLNS is used 

Optional 

File 

On-Orbit Astronaut File 

Optional 

File 

Standard or primary plot dump unit; required 
if subroutine BLKPLT is used 

Required 

File 

Scratch input/output 

Optional 

File 

Secondary plot, dump unit 
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TABLE 2-1.- CONCLUDED 


(c) Plotter 


LU 

Usage 


Purpose 

3 

Required 

File 

Scratch input/output 

5 

Required 

Reader 

Card input 

6 

Required 

Printer 

Printed output 

14 

Optional 

Tape 

Plotter output tape required for Cal Comp 
plots 

20 

Optional 

File 

Standard or primary plot data input unit 
as furnished by the simulation step 

29 

Optional * 

File 

Intermediate storage of plot data; required 
if the data is provided by tape and the 
PLOTER program is used 

n 

Optional 

File 

Secondary plot dump file, as furnished by 
the simulation step 


TABLE 2- II.» SOAP RUNSTREAM SETUP 







[Men multiple entries with the same number] 
[appear, only one is used for a simulation J 



ORUN.R/TPR 170JG, E/3033, ED4-L44444,60,800 
OHDG.P PAGE HEADING INDENTIFYING RUN 
OASG.AX ED4 -L44444*S0AP800. 

OUSE S., ED4-L44444*S0AP800. 

OASG.AX ED-L44444*S0APW0RK» 

OUSE SW., ED4-L44444*SOAPWORK. 

OUSE $5 . , ED4-L44444*S0AP805S1 . 

OUSE S5. , ED4-L44444*S0AP805D2 . 

0ASG,AX S5. 

OUSE SB., ED4-L44444*S0AP2B0DS . 

OUSE SB., ED4-L44444*S0AP2B0DD , 

OUSE SB., ED-L44444*S0AP3B0DD . 

OASG.AX SB. 

OASG.T 3. 4 F40/1/TRK/20 
@ASG,T 4., FlO/l/TRK/20 
OASG , T 8., F40/5/TRK/500 
@ASG,T 20., F40/20/TRK/500 
OASG,T 29. , F40/1/TRK/20 
OXQT S5. PREPRO 
OADD.L SW.SCHED/SOAP 
OADD.L SB.SCHED/MTV 
OADD.L SB.SCHED/MTVD 
OADD.L SB.SCHED/MECD 
OADD.L SB.SCHED/MEC2 
OADD.L SB.SCHED/SSO 
OADD.P 4. 

OED.U SW.MAPIN, SW.MAPINN 
ADD SB. MAP/MTV 
ADD SB.MAP/MTVD 
ADD SB.MAP/MECD 
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TABLE 2-1 I.- CONCLUDED 


| 21 E: ADD SB. MAP/SSO 

I 22: 0PREP SW. 

\ 23: 0MAP,S SW. MAP INN, TPF$. SOAP 

24: @ v .tJT SOAP 

^ 25: (SIMULATION INPUT DATA CARDS - see section 2.4) 

26: @JSC*CALLUP .TAPELABEL . 

27: @ASG,T0 $P-P$TAPE,U,,21__._FRM02/31 PLT DESC. 

28: OUSE 14., $P-P$TAPE, 

29: ©REWIND 14. 

30: OXQT S5.PL0TER 

31: (PLOTTER DATA CARDS - see section 2.5) 

32: ©FIN 




A description of each of the cards shown In table 2-H follows. A more 
detailed description of the control statements- and their options can be 
found in UNIVAC system documentation. In the following descriptions, 
each “card" represents one line of Instructions. 

Card 1: 


N 


Image: GRUN.P/O RunID,AcctID,ProjID,RT,PG/CD NAME 

Purpose: To initiate the run and provide identification and accounting 
information to the system. 


Format: 


mm 

p 


- Required mnemonics for the run card. 

- Priority code for the run. It consists of a single 

letter: R, S, U, V, or Z. R has the highest priority; 

Z is the system default value. 

0 - Run card options. 

RunID - Run identification field, which must consist of the box 
number, two of the user’s initials, and a letter designa- 
tor as required for uniqueness; e.g., 170NBB. Six 
characters is the maximum number. 

AcctID - Account identification field, which is composed of the 
first letter of the user's organization code and the 
user's project number separated by a slash; e.g., E/3096. 

ProjID - Project identification field, which consists of the user's 
organization code and employee number separated by a 
minus-sign; e.g., EH2-1.55555. 

RT - Maximum run time in minutes. 

PG/CD - Maximum number of pages (PG) and cards (CD) output for 
this run. If no cards are expected, the /CD should be 
omi tted . 

NAME - First six characters of the user's last name. The field 
begins in column 61. (optional) 
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N - A 6 In column 80, designating UNIVAC 1110-6 EXEC 8 
processing. 

Card 2 ; 

Image: (SHDG.P MESSAGE 

Purpose: To provide a heading message to be printed at the top of each page 
of printer output along with a page number and date. 

Format: OHDG - Required mnemonics for the heading control statement. 

P - Heading statement option specifying that the page count 

will begin with 1. 

MESSAGE - Heading that is to appear at the top of each page. The 
heading is limited to 96 characters, including blanks, 
NOTE: A period (.) cannot be used in the message because 
the portion of the message following the period may not 
be printed, 

Card 3 : 

Image: ?ASG,A ED4-L44444*S0AP800. 

Purpose: To assign the SOAP System to this run, 

Format: @ASG - Required mnemonics for the assign card. 

A - Assign card option: The A indicates that the assign 

is for an existing cataloged file. 

ED4-L44444* - Complete file name of the cataloged file containing 
S0AP800. the SOAP System. 

Card 4 : 

Image: @USE S. , ED4-L44444*S0AP800. 

• • % , 

Purpose: To link the internal file name S. to the external file name 

ED4-L44444*S0AP800. The MAPIN program references file S. as 
the standard program file for the SOAP system. 
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- Required mnemonics for the use control statement* 

- Internal file name* 

% 

- External file name. 


Format: OUSE 

S. 

ED4- 

Card 5 : 

Image: OASG.A ED4-L44444*S0APW0RK 

Purpose: To assign the file containing flight software. Environment, and 
utility routines used for most simulations* 

Format: OASG - Required mnemonics for the assign card. 

A - Assign card option Indicating an existing cataloged 

file. 

X - Indicates exclusive use of this file by this run. 

£04-144444* - Complete file name 
SOAPWORK 

Card 6 : 

Image: OUSE W. , ED4-L44444*S0APW0RK. 

Purpose: To link the internal file name W. to the external file name 
ED4-L44444*S0APW0RK, The MAPIN program refers to 
file W. as the standard program file. 

Format: OUSE - Required mnemonics for the use control statement. 

W, - Internal file name. 

ED4~ • • • " External file name. 
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Card 7 ; 
Image: 

Purpose: 

Format: 


Card 8 : 

Image: 

Purpose: 

Format: 


Card 9 : 
Image: 

Purpose: 




©USE,S5,ED4-L44444*S0AP805S1, 

©USE ,S5,ED4-L44444*SQAP805D2 , 

To link the internal file name S5. to the external file name 
ED-L44444*S0AP805S1 or ED4-t44444*S0AP805D2, S0AP805S1 Is the 
single precision, two-vehicle program file containing the pre- 
processor element. S0AP805D2 is the double precision, myltl- 
vehicle (up to 1 SSO and 2 MTVs) program file version file version 
of S0AP805S1. 

©USE - Required mnemonics for the use control statement. 

S5. - Internal file name. 

ED4-L44444*S0AP805Si, 1 

or > external file name. 
ED4~L44444*S0AP805D2 . J 

$ASG,A S5. 

To assign the file containing the PREPRO and PLOTER elements 
and various input/output routines. 

©ASG - Required mnemonics for the assign card. 

A - Assign card option indicating an existing cataloged file. 

S5. - Internal file name attached to either ED4-L.44444*S0AP805S1 

or ED4-L44444*S0AP805D2. 

©USE SB., ED4-L44444*S0AP2B0DS . 

©USE SB., ED4-L44444*S0AP2B0DD. 

©USE SB., ED4-L44444*S0AP3B0DD . 

To link the internal filename SB. to one of three external filename 
as follows: 
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A) ED4-L44444*S0AP2B0DS - Program file containing environment and 
flight software routines for a single precision two-vehicle 
simulation. 

B) ED4- L44444*S0AP2B0DD - Same as A) only in double precision. 

C) ED4-L44444*S0AP3B0DD - Same as B) only for 3 vehicles Instead of 2, 

Format: @USE - Required mnemonics for the use control statement, 

SB. - Internal file name 

ED4- L44444*S0AP2B0DS I 

> external file, use one 
ED4- L44444*S0AP2B0DD 1 

ED4-L44444*S0AP2B0DD ' 

Card 10 : 

Image : @ASG, A SB. 

Purpose: To assign the file containing the appropriate flight software and 
environment routines for the simulation. 

Format: @ASG - Required mnemonics for the assign card. 

A * Assign card option indicating and existing cataloged file. 

SB. - Internal filename attached to one of the following: 

ED4- L44444*S0AP2B0DS . 

ED4-L44444*S0AP2B0DD. 

ED4-L44444*S0AP3B0DD . 

Card l it 

Image: @ASG,T 3. .F40/1/TRK/10 

Purpose: To assign a FASTRAND-formatted mass storage file to this run, The 

file is used for Preprocessor output. 


Format: (MSG 

T 
3. 


- Required mnemonics for the assign card, 

- Assign card option denoting a temporary file assignment. 
-Internal file name assigned to the file, 


F40 - Code specifying that this assign statement is for a 

FASTRAND- formatted mass storage file and identifying 
the specific device required. 

1 - Integer specifying the minimum amount of mass storage 

required for this file, 

TRK - Increment or granule size for this assignment. One 

track (TRK) has 1792 words, 

10 - Integer specifying the maximum amount of mass storage 

required for this file (10 tracks 6 17,920 words). 

Card 12 : 

Image: (MSG,T 4, ,F40/ l/TRK/10 

Purpose: To assign a FASTRAND-formatted mass storage file to this run. The 

file is used for Preprocessor output. 


Format: @ASG - Required mnemonics for the assign card. 

T - Assign card option denoting a temporary file assign- 

ment. 

* 

4. - Intel. .a 1 file name assigned to the file. 

F40 - Code specifying that this assign statement is for a 

FASTRAND-formatted mass storage file and identifying 
the specific device required. 

1 - Integer specifying the minimum amount of mass storage 

required for this file. 

TRK - Increment or granule size for this assignment. One 

track (TRK) has 1792 words. 

10 - Integer specifying the maximum amount of mass storage 

required for this file (10 tracks = 17,920 words). 


Image: @ASG,T 8. .F40/5/TRK/500 

Purpose: To assign a FASTRAND- formatted mass storage file to this run. The 
file is used for state image dumps for recycle and/or restart 
purposes. 


Format: 

(?ASG 

- Required mnemonics for the assign card. 


T 

- Assign card option denoting a temporary file assignment 


8 . 

- Internal file name assigned to the file. 


F40 

- Code specifying that this assign statement is for a 
FASTRAND-formatted mass storage file and identifying 
the specific device required. 


5 

- Integer specifying the minimum amount of mass storage 
required for this file. 


TRK 

- Increment or granule size for this assignment. 


500 

- Integer specifying the maximum amount of mass storage 
required for this file. 

Card 14: 



Image: 

OASG.T 20. 

, F4Q/20/TRK/200 

Purpose: 

To assign a FASTRAND-formatted mass storage file to this run. The 
file is used for plot dump output. 

Format: 

@ASG 

- Required mnemonics for the assign card. 


T 

- Assign card option denoting a temporary file assignment 


20. 

- Internal name assigned to the file. 


F40 

- Code specifying that this assign statement is for a 
FASTRAND-formatted mass storage file and identifying 
the specific device required. 


20 

- Integer specifying the minimum amount of mass storage 
required for this file. 
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Card 15 ; 

Image: 

Purpose: 

Format: 


Card 16 : 
Image: 
Purpose: 
Format: 


TRK - Increment or granule size for this assignment. 

* 

200 - Integer specifying the maximum amount of mass storage 

required for this file. 


0ASG.T 29. ,F4Q/1/TRK/10 

To assign a FASTRAND- formatted mass storage file to this run. The 

file is used for scratch input and output. 

@ASG - Required mnemonics for the assign card. 

T - Assign card option denoting a temporary file assignment. 

29. - Internal name assigned to the file. 

F40 - Code specifying that this assign statement is for a 

FASTRAND-formatted mass storage file and identifying 
the specific device required. 

1 - Integer specifying the minimum amount of mass .storage 

required for this file. 

TRK - Increment or granule size for this assignment. 

10 - Integer specifying the maximum amount of mass storage 

required for this file. 

@XQT S5. PREPRO 

To execute the SOAP Preprocessor program. 

@XQT - Required mnemonics for the execute control statement. 

S5. - Name of the file containing the absolute element PREPRO. 

PREPRO - Name of the absolute element to be- executed. 


2-13 








Card 17, 18 (A, B, C, D. E) t These cards represent the data images for the 
Preprocessor program. The specific contents of this deck are discussed in 
section 2.3. 

Card 19 : 

Image: 0ADD.P 4. 

Purpose: To add the Images on logical unit 4 to the runstream at this point. 

f » 

Format: @ADD - Required mnemonics for the add control statement. 

P - Add statement option indicating that the ADO statement 

appears in the listing. 

4. •* External file name of the file to be added to the run- 

stream. 

See section 2.3.1 for the options affecting the output listing of routine 
SELECT and the Memory Allocation Processor (MAP) for the SOAP. 

Card 20 : 

Image: @ED,U SW.MAPIN, SW. MAP INN 

Purpose: To call the edit processor to update the existing element MAPIN and 

to name the new updated element MAPINN. MAPINN will contain a set 
of collector directive cards as specified by card 21. 

Format: @ED - Required mnemonics to call edit processor. 

U - Edit card option indicating an update of the element. 

SW. - Internal file containing the element MAPIN. 

MAPIN - Element containing some col 1 actor directive cards 
standard to all runs. 

MAPINN - Updated element containg all necessary collector directive 
cards to produce the absolute element SOAP. 


Card 21:8 


Image; ADD SB. MAP/MTV 
ADD SB.MAP/MTVD 
ADD SB.MAP/MECD 
ADD SB. MAP/MEC2 
ADD SB.MAP/SSO 

Purpose: To add to the existing set of collector directives in element MAPIN 
an additional set of directives corresponding to the vehicle(s) 
participating in the simulation. 

Format: ADD - Edit command that adds all lines of the element 

specified to the element being updated. 

SB. - Internal filename containing the map element. 

MAP/VEHICLE - Element name and version of lines to be added to 

the updated element. 

MTV 
MTVD 
MECD 
MEC2 
SSO 


substitute for vehicle, no parentheses 


: Card 22 : 

Image: G>PREP SW. 

Purpose: To create an entry point table for file SW. 

Format: @PREP - Required mnemonics to call FURPUR processor. 

SW. - Internal filename. 


Card 23 : 

Image: @MAP,S SW.MAPINN,TPF$.SOAP 

Purpose: To use the set of collector directives in element MAPINN and create 
an absolute executable element, SOAP, and put it in file TPF$. 

Format: @MAP - Required mnemonics to call the MAP processor. 

S - MAP processor option requesting a short listing of storage 
allocation and Octal addresses of the relocatables included 
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SW, - Internal file containing MAP INN. 

MAPINN - Set of collector directive instructions. 

TPF$ - Temporary file to contain the absolute element SOAP. 

SOAP - The absolute element to be executed. 

Card 24 ; 

Image: 0XQT SOAP 

Purpose: To execute the SOAP. 

Format: @XQT - Required mnemonics for the execute control statement. 

SOAP - Name of the absolute element to be executed. 


Card 25 : This card represents the input data set for the SOAP. The composition 
of this data set is discussed in section 2.4.3. 

Card 26 : 

Image: @JSC*CALLUP.TAPELABEL 

Purpose; To automatically generate gummed labels and the corresponding tape 
library transaction for new save tapes. This card must immediately 
precede the tape assign card. 

Format: @JSC*CALLUP. - Complete file name of the CALLUP processor file. 
TAPELABEL - Processor for tape labeling and saving. 


Card 27 : 

Image: OASG.T $P-P$TAPE. ,U,,21 . FRM02/31PLT label. 

Purpose: To create a CalComp plot label and to create a save plot tape. 


Format: @ASG 

T 

$P- 


P$TAPE. 


- Required mnemonics for the assign card. 

- Assign card option denoting a temporary file assign- 
ment. 

- The required first three characters of the filename 
portion of the file specified if the tape is to be 
saved is an "X" bin tape requiring peripheral processing 

- Tape file name. If D option is used, this file name 
must be used. 
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U - Code specifying that this assign statement is for a 

magnetic tape device and identifying the specific 
type of unit required. 

21 - Retention period in days. 

FRM02/31PLT - Description field, which must limit entire field to 
24 characters and contain type of form to be used 
for the plot (02), number of plots, label for tape. 

Card 28 : 

Image: @USE 14. , $P-P$TAPE. 

Purpose: To link the internal file name (14.) to the external file name 

($P-P$TAPE) so that FORTRAN references to logical unit 14 are 
equivalent to system references to PLOT. 

Format: OUSE ’-Required mnemonics for the use control statement* 

14. - Internal or FORTRAN file name. 

$P-P$TAPE. - External or system file name. 

Card 29 : 

Image: ©REWIND 14. 

Purpose: To rewind the 

Format: OREWIND 

14. 

Card 30 : 

Image: OXQT S5, PLOTER 

Purpose: To execute the SSFS plotter program. 

Format: OXQT - Required mnemonics for the execute control statement 

S5. - Name of the file that contains the absolute element 

PLOTER. 

PLOTER - Name of the absolute element to be executed. 
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tape file PLOT, to the load point. 

- Required mnemonics for the rewind control statement. 

- Internal file name of the tape file to be rewound. 





Card 31 : This card represents the data set for the SSFS plotter program* The 
composition of this data set is discussed in section 2.5, 

Card 32 : 

Image: @FIN 

Purpose: To indicate the end of the runstream. 

Format: @FIN - Required mnemonics for the finish control statement. 
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2.3 PREPROCESSOR DATA 


Input to the Preprocessor consists of six data sets s a SELECT/MAP 
specification card, Environment configuration specification cards, flight 
software configuration specification cards, subroutine scheduling cards, 
print/plot output definition cards, and subroutine MAXMIN output definition 
cards. 

On all Preprocessor data cards the capability exists to add a comment to each 
card if the comment is preceded by the $ character. For this reason, care 
should be taken not to use TPF$ as a file name to load from. A $ character 
placed in card column 1 designates the whole card as a comment. 

The Preprocessor uses three mass storage work files: 3, 4, and 29. Logical 
unit 29 is the scratch file used to decode and build card images. Logical 
units 3 and 4 contain output pertinent to the main simulation. On logical 
unit 3, the Preprocessor writes the following card images: 

- OED.IN CALMOD 

- (Subroutine CALMOD code) 

- @ED, IN DUMMY S 

- (Subroutine DUMMYS code) 

- @F0R,N CALMOD 

- @F0R,N DUMMYS 

- @ED,IN SELECT 

- (Subroutine SELECT code) 

- @F0R,N SELECT 

2.3.1 SELECT/MAP 

This data card must be inserted only as the first card of the preprocessor data 
deck to affect the printing of the subroutines SELECT CALMOD, and DUMMYS, and 
the MAP of the SOAP, The format of this card begins in column 1 and is 
SELECT(N) ,MAP(N) 
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The control characters "N, S, L, X" are placed In parentheses as user specifi- 
cations. An "N" specification will cause the most abbreviated listing to be 
generated. The "N" specification Is the default for SELECT. An "S'' specifi- 
cation will cause a moderately comprehensive listing. An "L" specification will 
cause the most comprehensive listing to be generated. The "X" specification 
applies only to the MAP, and must be specified so that no MAP will be generated. 
Following this card, the Environment configuration specification data set must 
be included to ensure the proper generation of routines CALMOD and DUMMY. This 
data Is stored In their respective MAP/VEHICLE elements In the SB. file, 

2.3.2 ENVIRONMENT CONFIGURATION 

The Environment configuration specification data set (see table 2-1 I I ) is used 
to define the Environment model(s) and the file(s) where the model (s) reside. 

If no models are specified, the Preprocessor will generate the CALMOD subroutine 
which contains a call to the active vehicle (ACTVEH) driver, and the DUMMY S 
subroutine which contains entry points for all other models, In all cases where 
a dummy ACTVEH is requested, a dummy version of subroutine ACTVEH will be loaded 
from the SOAP System library. 

In the case where an Environment model is to be included in an overlay segment 
other than the "ROOT" segment, an Environment specification card must be included 
for that model and the version of the model requested must be specified as 
"0VRI.AY . " In the particular segment where the model is to be loaded, the 
model will be specified with the model's real version name. Also, an additional 
specification can be made to request DCOML, the routine used to define the 
length of all the COMMON blocks to be loaded from a fil e other than "W . " 

Each image is free- field and only one specification should be made per card. 

A comment can be added to a card if the comment is preceded by a $ character. 

A $ character in column 1 designates the whole card as a comment. 

All models desired must be specified for a simulation configuration. A dummy 
entry point will be generated for all models not specifically requested. 
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All names are limited to six characters each, 

The minimum input for this data set is an *E0F card. 

The general form of the model selection equation is as follows*. 

Image: model j =fi 1 e^. versionj 

Purpose: To define the Environment configuration specification for a 

particular mission phase, 

Format: Any number of free-field images consisting of a model selection 

equation. The order of the equations is irrelevant, 

models - Field used to identify the type of model that is being 
specified by this equation, (This mnemonic is defined 
in table 2-IV, ) 

file.j. . * Optional field which, If used, allows the user to 
specify the file where this model resides, "file^" 
is limited to six characters. If the field is omitted, 
file W is assumed, If only the period (.) is used, the 
default file is W. 

CAUTION: do not use the name TPF$. 

vers 1 on j *- Version name of the model to be used in this simulation. 

Dummy models may be selected by specifying the version 
to be DUMMY or version number 00. For a dummy specifics 
tion, an entry point is generated in subroutine DUMMYS- 
Available models are described in this document. 

The format of the *E0F card is: 

Image: *E0F 

Purpose: To denote the end of a data set within a data stream. 

Format: *E0F > Required mnemonic for the end of a data set. 




TABLE 2-IV.- ENVIRONMENT MODEL MNEMONICS AND DEFINITIONS 


MODEL * 
(model^) 

VERSION 

(version^) 

MODEL DEFINITION 

ACTVEH 

AVEH 

Active Vehicle 

IMU 

IMU 

Inertial Measurement Unit 

RCS 

RCS 

Reaction Control System 

TVC 

TVC 

Thrust Vector Control 

THRUST 

THR 

Engine Thrust 

ACS 

ACS 

Aerodynamic Control System 

ACCEL 

ACC 

Accelerometer 

. SNSCTL 

SNSC 

Sensor Control 

MASPRO 

MASS 

Mass Properties 

AERO 

AERO 

Aerodynamic Forces 

GRAV 

GRAV 

Gravity 

ATMOS 

ATM 

Atmosphere 

SLOSH 

SLSH 

Fuel slosh i 

BEND 

BEND 

Vehicle Bending ! 

GYRO 

GYRO 

Gyroscopic Sensors 

PASVEH 

PVEH 

Passive Vehicle 

GUST 

GUST 

Wind Turbulance 

HORIZS 

HORS 

Horizon Sensor 

STARTR 

STAR 

Star Tracker 

RENDES 

RSEN 

Rendezvous Sensor 

ORBALT 

OALT 

Orbital Altimeter 

NAVA ID 

NAVA 

Navigation Aids 

GRADAR 

GRAD 

Ground Radar 

RADALT 

RALT 

Radar Altimeter 

GGTORX 

GGT 

Gravity Gradient Torque 

ORBPAR 

ORBP 

Orbital Parameters 

AIRDAT 

AIRD 

Air Data Sensor 

WIND 

WIND 

Wind Effects 

RNGNAV 

RNGN 

Range Navigation 

LNDAID 

LNDA 

Landing Aids 

BARALT 

BALT 

Barometric Altimeter 

LNDGK 

LNDG 

Landing Gear 


♦Not all models available except as dummies. 
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2.3.3 SCHEDULING 

This Preprocessor data subset specifies the time-dependent flight control sub- 
routines and print/plot routines to be executed and defines their execution 
characteristics, such as frequency, priority, and start and stop times. 

As in the previous data sets* a comment can be added to a card if the comment 
is preceded by a $ character, A $ character in column 1 designates the whole 
card as a comment. Unlike the other data sets, all the comments associated 
with this data set are printed before the scheduling information table is 
printed. 

Each image is free-field and unique in that each image is processed during 
execution as a single request. Therefore, multiple scheduling intervals for 
one routine can overlap. 

All names are limited to six characters each. 

The data image is defined as follows; 

Image : NAME, PRIOR, DELTAT , STARTT , STOPT , NAMEB , NAMEA.TIMEl ,TIME2 

Purpose; To define the execution characteristics of the time-dependent flight 

software and print/plot routines. 

Format: Free-field images with comments. 

NAME - Name of the flight control subroutine, entry point, 
print/plot routine, or other schedulable routine. 

PRIOR - Relative priority for the execution of subroutine NAME. 

A value of 0.0 represents the highest priority. If two 
or more routines are scheduled for execution at the same 
time, the one having the highest priority will be executed 
first. 

DELTAT - Time interval, in seconds, between executions of routine 
NAME. If the field is left blank or input as zero, the 
value of DELTAT will be set to 99999.999. 
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START! - Simulation time, in seconds relative to simulator 

initialization, at which the flight software routine 
(NAME) is to start execution. A value of -1.0 in this 
'• field Indicates that the routine will be scheduled by 
another routine during the simulation, 

STOPT - Time, in seconds, at which NAME is to cease executing. 

If the field is blank or zero, STOPT will be set to lxlO 10 , 


NAMEB - Name of the subroutine that is to be executed just prior to 
the execution of NAME, This field is optional. 

NAMEA - Name of the subroutine that is to be executed Immediately 
after the execution of NAME. The field is optional, 

TIME1 * Start time for NAMEA and/or NAMEB. If omitted, the start 
time will be STARTT. If TIME1 is specified, TIME2 must 
also be specified . 

TIME2 - Stop time for NAMEA and/or NAMEB. If TIM El and TIME2 are 
omitted, the stop time will be STOPT. 

Like the other data sets, it must be terminated by an end-of-file card. 

Image: *E0f 

Purpose: To denote the end of a data set within a data stream, 

Formats *E0F - Required mnemonics for the end-of-data set statement. 

To ensure that the calling linkage to a routine or entry point is generated, 
a schedule data card must be present for each routine or entry point that has 
been scheduled or is to be scheduled for execution on the SOAP, A start time 
of -1.0 indicates that the routine will be scheduled by some flight control 
algorithm during the simulation. More than one schedule data card may be 
furnished for the same routine as long as each has a start time greater than 
or equal to 0.0. Minimum input for this subset is a schedule data card for 

. -m j, ... ■ . ' . . • ' . . ■' ■ 

subroutine TERMIN (the termination routine for SOAP) and the end-of-file card. 

The TERMIN card is required to prevent a "run-away" simulation. The Preprocessor 
checks for the presence of the TERMIN card and inhibits the simulation if it is 
not found. 
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When an Initial simulation time bias Is given on the simulation Input card INIT, 
all scheduled start and stop times will be altere'd by the Initial time value, 


A restriction which must be observed is that schedule data cards for Preprocessor- 

written routines cannot specify whether other routines are to be run as "before" | 

or "after" routines. Preprocessor-written routines include print, plot, and MAXMIN 

requests which are entry points in SELECT* However, the user can include these 

routines as "before" or "after" routines on schedule data cards of non-Preproces- 

sor written routines. Routines scheduled as "before" or "after" routines should 

not exit by calling FREQEJ. If they do, each time one of these routines is called, 

the primary routines will be rescheduled. Samples of routine scheduling cards J 

are shown in table 2-V. 1 

■ ■ ■ i 

‘ 

Default values are assigned for the various times and priority if they are not 


specified. 

The defaults are: 


PRIORITY 

default is 

0.0 

DELTA TIME 

default is 

99999.999 

START TIME 

default is 

0.0 

STOP TIME 

default is 

1.0E+1Q 


an exception is made for routine SIDUMP (state image dump). Whatever the 
scheduled or default priority associated with scheduling SIDUMP, it is changed 
to 1.0E+38 when SIDUMP is actually scheduled on the waitlist. 

Another SOAP system routine that needs special attention is SIPRIN (state image 
print). Like SIDUMP, SIPRIN can be scheduled or called. SIPRIN when executed 
will print all the SOAP standard COMMON blocks listed in table 2-VI. Entry 
points that can be scheduled or called are also listed in table 2-VI. All of 
these entry points are in SIPRIN. 

Note that one subroutine is scheduled automatically for each simulation. Sub- 
routine KALNOW sets the immediate Environment update flag, LECC0M(740), to force 
a call to the Environment. f ^ 
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TABLE 2- VI.- SOAP COMMON BLOCKS 


COMMON 

block 

Definition 

S/D 

Entry poinl 
in SIPRIN 

LECCOM 

Executive COMMON 

S 

LECCOP 

FSCOM 

Flight Control (FC) COMMON 

s 

FSCOMP 

ENVCOM 

Environment Executive COMMON 

s 

ENVCOP 

ACTVEC 

. ACTVEH COMMON 

S,D 

ACTVEP 

IMUC 

IMU COMMON 

S,D 

IMUP 

RCSC 

RCS COMMON 

S,D 

RCSP 

TVCC 

TVC COMMON 

S,D 

TVCP 

THRUSC 

THRUST COMMON 

S,D 

THRUSP 

ACSC 

ACS COMMON 

S,D 

ACSP 

ACCELC 

ACCEL COMMON 

S,D 

ACC EL P 

SNSCTC 

SCSCTL COMMMON 

S,D 

SNSCTP 

MASPRC 

MASPRO COMMON 

S,D 

MASPRP 

AEROC 

AERO COMMON 

S,D 

AEROP 

GRAVC 

• GRAV COMMON 

S,D 

GRAVP 

ATMOSC 

ATMOS COMMON 

S,D 

ATMOSP 

SLOSHC 

SLOSH COMMON 

S,D 

SLOSHP 

BENDC 

BEND COMMON 

S,D 

BENDP 

GYROC 

GYRO COMMON 

S,D 

GYROP 

PASVEC 

PASVEH COMMON 

S,D 

PASVEP 

GUSTC 

GUST COMMON 

S,D 

GUSTP 

HORIZC 

HORIZS COMMON 

S,D 

HOR1ZP 

STARTC 

STARTR COMMON 

S»D 

STARTP 

rendec 

RENDES COMMON 

S,D 

RENDEP 

ORBALC 

ORBALT COMMON 

S,D 

ORBALP 

NAVA I C 

NAVA ID COMMON 

S,D 

NAVA IP 

GRADAC 

GRADAR COMMON 

S,D 

GRADAP 

RADALC 

RADALT COMMON 

S,D 

RADALP 

GGTORC 

GGTORK COMMON 

S,D 

GGTORP 

ORBPAC 

ORBPAR COMMON 

S ,D 

ORBPAP 

AIRDAC 

AIRDAT COMMON 

S,D 

AIRDAP 


S = Single precision 
D * Double precision 



WINDC 

MIND COMMON 

RN6NAC 

RNGNAV COMMON 

LNDAIC 

LNDAID COMMON 

8ARALC 

BARALT COMMON 

FSC1 

Additional FC COMMON 

FSC2 

Additional FC COMMON 

FSC3 

Additional FC COMMON 

FSC4 

Additional FC COMMON 

FSC5 

Additional FC COMMON 

LNDGRC 

Landing Gear COMMON 

DATBUS 

Environment- co-flight 
control Bus COMMON 

DPDBC 

Environment- to-fl ight 
control Bus COMMON 

DPFSC ! 

FC COMMON 

DPFSClj 

Additional FC COMMON 

DAVEHCi 

ACTVEH COMMON 

DIMUG 

IMU COMMON 

DMASSC 

MASPRO COMMON 

DGRAVC. 

GRAV COMMON 

FS1C 

VEH#1 FC COMMON 

FS2C 

VEH#2 FC COMMON 

RCS1C 

VEH#1 RCS COMMON 

RCS2C 

VEH#2 RCS COMMON 

ACTV1C 

VEH#1 ACTVEH COMMON 

AGTV2C 

VEH#2 ACTVEH COMMON 

MASP1C ■ 

VEH#1 MASPRO COMMON 

MASP2C 

VEH#2 MASPRO COMMON 

H0RZ1C 

VEH#1 HORIZ COMMON 

H0RZ2C 

VEH#2 HORIZ COMMON 

AER1C 

VEH#1 AERO COMMON 

AER2C 

VEH#2 AERO COMMON 
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If there are no print/plot routine requests, the *E0F terminating the scheduling 
data set should be followed immediately by a card image with "QUIT" beginning 
In column 1. 

2.3.4 PRINT/PLOT 

This data subset allows the user to build print and/or plot dump subprograms 
which may either be scheduled or called. The user may specify up to 90 print 
routines and up to 10 plot dump routines. Any unique name (up to six charac- 
ters) may be used for the name of the print routine. Plot dump routine names 
are restricted to one of the following: PLOTR, PL0TR1, ...» PL0TR9. Variables 

to be printed and/or plotted must be selected from any location of any labeled 
COMMON block. The user specifies mnemonics (up to six characters) for each 
parameter to be printed or plotted. Note that no plotting actually takes 
place during the simulation, but that variables to be plotted are written on an 
output device (tape or drum) to be used by one of the SOAP postprocessors. 

* 

Print or plot dump routine data consists of two parts: a card defining the 

name of the routine and a card or cards indicating the variables to be printed 
or plotted. Each print or plot dump routine set must be terminated with an 
end of file (*E0F). The format of the card Images follows. 

Image: NAME 

Purpose: To define the name of the print or plot dump routine which will be 
built by the Preprocessor. The routine may be scheduled and/or 
called. 

Format: Fixed-field (A6) format. NAME must be left-justified in the field; 

i.e. , it must start in column 1. 

NAME - Name of the routine. A print routine may be named with 

any unique mnemonics (up to six characters). Plot routine 
names must be one of the following: PLOTR, PL0TR1 , . .., 

PL0TR9. If NAME is the word QUIT, it denotes the end of the 
Preprocessor output routine requests and the end of the 
Preprocessor input data set. 
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Image: VARNAM = GOMMON(I) OP t SF,FMT $ COMMENT 

Purpose: To define the parameters to be printed or plotted. 

Format: Free-field format. 

VARNAM - Mnemonics (up to six characters) used to identify the 

parameter in the printed output or used as plotter input. 

COMMON - Name of the labeled COMMON block in which the parameter 
is located. 

I - Location in the COMMON block where the output parameter 

is stored. This code is limited to the range from 1 to 
9999. 

OP - Optional operation code used to scale the value in COMMON(I), 

OP is one of the four arithmetic operators: or /. 

+ * 

- - Optional sign + or - can be assigned to the scale factor. 


SF - Optional field that allows the user to input a scale factor 
which will be applied to the data before it is output. Only 
the first 18 characters of SF are used. SF must be the 
same type as the data. In addition, 19 mnemonics are avail- 
able for SF, providing primarily English-metric conversion 
factors. The mnemonics and corresponding conversion factors 
are given in table 2- VI I . Also, the scale factor can be 
another COMMON name and location. In this case, no effort 
is made to decode the name and location. Therefore, the 
name must appear at least once in the "COMMON! I)" position 
in the image above. This appearance can be in any of the 
print/plot images specified. When using the COMMON blocks 
LECCOM, IMUC, MASPRC, NAVAIC, LNDAIC, and LNDGRC as a scale 
factor, the first letter must be replaced with a Q. This 
is true only when used as a scale factor. 

,FMT - Optional field that defines the format of the output param- 
eter. If FMT is specified, the comma (,) is required. FMT 
may be E or F for floating point, I for integer, or B for 
octal . If FMT is omitted, the standard FORTRAN convention 
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TABLE 2— VI I . — SCALE FACTOR MNEMONICS AND RESULTING VALUES 
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Image: 

Purpose: 

Format: 


$ 


Is used; that Is, floating-point data is assumed for all 
variables unless they begin with I, J, K, L, M, or N, in 
which case the integer format is assumed. 


Optional field that terminates the data scan on the card. 
The remainder of the card may be used for comments, The 
$ can be incorporated into card column 1 designating the 
whole card as a comment. The $ in column 1 must not be 
used before the NAME card, 


COMMENT - Optional comments field. The $ must appear prior to the 
comment. 


♦EOF 

To denote the end of file or the end of this print or plot routine 
definition. 

♦EOF - Required mnemonics for the end-of-file statement.. 


The print/plot routine definition subset is terminated by the name QUIT 
starting in column 1 and following the last *E0F. The QUIT card also repre 
sents the minimum data requirements for the subset, 


Examples of output routine definition cards are shown in table 2-VIII. 

When a Preprocessor-generated print routine is executed, the output will be 
printed with a title line consisting of the routine name and the current Simula 
tlon time. The variables requested will be printed in a NAME-VALUE sequence, 
five to a line, In the order and type specified during input. 

2.3.5 MAXMIN 

This data subset allows the user to build a routine called MAXMIN through input 
data cards using the same card format used to define the print/plot routines. 
MAXMIN Is a schedulable and/or callable routine designed to determine the 
maximum and minimum values for the variables specified. Upon completing 
the run and entering subroutine TERMIN, entry point MMPRNT in SELECT 


TABLE 2- VIII.- SAMPLE PRINT/PLOT ROUTINE DEFINITION CARDS 


Card column 1 


PLOTR 


♦EOF 


TIME 

ERRORX 

ERRORY 

ERRORZ 


LECCOM(l 

FSCOM(3Q 

FSCOM(31 

FSCOM(32 


♦RTD.E 

i*RTD,E 

i*RTD,E 


ACCPRT 


AVCAX = ACTVEC (85) 

$AV CONT ACC X 

AVCAY = ACTVEC 86) 

$AV CONT ACC Y 

AVCAZ = ACTVEC (87) 

$AV CONT ACC Z 

♦EOF 


PRTACC 



AVCAX = ACT V EC (85) • E 
AVCAY = ACTVEC 86), E 
AVCAZ = ACTVEC (87 ),E 

♦EOF 


MISC 

$ 

TIME * 
D 

123456 

JOB 

CTABLE 

I 

J 

A 

B 

L 

M 

Y 

BIG 

$ 

51 

52 

53 

54 

55 

56 

57 

58 

♦EOF 


SHOW MISCELLANEOUS CAPABILITY 
LECCOM(l) 

ENVCOM ( 324) ,1 

FSCOM( 8 ,B 

FSCOM( 2) ,B 

LECCOM( 4)*2.0,E 
LECCOM( 631 *10 ,1 
LECCOM ( 632 )*10 ,B 
ACTVEC (5) 

RCSC(5) 

MASPRC(5) ,E 

GRAVC(5),E 

GGT0RC(5) 

FSC1 (297)*. 23456789 
SHOW OPERATOR CAPABILITY 
IMUC(1)/-RTD 
RCSC ( 1 )*-PTK 
ACSC 1 /-16.13 
ACCC(1)*-14. 23763 
SL0SHC(20)+26.7 
BENDC (20)-0.0043 
FSCOM (20) * F SCOM (18) 

ACTVEC (20 )/ -ACSC (TO) 


■t 


* 


1 $ called to print the summary run The generation of this 

also be called or scheduled to print dur 9 terminated by 

routine Is Initiated when NAME Is MAXMIN. This data 

a *E0F . 


„ nonprated entry point MMPRNT is called, a summary print 
rjrSSTS of the print and then Hstln, the variable names 
LTe L, J and minima values and the times they occurred. 



2.4 SIMULATION DATA SET 


This section defines the SOAP simulation data card sets. 

2.4.1 SIMULATION TYPE 

Tables 2-IX and 2-X illustrate the two types of SOAP simulations. Initial and 
restart. The Initial simulation Is from cards only. 

Card 1 : 

Image: TYPE S.F 

Purpose: To define the SOAP simulation type. Do not add any type of comment 
to this card. 

Format: Free-field format. 

TYPE - Word INIT for an initial simulation and the word RESTAR 
for a restart simulation. 

S.F - When TYPE = INIT, S.F Is the initial simulation time 
in seconds. When this time is nonzero, the scheduling 
of all SOAP routines Is shifted by this value. 

When TYPE = RESTAR, S.F defines the particular state Image 
to be used for restart. S is the state image number and F 
is the particular file on the state image file In which it 
is located. 

2.4.2 RESTART SIMULATION 

The restart simulation initialization consists of reading a state image 
(labeled COMMON images) stored on magnetic tape or mass storage from a previous 
simulation and modifying it by additional card input. No routines are read. 

The state image input tape or file must be assigned to logical unit 7 for 
restart runs. Any new state images are output to logical unit 8. 
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TABLE 2- IX. - SAMPLE SOAP INITIAL INPUT DECK 


Column 1 


IN IT 40.0 

LECCOM (956) « 3000 $ FSCOM INPUT AND ... 

LECCOM (9 57) - 1000 $ ENVCOM INPUT AND ... 
ACTVEH (163) « 1 $ ATTFLG FOR IC PASS 

(Additional input data cards) 

♦EOF 

si l $ second trajectory (This Is a new simulation) 

FSCOM (384) * 1 

FSCOM (1485) * 0.0 

♦EOF 

NORMAL 

♦EOF 

TABLE 2-X. - SAMPLE SOAP RESTART DECK 


Card Column 1 

1 RESTAR 7.1 

2 LECCOM (621) = 19 $ PLOTR1 LU=19 

3 FSCOM (105) = PLOTR1 $ PLOTTER NAME 

(Additional input data cards) 

4 FSCOM (115) * 0.5 $ DELTA T FOR PLOTRl 


.. ... 



When setting up a restart simulation, the user must ensure that the COMMON 
assignments and lengths In all the routines to.be used In the restart are the 
same as those defined In the Initial simulation. If they are, the same routines, 
new routines, or edited routines can be used, 

Example: A sign error was found to exist In an atmospheric model equation for 
altitudes greater than h. If this model could be corrected without a COMMON 
block location change, a restart simulation could be made beginning at some 
altitude below h using the edited model. 

Therefore, the preprocessor configuration specifications may be changed for a 
restarted run. 

A complete set of preprocessor schedule cards must be defined for the restart 
simulation. This includes all before and after routines. The preprocessor 
uses these cards to generate the calling linkage in subroutine SELECT to all 
scheduled routines. In the typical case, the same set of schedule cards defined 
in the initial simulation Is used, However, start times, stop times, delta 
times, and priorities may be changed. A new print routine can be defined by 
the preprocessor and used as a "before" and/or "after" scheduled routine for 
debug purposes. 

If a simulation start time was specified on the initial simulation, then the 
user on the restart simulation must account for this time in the start and 
stop times of his restart schedule cards. ' 

The input data cards for the initial simulation are not used. Only the data 
needed to change the run is required. If the reason for this restart run 
is a program change, data input is not required. 

2.4.3 DATA INPUT 

The data input cards for an initial or a restart simulation have the same 
format (see tables 2- IX and 2-X). 


Image: COMMON (INDEX) * DATA W } SF 1 0P 2 SF 2 0P 3 SF 3 0P 4 SF 4 0P g SF g 

$ COMMENT 

Purpose: To scale and input data into SSFS single and double precision COMMON 

blocks. 

Format: Free-field format. , 

COMMON » Any one of the COMMON blocks defined in table 2-1 V. 

When defining COMMON block lengths in input cards, the 
specification should appear before any input reference 
is made to the particular COMMON block. The COMMON 
block must be dimensioned in a routine equal to or 
greater than the specification being made. 

INDEX - Location of the COMMON block to be initialized. 

DATA - Integer or an octal, floating-point, or Hollerith number 
to be stored in the COMMON block. An integer is signified 
by the absence of a decimal point; it may not exceed 
34,359,638,367. An octal number is signified by a B at 
the end of the number; it must be 12 digits or fewer. 

A floating-point number is signified by a decimal point; 
it is restricted to a card field of 24 columns or fewer 
and may be of normal or exponent form (e.g. , 100.0, 

1.0E+2, or 2.0D+3). A Hollerith number is signified by 
an initial alphabetic character; it must be six characters 
or fewer and is stored left-justified with blank fill* 

OP^ - Optional operation code used to scale DATA prior to loading. 
For each OP^ used (a maximum of five is allowed), an SF| 
must be supplied. OP^ is one of the four arithmetic 
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operators: +, *, or /. The evaluation of the input 

expression is in a strict left-to-right order, with all 
operators having the same weight. 

SFj - Input following an OP^ which is used to scale DATA prior 
to loading. SF must be of the same type as DATA. For 
numeric data, SF,. has the same restrictions imposed upon 
it that DATA does. SF,. may also be one of the 19 mnemonic 
conversion factors defined in table 2-VIL The same 
SF^ mnemonics are used for double precision. In this case, 
the system keys off of the COMMON block type. 

$ - Optional code used to terminate an input card, leaving 

the rest of the card available for comments. If $ is 
placed in column 1, the entire card is considered to be 
a comment card and is centered on the page. 

Card 5 : The input data set must be terminated by an *E0F card. 

Examples of data input cards are shown in table 2-XI. The number of inputs 
or input phrases per card is not restricted, but a data field may not be 
split between cards. Several data fields may be contained on one card (i .e. , 
for loading into consecutive COMMON block locations) with the data fields 
separated by commas. Data input may be continued from one card to the next 
with the comma at the end of the statement. For loading into several con- 
secutive locations, a particular location may be left unchanged by leaving a 
blank field on the card (i.e. , two consecutive commas in adjacent columns or 
separated by blanks). The data set is terminated with an *E0F card. When 
specifying the length of a COMMON block, the input value must be an integer. 

2.4.4 RECYCLE 

A recycle is initiated during a simulation run by scheduling subroutine RECYCL 
or by a call to RECYCL by the flight software. Initiation of RECYCL results 
in a pause in the simulation run to (1) read input cards and continue the 
simulation from that point or (2) "roll back" to a previous point in the sim- 
ulation run, read input cards, and continue the simulation from that point. 




2-40 


0 % 
• H 
in 


it £ 
25 


o 

o * 
o o 

H O 
o 


< H 

</v 


y 

o 



in 

O 

H 


• 

S5 •• 

> H 

<N 







H 

d 



p 







H 



H 






* 



© 

to 






• 

* 









ro 

VO 


21 

05 





cn 

H 

H 


O 

w 


0) 



w 




a 

© 


H 



01 

% 

«b 


0) 

w 


25 



p 

• 

in 


Cu 

EH 


W 



2 

d 

H 



£ 


SB 



H 

H 



* 

H 


O 



SB 


to 


m 



cu 




% 






X 




• 

H 


ii 



w 




H 




<0* 





t/> 

H 

% 










m 


cn 







* 

H 


00 


>4 

</K 




§ 



•w* 


>4 





© 

* 


SB 


04 


25 


o 

H 

Cl 


O 

in 

H 


H 


• 


H 


u 

» 

(H 


O 


Cl 

v 



CO 

* 

y 


CM 


1 

* 


• 

CM 






* 

o> 

H 

0V 


1 

Sc 


CD 


in 


H 

as 

% 



Cl 

25 

05 

• 

* 


as 

<N 

% 


1 

M 

« 

00 

• 

* 

a\ 


m 

25 

w 

►3 H 

CD 

1 

00 

O 

as 

II 

l 


H O 
O * 
0* © 
© 

O Os 
SB * 
H © 
H • 


Eh O 
O b 


if r> 
co 
I 

# • 
o vo 


H © 
o> 
* as 
as as 


© 

* 0> 
d cr> 


CM 

oo * 
v* N 
£ + 
O * 

a m 
ox i 

* .1 


• 

»— 4 

o 


© 



5 





H 




* 


vo 

a> 

H 

r»»» 

X 


Eh 

O 


H 

u 


05 



1 




© 



a> 


d 

1 

Eh 


o 

© 




W 



W 




• 

* 

% 


It 

oo 

CN 

© 

2 

H 

r-| 


z 


CD 


CO* 

© 




Cl 


in 

as 


w 


04 

£ 




» 


W 

1)1 


• 

p 

• 


« 





£ 

UJ 

25 


% 

«* 

w 



EH 

d 


© 

H 

© 

d 

in 


% 

* 

H 

o 

•pi 

H 

P 

o 

o 




25 

H 


CO 

© 

o> 

so 

• 

• 

•r 

0 

CO 

o 

£g 


o 

o 

H 


* 


H 




Ch 

© 

m 

r- 

m 


as 


cn 

P 

< 

>4 

o 



Cl 

© 


% 

© 


© 

as 

CO 

1 


% 

as 

s 

Cm 




H 

* 

SB 


■ « 


P 

• 

4h 

d 

as 

m 


«h 

n 

as 

o 



s 

0 


© 


It 

O 

co- 

H 

H 

H 

VO 

as 

vo 


t 


as 

u 



P 

z 

II 

r+ 



© 


H 


© 

in 

as 

a> 

o 

d 

«h 

as 

cn 




H 





© 



* 

4* 


as 

tn 

• 


(N 

as 

Cm 



P 

X 


. % ■ 

O 

VO 

H 


'* 


W 

ro 

as 

m 


<k 


as 




y 

H 

m 

© 


d 


ni 

P 

© 

© 

Cl 

as 


rH 

• 

% 

as 

% 



w 

?H 

m 

H 


av 

it 

« 

o 

... • 

• 

rH 

as 

n 

1 

H 

H 

as 

o 



H 

3 

ov 



«w 


Cl 

H 

VO 

d 













© 


£ 

© 

II 

II 

it 

ii 









Eh 

£ 

H 


O 

© 






^-s 



S"*« 



#— <* 

D 

o 



a 

© 





H 

CN 

cn 

•r 

CO 

H 

in 

o 

CU 

u 


cn 


H 

Cl 


vo 

d 

H 

H 

H 

H 

H 


d 

© 


u 



. 2* 

w 




Sw# 

V— » 


w 

w 

•>w« 

w 

Wf 

•w 


w 



01 

£ 

£ 

£ 

£ 

£ 

£ 

£ 

£ 

£ 

£ 

£ 

£ 

£ 


y 




O 

O 

O 

O 

o 

O 

O 

a 

O 

O 

O 

o 

o 




Cm 


O 

U 

o 

a 

a 

a 

a 

u 

a 

a 

a 

u 

u 






cn 

cn 

cn 

cn 

cn 

cn 

cn 

cn 

cn 

w 

cn 

p 

© 






cu 

Oh 

Cm 

Cm 

Cm 

Cm 

Cm 

CM 

Cm 

CM 

CM 

Cm 

Cm 


Hcim<*invod»o>©rHd 


in vo d cd 6 o 

Hf H H H H rH Cl 



TABLE 2-XI CONCLUDED 


o 


H O 

w • 

U cr> 

S> * 

o! o 

O * 


h o 


w 




U 

CN 















« 

% 











E 




04 

O 


% 









9 




CO 

« 


VD 



m 






H 

CO 



rfj 

H 


in 



• 






0 

U 



a 

II 


co- 



o 






a 

H 





o 











52; 



* 

O 

• 




EH 







O 



o 


o 




D 






M 

a 



1 


CN 

rn 



O 






03 

« 



H 


H 

4* 



cq 






U 

§ 




a 


CN 


CO* 

< 






w 

i § 



II 

o 

II 

+ 










j </> 




u 


H 



CO- 






52 





CO 

to* 

+ 









O 




H 

Cm 


CN 




Eh 

Eh 




H 





«n 


+ 


a 


CD 

D 








U 



H 


EH 


0< 

04 




H 




CO 

H 



o 

Q 


a 

S2 




52 




o 

* 

o 

II 

H 

* 

o 

H 

H 



1 

H 




a 

ro 

• 



a 

• 





Q 

Cm 





H 

ro 


CO- 

EH 

CN 

a 

a 



m 

W 




% 

V 

« 

ON 


Cm 

X 

EH 

Eh 



GO 

Q 




© 

(N 

© 



# 

a 

H 

H 



H 





• 

H 

• 

r> 


in 

w 

a 

a 



CN 

Q 




H 

* 

CN 

nn 


. • 

H 

w 

w 



rn 

2 






\ 

a 



X 










II 


o 

o 


V 

w 



CN 


• 

u 





hi 

• 

u 


GD 

a 

o 

o 

+ 


o 


s 

a 

X 

#•— * 

© 

CN 

CO 

CN 

• . 

H 

a 

a 

Q 


<h 

Eh 

EH 

Eh 

Eh 

H 

H 

1 

Pm 

*K 

r* 

X 



cn 



D 

H 

Cm 

CO 



o 


m 

H 

04 

CO- 

CO* 

• 


o 

04 




U 

% 

• 


\ 

* 

EH 





• 

2: 

4c 

* 

* 

OS 

00 

rH 

CN 

o 

CN 

CO 



rH 

a 

o 

H 




o 


4* 

* 

H 

CO 

X 



X 

Eh 

o 


O 

O 

O 

EH 

r* 

o 

CN 

1 


CO 




Cm 


rt 

• 

• 

• 

o 

■ % 

• 

«R 

in 

m 

EH 

Cm 


-f 

« 


Eh 

H 

H 

rH 

C5 

CD 

as 

H 

4- 

ID 

04 

W 


Q 

CN 

CN 

< 





•« 

* 

* 

CN 

r> 

X 

Q 


CN 

+ 

4- 

Q 




- *, 


O 

r- 

* 

o 

o 

U 



Q 

a 





o 

•h 

• 

* 

m 

• 

• 

a 

>1 

vD 

o 

o 

Q 





* 

© 

CN 

H 

in 

H 

rt 

X 

• 

• 

• 

►4 




H 

CN 








in 


in 

W 

11 

II 

II 



ii 

II 

II 

II 

II 

ii 

H 




H 




II 

II 








W 

II 

II 

Cm 

•*—* 




^ ^ 

<•> 

«—N 

r-N ■ 

*—v 

.*-> 

r-N. 

r* — » 





o 

H 

(N 


CN 


00 


o 

r> 

o> 

o 




« 

o 

o 

O 

H 

O 

CN 

*3* 

in 

GO 

ID 

GO 

r- 

o 

H 

CN 

W 

H 

H 

H 



r-> 

r> 

r* 


r> 

r> 

f^* 

H 

rH 

rH 

2 

k '-»* 



O 

•w* 

-w* 

■*W» 



N*4. 

w 

w 




CM 

a 

a 

a 


a 

a 

a 

a 

a 

a 

a 

a 

U 

u 

o 


o 

o 

o 

s 

o 

o 

o 

o 

o 

o 

o 

o 

CO 

CO 

CO 


u 

u 

u 

Eh 

u 

u 

u 

u 

u 

U 

o 

o 

&4 

Cm 

Cm 


i w 

co 

GO 

U 

C0 

CO 

CO 

CO 

CO 

CO 

CO 

C0 

04 

04 

04 


1 Cm 

Cm 

Cm 

< 

Cm 

Cm 

Cm 

Pm 

Cm 

Cm 

Cm 

Cm 

a 

a 

Q 


H 

CN 

00 


in 

GD 


00 

as 

o 

H 

CN 

ro 

^9* 

in 

CM 

CN 

CN 

CN 

CN 

CN 

CN 

CN 

CN 

CO 

CO 

ro 

ro 

ro 

co 






The SSFS has three types of recycle capabilities (see card 6, table 2-IX); 
NORMAL* SI, and SISAVE. 


a. The NORMAL recycle Instructs the SOAP to read input cards and resume the 
simulation run from that point, 


The SI recycle instructs the SOAP to read the state image specified on 
the SI input card. Input cards following the SI card are read, and the 
simulation run is continued from that point of the state image just read 
in. State image dumps that occur after the recycle are written over the 
old state images that follow the one used to reinitialize the simulation. 


The SISAVE recycle instructs the SOAP to dump a final (for this trajectory) 
state image on the tape, write an end-of-file code, and "roll back" to 
the state image specified on the SISAVE card. The state image is read 
into core; the tape is repositioned past the end-of-file; and the simulation 
reads additional input cards, dumps another state image, and resumes the 
simulation from the point of reinitialization. State image tapes produced 
by the SISAVE capability contain one file for the original trajectory, 
plus additional files for each SISAVE recycle. 


The SISAVE recycle capability was designed primarily to maintain a permanent 
record of a multiple trajectory run, so that any trajectory could be rerun to 
diagnose problems without having to repeat the complete run. The SI recycle 
is used more frequently for multiple trajectory runs, since most users are 
not concerned with saving a permanent record of the trajectories, and SI is 
more economical. The state image output file must be assigned to logical 
unit 8. 


For multiple trajectory runs, a NORMAL recycle deck (NORMAL, *E0F) must 
follow the last SI or SISAVE recycle deck to terminate the run normally. This 
is required because an SI or SISAVE recycle causes a backup in time and causes 
RECYCL to be rescheduled when that time is reached again. 


The optional recycle decks follow the simulation input deck shown in table 2-IX. 
All three types of recycles have the same input card structure or format. 
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Card 6: 


fl 


Image: TYPE S 

Purpose: To define the type of recycle. 

Format: Free-field format. 

TYPE - One of three words: NORMAL, SI, or SISAVE. 

S - Code used only when TYPE is SI or SISAVE. S is an integer 

corresponding to the state image at which reinitialization is 
to occur. 

Cards 7 and 8 : Cards 7 and 8 have the same format and purpose as input data 

cards 2, ...» 4 in section 2.4.3. 

Card 9: Each recycle data set is terminated with an *E0F card. 




2.5 POSTPROCESSING 


Postprocessing of the SOAP consists of plotting, using one of the two SSFS 
plotter packages: the production plotter or the high-speed plotter. The 

two plotters are denoted by their absolute element names, PLOTER and SPDPLT, 
respectively. With PLOTER, the user has the option of selecting CalComp, 
printer, and/or SD-4060 plots. With SPDPLT, the user can select CalComp and/or 
printer plots. 

PLOTER or SPDPLT may be run after a simulation, in an in-line mode, or in a 
subsequent run to produce plots of data dumped onto a tape or mass storage 
file during the simulation. 

2.5,1 PRODUCTION PLOTTER 

With the production plotter (PLOTER), the user has the option of specifying 
the following: 

a. One, two, or three variables on the same graph with a single abscissa 
parameter 

b. Logical unit and file from which the data is to be plotted 

c. Any combination of CalComp, SD-4060, or line-printer plots 

d. Frequency of the data to be plotted 

e. Number of points between symbols for CalComp plots for the set and/or 
for individual plots 

f. Graph and axis titles 

g. tabular listings of the parameters to be plotted 

h. Solid line or point plots 

i. Number of frames (limited to four) in which the data is to be plotted 
and also whether multiple frames are to have the same scale 

j. Standard CalComp scaling or an alternate SCALIE scaling routine, which 
may produce larger plots for the plot set 

k. Axis scaling for an individual plot, which overrides the selection in 
the previous option listed 


l. Plots from a portion of the input data file related to a specified time 
interval 

m. Multiple plot sets with minimum input requirements 

n. Choice of symbol representation for each variable 

o. Symbols may be rotated to display an additional variable 

The limitations in PLOTER include: 

a. Maximum of 2000 points per plot for CalComp plots. Additional points are 
processed as separate plots. 

b. Maximum of 200 points per frame for SD-4060 or printer plots. This is 
limited internally by recomputing the frequency of the points. 

c. Maximum of four frames per plot. 

2.5.2 HIGH-SPEED PLOTTER 

The high-speed plotter (SPDPLT) has the same capabilities and limitations as 
PLOTER with the following exceptions: 

a. No SD-4060 film plots 

b. No optional graph and axis titles 

c. One frame per plot request 

d. No more than the first 1350 points from the input file plotted for 
CalComp plots 

e. Limit of 100 prints internally for line printer plots 

f. Maximum of 30 plot variables 

g. No user-defined input scaling, although selection of SCALE or SCALIE 
for the plot set is permitted 

h. No capability to plot a time interval portion of the data file 



2.5.3 PLOTTER INPUT 

Despite all the options available In the plotters, the programs are designed 
to require a minimum of Input. In the simple plotter Input deck shown In 
table 2-XII, the plot request set control card contains only the plot Identifi- 
cation field, while the Individual plot request cards list only the variable 

names. This type of deck Is most often submitted and results In valid sets 
of printer and CalComp plots using standard options and nominal data values. 
For example, the blank plot request set control card implies that the data 
to be plotted is to be obtained from a mass storage file on file 1 of assigned 
logical unit 20. Every point is to be plotted to produce CalComp and line 
printer plots, using the SCALE routine for axis scaling. Using this infor- 
mation* plots are then to be processed for variables defined on the individual 
plot requests. 


TABLE 2-XII. - SAMPLE PLOTTER INPUT DECK 1 



Col . 


Col . 


Col. 


Col. 


Col. 

Card 

1- 


10-15 


20-25 


30-35 


m 

o 

1 

PLT01 









2 



TIME 


AIRSPD 





3 



TIME 


AVACX 


AVACY 


AVACZ 

4 



LAT 


LONG 








(Additiona 

.1 inpu 

,t ca 

irds) 



5 

*EOF 



_ 


u 





Additional input is required when nonstandard options are desired, as shown 
in the sample input deck in table 2-XIII. 

The input sets for the two plotters are compatible, with the high-speed plotter 
(SPDPLT) ignoring some of the options. The following sections describe the 
input cards for the production plotter (PLOTER); SPDPLT limitations are noted. 
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Image: IDPLT IFILE NPTS I FREQ NSET LU ISTORE I PART NOCAL NOPRT SD4060 ISCALE 

Purpose: To define the general input/output requirements for the following 

plot request set. 
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Format: 


Fixed-field. 
Card column 

Name 

• Format 

1-5 

IDPLT 

A5 

6-10 

IFILE 

15 

11-15 

NPTS 

15 

16-20 

IFREQ 

15 

25 

NSET 

A1 

29-30 

LU 

12 

31 

I STORE 

A1 

35 

I PART 

A1 

40-45 

N0CAL 

A6 

50-55 

NOPRT 

A6 

60-65 

SD4060 

A6 

80 

ISCALE 

A1 


IDPLT - Plot identification. If blank, no plot ID will appear 
on plots or printed output. 

IFILE - Defines the file from which the data is to be plotted for 
the following plot request set. If blank or 0, file 1 is 
assumed. 

NPTS - Is applicable only to CalComp plots. It specifies the 
number of data points between symbols for the set. If 
blank, it is assumed to be 50. This option may be over- 
ridden on individual plot requests. 

IFREQ - Defines the frequency at which the data for the following 
plot request set is to be plotted. If IFREQ is blank, 

0, or 1, all data is plotted; if 2, every other data point 
is plotted; if 3, every third data point is plotted; and 
so on. SPDPLT ignores this field. 


2-49 


NSET - If blank, individual plot request cards are required for 
all plots in the set, and this set of plot request cards 
becomes a standard set. If nonblank, the last standard set 
of plot variables is plotted as specified on this plot 
request card. Additional plot request cards may be 
supplied to supplement the standard set of plot requests 
for this set of plots only. 

NU - Specifies the unit that contains the data to be plotted 

for the following request set. If blank or 0, logical unit 
20 is assumed. 

ISTORE - If nonblank, the data to be plotted is obtained from tape. 

If blank, the data is obtained from the mass storage file. 

IPART - If nonblank, only the time interval portion of the plot 

data file is plotted, and a second set control card is 
required to specify the minimum and maximum plot times. 

This option is not valid in SPDPLT. 

NOCAL - If blank, a Cal Comp plot is produced for each plot in 

the succeeding plot request set. If nonblank, no Cal Comp 
plots are produced unless one of the plot requests in the 
succeeding set indicates otherwise. 

NOPRT - If blank, a line printer plot is produced for each plot 
in the succeeding plot request set. If nonblank, no line 
printer plots are produced unless one of the plot requests 
in the succeeding set indicates otherwise. 

SD4060 - If nonblank, an SD-4060 microfilm plot is produced for each 
plot in the following plot request set. If blank, no 
microfilm plots are produced unless one of the plot requests 
in the following set Indicates otherwise. This field is 
ignored by SPDPLT. 
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ISCALE - If blank, the standard Cal Comp SCALE routine Is used for 
the plots. If nonblank, the SCALIE routine Is used for 
axis scaling. SCALIE produces larger curves than SCALE, 
but the graph scale may not be as easy to Interpolate. 

This option may be overriden In PLOTER (but not In SPOPLT) 
by user-defined scaling for Individual plots. 

2,5, 3. 2 Optional Plot Set Control Card 

Card 2 (table 2-XIIIh This card Is supplied only when column 35 of card 1 

Is nonblank. 

Image: TMIN TMAX 

Purpose: To define a limited time Interval for the plot request set. 


Format: Fixed- field. 



Card column 

Name 

Format 

1-10 

TMIN 

F10.1 

11-20 

TMAX 

F10.1 


TMIN - Minimum plot time; may be greater than or equal to the 
Initial simulation time. 

TMAX - Maximum plot time; may be less than or equal to the 
simulation termination time. 

SPDPLT reads but ignores this card when column 35 on card 1 is nonblank. 

When this option is used In PLOTER, the following conditions must be met: 

a. TIME must be the first variable in the plot dump file. 

b. Logical unit 29 must be assigned as shown in table 2-II (card 15) 


2-51 



2. 5, 3. 3 Individual Plot Request Cards 


Card 3 (table 2-XIIII : 

Image: A I X I YA I YB I YC XAXS YAXS NSYM NF FS C P S 
Purpose: To define a plot to be produced. 

Format; Fixed-field, 


Card column 

Name 

Format 

2 

A 

A1 

9 

I 

A1 

10-15 

X 

A6 

19 

I 

A1 

20-25 

YA 

A6 

29 

I 

AT 

30-35 

YB 

A6 

39 

I 

A1 

40-45 

YC 

A6 

57-60 

XAXS 

F4.1 

62-65 

YAXS 

FA.l 

67-70 

NSYM 

14 

72-73 

NF 

12 

75-76 

FS 

A2 

78 

C 

A1 

79 

P 

A1 

80 

S 

AT 


A - For PLOTER, if A is an asterisk (*), then there Is another 
card Immediately following with format as specified In the 
description of card 7 following, which gives data determi- 
ning what symbol Is to represent a certain variable and 
at what angle the symbol Is to be printed. 
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I - For PLOTER, if I is an asterisk (*), a tabular 
listing of the data plotted ,is produced for the 
variable which follows the asterisk. For SPDPLT, if 
any I is an asterisk, a listing is produced for all 
variables specified on the plot request card. For both 
plotters, if the I preceding YA, YB, or YC is a slash 
(/), the corresponding plot (curve) consists of an X 
at each point instead of the normal line connecting 
the points. 

X - Name of the X-axis variable. 

YA - Name of the Y-axis variable. 

YB - Optional name of the second Y-axis variable. 

YC - Optional name of the third Y-axis variable. 

XAXS - Applicable only to Cal Comp plots. It defines the 

X-axis plot frame length in inches. If blank, 7 inches 
are assumed. 

YAXS - Applicable only to CalComp plots. It defines the 

Y-axis plot frame length in inches. If blank, 5 inches 
are assumed. YAXS is limited to 9 inches. 

NSYM - Applicable only to CalComp plots. This code specifies 
the number of data points between symbols. If blank, 
it is assumed to be NPTS, as specified in columns 11 
through 15 of card 1 of the data set. 

(NOTE: The remaining inputs are ignored by SPDPLT.) 

NF - Specifies the number of frames in wt^ch the data is 
to be plotted. If blank, it is assumed to be 1. NF 
is limited to 4. 

FS - Used only if NF is greater than 1. If nonblank, it 
specified that all frames have the same maximum and 
minimum Y-axis values. If blank, the scales for each 
frame are independent of the other frames. 
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C - If nonblank, a Cal Comp plot is produced for this plot 
in addition to the types requested on the plotter 
control card (card 1). If blank, the plotter control 
card specifies the types of plots to be produced for 
this request. 

P - If nonblank, a printer plot is produced for this plot 
in addition to the types requested on the plotter 
control card (card 1). If blank, the plotter control 
card specifies the types of plots to be produced for 
this request. 

S - If nonblank, an SD-4060 plot is produced for this plot 

in addition to the types requested on the plotter control 
card (card 1). If blank, the plotter control card 
specifies the types of plots to be produced for this request. 


Card 7 (table 2-XIIIl: 

Image: ISYAANG ISYBANG ISYCANG 

Purpose: To allow user specification of symbol to represent any of three 
Y variables (see below for choice of symbols) 

Format: Fixed-field 


Card column 

Name 

Format 

9 

IS 

11 

10-15 

YAANG 

A6 

19 

IS 

11 

20-25 

YBANG ' 

A6 

29 

IS 

11 

30-35 

YCANG 

A6 


IS - Integer defining choice of symbol 

YAANG - Name of variable angle at which the first Y-axis variable symbol 
is drawn. 
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YBAN6 - Same as YAANG for second Y-axis symbol. 
YCANG - Same as YAANG for third Y-axis symbol, 


2. 5. 3.4 Optional Plot Request Cards 

Any combination of the optional cards defined in this section may be inserted 
immediately following an Individual plot request in the PLOTER input deck. 

Axis scaling card ; 

Image: S XMIN XMAX YMIN YMAX 

or 

S' XMIN DELX YMIN DELY 

Purpose: To permit user input of axis scaling parameters. This card is ignored 


by 5PDPLT. 



Format: Fixed-field. 



Card column 

Name 

Format 

1 

S 

A1 

2-12 

XMIN 

Ell . 1 

13-24 

XMAX 

El 2.1 


or 


25-36 

DELX 

El 2.1 
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Card column 

Name 


Format 

37-48 

YMIN 

* 

E12.1 

49-60 

YMAX 

or 


El 2.1 

61-72 

DELY 


E12.1 


S - Character S, which must appear in column 1 of the 
axis scaling card 

XMIN - X-axis minimum value 

XMAX - X-axis maximum value 

or 

DELX - X-axis units per inch of plot length 
YMIN - Y-axis minimum value 

YMAX - Y-axis maximum value 

or 

DELY - Y-axis units per inch of plot height 

When an S appears in column 1, minimum values of X and Y must be specified. 
Either maximum values of X and Y or unit increments per inch on the X and Y 
axes must be specified consistent with the dimensions specified in XAXS and 
YAXS of the plot request card. 

Plot title card : 

Image: G LONGITUDE VS. LATITUDE 

Purpose: To permit the specification of a title for a plot. This title 

card is applicable to the plot request card immediately preceding 
this card (i.e., card 4 in table 2-XII. This card is optional for 
PLOTER and is ignored by SPDPLT. 
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Format: Fixed-field, 



Card column 

Name 

Format 

1 

G 

A1 

2-80 

TITLE 

13A6.A1 


6 - Character 6, which must appear In column 1 of the plot 

title card. 

TITLE - Title of the plot. 

Axis title cards : 

Image: A LONGITUDE (DEGREES) LATITUDE (DEGREES) 

Purpose: To permit the specification of a title for the X and Y axes of the 
plot. Two axis title cards may be used to define an X-axis and up 
to three Y-axes. The axis title card is applicable to the plot request 
card immediately preceding it. The X-axis title is defined on the 
first half of the axis title card. The first Y-axis title is defined 
on the last half of the card. If two or three Y-axis variables are 
defined, a second axis title card is permitted, with the first half 
of the card defining the second Y-axis parameter and the second half 
defining the third Y-axis parameter. This card is optional for 
PLOTER and is ignored by SPDPLT. 


Format: Fixed-field. 



Card column 

Name 

Format 

1 

A 

A1 

2-37 

AXIS! 

6A6 

41-76 

AXIS2 

6A6 


A - Character A, which must appear in column 1 of the axis 
title card. 

AXIS1 - Title for the X-axis or the second Y-axis. 

AXIS2 - Title for the first or third Y-axis. 
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End-of-file card: 


Image: *E0F 

Purpose: To denote the end of file or the end of a data set within a data 
stream. 

Format; *E0F - Required mnemonics for the end-of-file card. 

Any number of plot request sets may be inserted in the input deck, with an 
end-of-file card between each set. Two end-of-file cards should follow the 
last data set. 

2.5. 3. 5 Logical Unit Assignments 

For most plotter runs, logical units 3, 14, 20, and 29 must be assigned by 
inserting corresponding @ASG and @USE statements in the input deck. Table 2-1 
describes the input/output logical unit requirements for the plotters. 

2.5.3. 6 Diagnostics 

If one or more of the dependent variables requested as YA, YB, and YC are not 
defined in the plot .dump routine (section 2.5), the undefined variables 
cannot be plotted, but a message is printed at the bottom of the plot frame. 
For example, 

PP IS NOT DEFINED PLOT IGNORED 

If the independent variable is undefined in the plot dump routine, then a 
blank frame is printed with the following message at the bottom: 

INDEPENDENT VARIABLE IS NOT DEFINED PLOT SKIPPED 

2.5. 3. 7 Plot Output Data Format 

There is a requirement for the SOAP to generate a data file or tape to be 
used as input to the SOAP plot processors and/or other programs. Through 
the SOAP Preprocessor scheduling and plot routine generation capability, it 
is possible to output any information stored in labeled COMMON blocks at the 
scheduled plot delta time. 
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The SOAP plot output data is written using the standard computer system 
format onto a file or tape in the following manner: 

a. Symbol record 

b. Data records 

c. EOF records 

The symbol record consists of the variable names as assigned in the SOAP 
Preprocessor plot routine data set and the number of plot variables (N). 

The symbol record is 256 words long and is arranged in the following manner: 

a. Locations 1 through N (where N < 254) - variable names assigned to plot 
data values. 

b. Location N + 1 - number of variables (N) in the plot data set. 

c. Location 256 - Number of variables (N) in the plot data set. 

The data records consist of the values for the plot data set. These records 
are N words long and are written out each time the plot routine is executed. 

The EOF record is 256 words long and is used as a file mark. It consists of 
a coded end-of-file (left-justified Hollerith data - EOFbbb) in the first 
word. The values in the other 255 words are not important. The SSFS has 
multifile capability, and each file is terminated by an EOF record. Two EOF 
records are written at the end of the last file. It is necessary that this 
last record be a maximum of 256 words long because of the variable lengths 
possible for the Data records. 

Table 2-XIV is a sample format of the plot output data and table 2- XV is 

the coding necessary to read the plot output data if used as input to another 

program. 

It should be noted that if an SSFS simulation is made with the plot output 
data going to a file and later copied to tape, the tape cannot be read 
directly. It would be necessary to copy the tape back to a file, which can 
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TABLE 2-XIV. - PLOT FILE FORMAT 


TIME^ 

RX^^ 



0.0 

6.5E+6 

o 

• 

o 

i 

O 

• 

o 

BD 

6.4E+6 

1.0E+5 

-l.OE+5 






490.0 -4 . 1E+5 3.1E+6 -3.1E+6 

500.0 -4.2E+5 3.2E+6 -3.2E+6 

EOFMMjtf 


TABLE 2- XV. - CODE TO READ PLOT FILE 


DIMENSION BUFF (256), IBUFF(256) 
EQUIVALENCE (BU.?F(1) , XBUFF(l)) 
DATA IEOF/3HEOF/LU/20/ 


! READ SYMBOL RECORD 

READ (LU) IBUFF 
NDP ■ IBUFF (256) 

! * READ DATA RECORD 

,100 READ (LU) (BUFF (I), 1-1, NDP) 

! IS THIS EOF RECORD 

IF (IBUFF (1) .EQ.IEOF) GO TO 200 

! *** USE DATA *** 

GO TO 100 

: • 

*** ALL DATA READ *** 

200 CONTINUE 
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then be read. If the data goes directly to tape during a simulation, the 
data can be read directly from tape. 


Also, any tape or file generated by other programs in the above format could 
use the SOAP plot capabilities. 


i 

i 





3. UTILITY ROUTINES 


These routines comprise SOAP preprocessor, executive, and postprocessor rou- 
ines, and general purpose routines which the user may include In a flight 
control model. Typically the user need not be at all concerned with the 
systems features and need only select those items from sections 4 and 5 
which are desired. The only caution which must be included in this recommen- 
dation Is that some mathematical functions have both single precision and 
double-precision versions. The double-precision versions support double- 
precision environment routines, while the single-precision versions may be 
used in the same simulation by flight control routines. The user need only 
be aware of the distinction and should not attempt to substitute. 

3.1 PREPROCESSOR ROUTINES 

The preprocessor writes the routines SELECT, CALMOD, and DUMMYS based on the 
environment input cards. It creates the appropriate wait list based on 
schedule cards and on print/plot requests. Typically, the user will neither 
call nor modify these routines directly, but will simply use the Preprocessor 
in accord with the runstream instructions in section 2. 

Preprocessor routines are: 

ENVRD 

KARDKR 

PREENV 

PRMAIN 

PROUT 

A functional diagram is given in figure 3-1 and a discussion of each routine 
follows. 

ENVRD 

Subroutine ENVRD first checks to see if the optional input card for controlling 
the listings of the FORTRAN subroutines (CALMOD, DUMMYS, SELECT) generated 
during preprocessor execution. If it is, the appropriate flags are set to 
obtain the desired listings. By default the most abbreviated listing of the 
preprocessor-generated FORTRAN subroutines are given. Subroutine ENVRD then 
reads and decodes a set of free- field input cards specifying the Environment 
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PRMAIN - Preprocessor MAIN program. 

ENVRD - Reads and decodes the Environment configuration 
selection cards. 

PREENV - Matches Environment model selection with desired 
configuration. The FORTRAN subroutines CALMOD and 
DUMMYS are generated on logical unit 3 (LU3). 

r ™. PREENV writes compilation and editing instruc- 
tions for CALMOD, etc., on LU4. 

PROUT - Builds FORTRAN subroutine SELECT from the schedule and 
print/plot data sets, and writes SELECT on LU3. 

KARDKR - Decodes print/plot data cards for subroutine PROUT. 


Figure 3-1 . - Preprocessor functi onal di agram . 




rSf • • . 
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models to be used. The cards consist of phrases of the form: 

MODEL - FILE. VERSION 

where 

MODEL is the name of an Environment function (e.g., ACTVEH). The element 
specified is typically the location of relocatable code for a driver 
program, 

FILE is the name of the program file 

VERSION is the version of the particular Environment model to be used to 
perform that function (e.g. , AVEH1 ) 

One input phrase to ENVRD is unique. This phrase selects the COMMON block 
definition routine (DCOM = ). In processing the input, ENVRD stores the 
desired Environment configuration specified into two arrays with the object 
code names and versions (the object code name is of the form "NAME/VERSION") 
of the particular Environment model subroutines specified. After processing 
the input cards, ENVRD returns control to PRMAIN. 

KARDKR 

Subroutine KARDKR is called by PROUT to process each output routine variable 
definition card. The free-field cards (one per variable) specify a COMMON 
block name, a location within the COMMON block, a Hollerith name for the 
variable and an optional format and scale factor. KARDKR decodes the Hollerith 
card image and returns the information to PROUT. This includes identifying 
mnemonic scale factors, if specified. Whenever a bad card is detected, KARDKR 
sets a flag for PROUT. 

PREENV 

PREENV processes the input configuration cards decoded by ENVRD specifying 
the particular Environment model subroutines which are to be used during the 
simulation. During a simulation job, a secured file containing Environment 
model subroutines is assigned to the run. This file (designated by the single 
Tetter W) may contain several different subroutines performing the same 
Environment function and therefore being referenced by the same name. PREENV 
makes it possible for the computer system allocator to choose the correct 
Environment model subroutine to load into core* The environment model specified 
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calls any number of differently named subroutines, one for each vehicle, 

PREENV begins writing FORTRAN editing and compilation directives on LU4, 

These directives will compile the FORTRAN routines (CALMOD, DUMMYS, and 
SELECT) generated on LU3, The routines are first loaded using the ED proces- 
sor and then compiled. This sequence allows the SOAP to be executed inter- 
active from the terminal or in a batch mode. PREENV then uses the information 
supplied by ENVRD to begin writing an allocation element onto LU4. This 
portion of the element will supply the allocator with the necessary Instructions 
to call the desired Environment model subroutines and their associated Block 
Data subprograms Into core for the simulation. PREENV then writes the FORTRAN 
code onto LU3 for subroutines CALMOD and DUMMYS. PREENV then returns control 
to PRMAIN. 

PRMAIN 

PRMAIN is the main program which controls the flow of Preprocessor execution. 
Initially, it calls subroutine ENVRD for the selection of the particular 
Environment model subroutines which are to be used during the simulation, 

It then calls PREENV to define the Environment load directives from the input 
decoded by ENVRD and to build two routines for Environment interface. After 
returning, PRMAIN calls PROUT, which calls KARDKR, to build a routine for 
flight control routine scheduling „nd printing and/or plotting of requested 
simulation parameters. After calling PROUT, PRMAIN prints the duration 
(in seconds) of preprocessor execution and terminates. 

PROUT 

PROUT processes three sets of Preprocessor input cards, 
a. Schedule data cards specify the names of the flight software routines 
and output routines to be executed during the simulation. Accompanying 
each routine name is a priority, delta time, start time, and stop time. 
There may also be Specified the name of a routine to be executed before 
or after the scheduled routines call , with start and stop times. Using 


the schedule data cards, PROUT bu11ds,a table of data to be used later 
as Input (via data statements) to the SOAP Executive, 

b, The print and plot output routine data cards provide the names of output 
routines which the preprocessor Is to build and define the variables to 
be output during the execution of each routine. 

c, The maximum/minimum routine data cards provide the name of the output 
routine (MAXMIN) which the Preprocessor Is to build and define the vari- 
ables to be used during the execution of the routine. An entry point is 
generated when using this routine which Is called at the end of the 
simulation to print a summary of maximum/mi nlmum values of the variables 
specified. 

As output, PROUT writes the FORTRAN subroutine SELECT onto LU3. PROIJT begins 
by reading the schedule data cards and storing the information. In the pro- 
cess, PROUT replaces all stop times and delta times which are less than or 
equal to zero with values of 10^° and 10 3 ^, respectively, to prevent looping 
during the simulation* 

After finishing (reading an end-of-file (EOF) card) the schedule cards, PROUT 
processes the output routine data cards. These cards occur In groups, each 
of which contains an output routine name card, followed by one or more 
variable definition cards, and ending with the EOF card. PROUT uses the 
input cards to generate several arrays which later will be used in writing 
portions of SELECT. The processing of each variable definition card is 
accomplished by a call to KARDKR. The end of the output routine data cards 
is signified by the word "QUIT" in columns one through four of an input card. 
After this card has been read, PROUT begins the Writing of subroutine SELECT 
onto LU3. SELECT contains branching logic used to transfer control to the 
first routine on the wait list during simulator execution and all output 
routines specified by the Preprocessor input, each routine having its own 
entry point within SELECT. If one of the words PLOTR, PL0TR1, PL0TR2, ..., 
PL0TR8, or PL0TR9 appear as an output routine name (signifying a plot dump 
routine), a call to PLTOUT is inserted in the SELECT logic for plot dumps. 

If the output routine name specified is MAXMIN, the maximum/minimum logic 
is inserted in the SELECT logic. For all print routine names, a call to a 
DATOUT is included in the output routine logic for the actual printing. 

After writing SELECT, PROUT returns control to PRMAIN. 
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3,2 SOAP EXECUTIVE SYSTEM 


Figures 3-2 y 3-3, and 3-4 depict SOAP operation, Commands are carried out, 
flight control routines are executed, and environmental Information Is up- 
dated when the time designated for such an operation and the simulation time 
coincide. The wait list and command file hold lists of operations to be 
performed and commands to be executed, respectively. As each event takes 
place, It Is removed from the appropriate list. In order to cause an entry 
to be made In the command file, CMDFIL Is called by the user. To cause a 
routine to be scheduled or rescheduled, calls to FREQEJ (or FREQUP), CHGBLD, 
REWTLT , or SCHED are used. 

The user will not have a need to call most of these routines directly. The 
exceptions include those routines mentioned above, ENDJOB (to remove a routine 
from the wait list), and 'routines such as TERMIN, SIPRINT, and SIDUMP. 

TERMIN ends simulation' execution and outputs lists of all COMMON block locations, 
SIPRINT causes a state image, i.e., COMMON block location printout. SIDUMP 
is similar to SIPRINT but a variety of output devices are available. 

When constructing flight control routines, the user must remember that any 
routine which is to be scheduled for execution at regular intervals must be 
Identified to the system by a schedule sard (see section 2). After execution 
of such a routine, it must be rescheduled by a call to FREQEJ or FREQUP 
within the subroutine . Routines which are only called by other routines have 
no special requirements. 

The Input/ output routines INPUT and PINOUT have versions designed to deal 
with c'ouble-preclslon blocks. The names are the same, the difference is 
specified in the mapping stage of execution. Listed below are the names of 
subroutines Involved in the SOAP executive; descriptions follow. 


BLKPLT 

DUMMYS 

FREQEJ 

!$PPRT 

PLTOUT 

RESTRT 

S0RT3V 

CALENV 

ENVC 

FREQUP 

JSPPRT 

PNNOUT 

REWTLT 

SREAD 

CALMOD 

ENVCAL 

HDLNS 

KALNOW 

PRTSDT 

SANDF 

SUPTIM 

CHGBLD 

ENDJOB 

HEDLNP 

MAIN 

PRWAIT 

SCHED 

SWRITE 

CMDEXE 

ENVSKL 

INOUT 

NEXTIM 

QZKILL 

SELECT 

SYSPRT 

CMDFIL 

EXECT 

INPUT 

NINOUT 

RECOV 

SIDUMP 

TERMIN 

DPOUT 

EXECTR 

INPUT2 

PINOUT 

RECYCL 

SIPRIN 

TIMER 


1 



8 ( 

H : l 


-I 


Figure 3-2.- Relationship of SOAP operational blocks. 
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Figure 3-3.- Executive functional diagram. 

(subroutine names are underlined) 

































BLKPLT 

Subroutine BLKPLT (NPOINT, NDATA , I BCD, ARRAY) 

BLKPLT may be called by the flight software to dump one file of data for 
subsequent plotting. The data to be plotted must be stored in ARRAY* 
dimensioned ARRAY (NPOINT, NDATA). ARRAY thus contains NPOINT data points 
for NDATA variables. IBCD is an array of length NDATA containing BCD 
(e.g., TIME) names corresponding to the variables to be plotted. These 
names will be used by the plot processor to identify requested plots. After 
dumping the requested data, BLKPLT writes an EOF. BLKPLT uses LU20 for output. 


CAI.ENV 

CALENV is used by all areas of the simulator for updating the Environment. 
Scheduling or calling CALENV ‘is the only way in which the Environment can 
properly be updated. If an update is needed but not immediately necessary, 
LECC0M( 704 ) may be set to a 1 (integer). The Executive will then update the 
Environment (via CALENV) before the next flight control Routine is executed. / 

If an immediate update is required, the flight control should call CALENV 
directly. In either case, LECC0M(869) should be set to a 1 (floating-point) 
if a full Environment update is desired. Normally, the Environment update 
consists of only the active vehicle dynamics and the dynamics sensors. With 
LECC0M(869) set to a 1, all of the Environment is updated, including the non- 
dynamic sensors, ■ 

If CALENV is entered and the Environment has previously been updated to the 
current time and there are no commands to be incorporated, control is returned 
to the calling program. . 

CALMOD 

CALMOD makes direct calls to the Environment model routines. 


CHGBLD 

Subroutine CHGBLD (JOBNAM, PRIOR, DELTAT, TSTOP) 

CHGBLD may be called by the flight software to alter the Preprocessor data 
sets for scheduled routines. In this manner, the flight control can define 
its own future scheduling so that subsequent calls to FREQUP, FREQEJ, SANDF, 
or WLINFO will use the new data. JOBNAM (six-character or less binary code 
decimal (BCD) name) identifies the routine whose data set is to be altered. 

CMDEXE 

Subroutine CMDEXE handles the command execution function. Its operation is 
shown functionally in figure 3-5. Actually, there is logic in both ENVSKL 
and CMDEXE (not shown in functional flowcharts) for separate calls to CMDEXE 
for Dynamics commands and for Sensor commands. 

CMDFIL 

Subroutine CMDFIL (CODE, DATA) 

CMDFIL is called by the flight control to file away commands for later execu- ' ' 
tion by the Environment. The commands represent the flight computer output 
interface with the flight hardware. For the SOAP, 1 the commands may range from 
a simulation of output channel bit setting (the most detailed simulation) up 
to the most function level (e.g., a thrust pointing vector). 

CMDFIL is called with two arguments, CODE and DATA. CODE is the actual command 
code recognized by the Environment. DATA is a one-, two-, three- piece set of 
data representing the information associated with CODE. DATA might be a 1 or 
a 0 for a bit on or off, or it might represent a pointing vector. Both CODE and 

DATA are defined for each Environment model. 

i 

When called, CMDFIL first checks to ensure that it has room to file the edmmand. 
At the present time, a maximum of 49 commands may be accumulated. If the table 
is full, CALENV is called to update the Environment and to incorporate any com- 
mands on the list scheduled to be executed before the current flight software 
time. This will make room on the push list for the new cor<mand. Before the 
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DETERMINE WHICH MODEL 
RESPONDS TO THE NEXT 
COMMAND 


COUNT NUMBER OF COMMANDS 
FOR THIS MODEL AT THIS 
TIME AND STORE IN 
NCMMDS 


CALL MODEL (VIA CALMOD) 
TO INCORPORATE COMMANDS 


REMOVE NCMMDS NUMBER 
OF COMMANDS FROM LIST 
AND PUSH UP 


NEXT COMMAND AT 
's s THIS TIME^x 


NCMMDS = 0 


RETURN 


Figure 3-5-- Command execution functional diagram. 





command is filed on the push list, any time delays associated with the command 
are added to determine the command execution time. The time delay is defined 
by the Environment, along with a flagword. The fjagword specifies whether the 
command invalidates further Environment sensor extrapolation and whether an im- 
mediate Environment update' is required by the issuance of the command* 


If an immediate update is flagged, it is to be forced at the command execution 
time. The routine KALNOW, which sets the update flag, LECC0M(704), tol is 
scheduled at the highest priority at the command execution time. The 
correspondence of command codes with models used in the SOAP can be found in 
table 1-1. 

DPOUT 

Subroutine DPOUT (DPA, NF, NL, NDPCB) 

DPOUT is a double-precision print output routine used to print the double- 
precision named NDPCB array DPA from location NF to NL irt the double preci- 
sion common block named NDPCB. 


DUMMYS 

DUMMYS supplies entry points with names corresponding to all Environment models 
which are not specifically required for the simulation at hand. 


ENVC 

I 

! Subroutine ENVC is an entry point in ENVCAL and performs identically to ENVCAL. 
ENVCAL 

Subroutine ENVCAL is a schedulable routine which when called performs an im- */ 
mediate environment update by calling subroutine CALENV. 


ENDJOB 

Subroutine ENDJOB is an entry point in SELECT and is used to terminate the 
execution of flight software and to return simulation control to CONTRL. 

This special termination routine is used to ensure that CONTRL regains control 
after the flight software execution, since the normal FORTRAN subroutine 
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linkages (using the RETURN statement) are sometimes destroyed during a simula- 
tion. All flight control routines which are scheduled or which complete the 
execution of a scheduled routine should call ENDJ.OB or FREQEJ, which calls ENDJOB, 

ENVSKL 

Subroutine ENVSKL is the Environment Executive routine and is the central 
control routine, Figure 3-4 represents a functional flow of ENVSKL, At 
the start of a simulation run (initialization pass), ENVSKL calls each of 
the functions which may be modeled. If a particular function is not modeled 
in a run, its call will be answered by a dummy routine supplied by the Environment 

system. 

On subsequent calls, ENVSKL updates the Dynamic models, and if the full update 
flag, REASON, is set to 1, it also updates the Sensor models. The update con- 
tinues, stopping to incorporate scheduled commands, until Environment time is 
equal to the flight control time, 

EXECT. 

EXECT is called by the MAIN program of the SOAP and accomplishes the simulator 
initialization. Initialization may involve a complete simulator initialization 
or a restart run (i.e., a resumption of an earlier run). EXECT calls INPUT 
to read the SOAP input cards, CALENV to initialize the Environment, SIDUMP 
for the optional initial State Image dumps, and SIPRIN for an initial State 
Image print for diagnostic purposes. After initialization, EXECT returns 
control to MAIN, which calls CONTRL for the actual simulation. EXECT should 
never be called by the flight control routines. 

EXECTR 

EXECTR is an entry point in CONTRL and is called first by MAIN to start the 
simulation. EXECTR defines the first job and the time of the first job, and 
calls SELECT to execute the job. 
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FREQEJ 


Subroutine FREQEJ Is an entry point In FREQUP and Is used by flight control 
routines for termination and rescheduling, FREQEJ should be called to complete 
the execution of any scheduled flight control routine or any nonscheduled 
routine which completes the execution of a scheduled one. FREQEJ reschedules 
the currently scheduled (and executing) waitlist routine using information 
from the appropriate Preprocessor data set. After the routine has been 
rescheduled, FREQEJ calls ENDJOB to return control to the simulator control 
loop. 

FREQUP 

Subroutine FREQUP performs the same function as FREQEJ, except that control 
is returned to the calling routine, For both FREQEJ and FREQUP, the 
rescheduling of a routine is skipped if the simulation time has reached the 
Preprocessor data set stoptime. 

HDLNS 

Subroutine HDLNS (HDLTIM , NDATA, MESAGE) 

HDLNS may be called during a simulation to print messages and data and to 
store the messages and the data for a summary reprint at the end of the 
simulation, HDLTIM is a time to be printed with MESAGE. MESAGE is assumed 
to be a 8(6H) print line. NDATA specifies the number of COMMON block SCRTCH 
locations (0 to 10) to be printed on the line following the message. HDLNS 
also writes the message and the data onto unit 11. However, if NDATA- is nega- 
tive, no data is expected, sod MESAGE is written onto unit 11 without being 
printed. 

HEDLNP 

HEDLNP is an entry point in subroutine HDLNS. HEDLNP is called by TERMIN 
for the reprint of all data and messages at the end of a simulation. 
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INOUT 

Subroutine INOUT (ARRAY, ISTART, I STOP) 

INOUT Is an entry point In subroutine INOUT and prints -the variables ARRAY 
(ISTART) through ARRAY (ISTOP), These variables are printed using a floating- 
point format, five words to a line, omitting lines where all variables are 
zero. 

INPUT 

INPUT Is a subroutine called by EXECT and RECYCL to read free-fleld simulator 
Input cards. INPUT loads any of the COMMON blocks whose names are recognized 
by SOAP. INPUT may be called by the flight control , but all possible simulator 
initialization should be accomplished with the EXECT or RECYCL calls. 

INPUT2 

INPUT2 is an entry point in subroutine INPUT and is called by EXECT and RECYCL 
to read free-fi eld simulator input cards for restart and recycle simulations, 
respectively. ' . 

ISPPRT 

Subroutine ISPPRT is called by several Executive routines for diagnostic 
prints. It prints the current top wait list entry, the current flight time, 
and the elapsed computer time since the simulation started. This routine 
is primarily used for checkout runs to verify proper flight control sequencing. 

JSPPRT 

JSPPRT is an entry point in subroutine ISPPRT. JSPPRT is similar to ISPPRT 
except that the Executive routines call it prior to an Environment update. 
Instead of the wait list entry being printed, the mnemonic CALENV is printed 
to identify the Environment update. 

KALNOW 

KALNOW is a scheduled subroutine which sets the immediate update flag, LECC0M(704) 
to 1, forcing a call to the Environment. It is scheduled by CMDFIL. 
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MAIN 


MAIN Is the main program of the SOAP simulation. The simulation is Initialized 
by MAIN calling EXECT. The simulation is started by’ calling EXECTR. Once 
EXECTR Is called, control never returns to MAIN. 

t, 

NEXTIM 

Subroutine NEXTIM (T, DT, PDT,.. PC, TREF, TNEXT) 

NEXTIM is used by the system scheduling logic to calculate the next time (TNEXT) 
a task is to be executed. The technique used is to multiply the delta time step 
(DT) by a pass counter (PC), The pass counter for a task is incremented by 1.0 
each time NEXTIM is called, If the time step has changed then the previous time 
step (PDT) and the time of that change (TREF) are saved, A! the time of the 
change, the pass counter (PC) is reset to 0.0. TNEXT is then computed as above 
and the TREF is Included. 

NINOUT 

Subroutine NINOUT (ARRAY, ISTART, I STOP) 

NINOUT prints the variables ARRAY (ISTART) through ARRAY (ISTOP) much like 
INOUT. However, NINOUT uses a variable format for the print instead of the 
standard floating point of INOUT. For each variable, NINOUT determines the 
type of format, to use. A floating-point format is used if the data is 
floating point, and integer if it is integer. Otherwise, a variable octal 
format is used with the Hollerith representation printed with the octal 
number. The octal format is the minimum format required to print the. number. 

As in INOUT, the variables are printed five to a line with lines consisting 
of all zeroes being omitted. 

PINOUT 

Subroutine PINOUT (ARRAY, ISTART, ISTOP, IBLOCK) 

PINOUT is another print output routine, PINOUT first prints the array 
identifier IBLOCK, a* six-character or less BCD name, INOUT is then called 
to print ARRAY (ISTART) through ARRAY (ISTOP). 
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PITOUT * 

Subroutine PLTOUT (N, IBCD, VAL, IFOR. LU, IPASS) 

PLTOUT is used by the preprocessor-generated plot routines to dump data for 
subsequent plotting and may also be called by the flight software for the 
same purpose. N Is the number of variables to be dumped and LU is the logical 
unit onto which the data is to be dumped (default is 20). IBCD. VAL. and IFOR 
are three arrays of length N. VAL contains the data to be plotted} IBCD the 
BCD (six characters or less) identifiers for the data; and IFOR the codes for the 
format of the data which will be dumped, The IFOR codes are one for floating 
point, two for Integer, and three for octal. Prior to dumping the data, 

PLTOUT converts integer and octal data to floating point* If a flight soft- 
ware routine uses PLTOUT, it should set IPASS to a one for the first pass. 

PLTOUT will then dump the IBCD array onto LU for later use by the Plot 
Processor and reset IPASS to a zero. An additional feature for PLOTR is 
that an Environment update will be made prior to the plot dump if a simulator 
switch Is set (LECCOM (579)=1). For the PLOTR! through PL0TR9 routines, LECCOM 
(621) through LECCOM (629) contain the output logical units (default is 20) 
and LECCOM (611) through LECCOM (619) contain the Environment update switches. 

PNNOUT 

Subroutine PNNOUT (ARRAY, I START, ISTOP, IBLOCK) 

PNNOUT is an entry point in subroutine PINOUT and performs the same function 
as PINOUT, except subroutine NINOUT is used to print the variables. 

PRTSDT- 

PRTSDT is an entry point in subroutine SYSPRT and is called to print the 
schedule data tables, which consist of routine name, start time, stop time, 
delta time, priority, time of last delta time change, and schedule pass 
counter. PRTSDT can be scheduled. 


PRWAIT 

I 

PRWAIT is an entry point in subroutine SYSPRT and Is called to print the 
waitlist, which consists of scheduled routine's position on list (ELEMENT), 
time of next scheduled execution (CTABLE), name of routine scheduled (JOBTAB), 
priority (PRIOR), and location of job in scheduled data table. PRWAIT can 
be scheduled. 

QZKILL 

QZKILL is an entry point in subroutine RECOV. QZKILL is called by the UNIVAC 
1110 system routine NERR$ whenever a fatal system error occurs (except for 
print limit or time limit). QZKILL calls RESTRT for either run termination 

or a RECYCL for .another trajectory. 

» 

RECOV * 

Subroutine RECOV is written in assembly language and is called by EXECT at the 
start of a simulation. RECOV calls UNIVAC 1110 systems routines to set up* 
the call to QZKILL in case a fatal system error occurs. RECOV is also called 
by RESTRT if the user has specified another trajectory to be run. 

RECYCL 

Subroutine RECYCL may be scheduled or called by the flight control routines to read 
additional simulator Inputs at some time during a simulation and/or to 
recycle the simulation back to a previous State Image, Initially, RECYCL 
reads a data card which defines the type of action it is to perform. In one 
case, RECYCL calls INPUT to read additional data cards and then resumes the 
simulation. The other two cases involve a simulator restart at some previously 
dumped State Image. The difference between the latter two cases is that in 
one, new State Images are written over the old ones after the restart, while 
in the second case, a new simulator State Image file is begun. The former 
case would be used for recycling to an earlier State Image to correct an 
error, while the latter case represents a series of test cases run from one 
nominal starting point. In either case, RECYCL calls INPUT for additional 
simulator input after the desired State Image is obtained. For all three 
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types of RECYCLE action, optional State Image dumps are made prior to reading 
the first Input card and after INPUT has been executed. Also, RECYCL calls 
FREQEJ to be rescheduled If appropriate, 

* 

RESTRT 

Subroutine RESTRT is called by QZKILL whenever a fatal UNIVAC 1110 EXEC 8 
error occurs. Normally RESTRT calls TERMIN to terminate the simulation properly. 
However, If LECCOM (580) is set to a 1, RESTRT calls RECOV to reestablish the 
termination linkages and then calls RECYCL for another trajectory. 

REWTLT 

Subroutine REWTLT (JOB, T, ISWITC, NENTRY) 

REWTLT provides the flight control with the capability of altering and/or 
obtaining information from the waitlist, depending upon the setting of ISWITC, 
JOB is a sir-character or less BCD name of the flight control routine to be 
considered, if ISWITC = 1, the presently scheduled waitlist time for JOB is 
replaced with T. If JOB does not appear on the waitlist, it is scheduled to 
occur at T, If ISWITC = 2, JOB is removed from the waitlist. If ISWITC « 3, 
the presently scheduled time of JOB is returned to T, and the waitlist is 
not altered. If JOB does not appear on the waitlist, T is set to -1.0. On 
every call to REWTLT, NENTRY is returned to the caller as the present 
position of JOB on the waitlist or set to a -1 if JOB does not appear on 
the waitlist. All three options are applied to the first occurrence of JOB 
in the waitlist. 

SANDF 

Subroutine SANDF (IFLAG, PRIOR, UPDT, TSTOP) 

SANDF provides the Executive and flight; control with a means of obtaining 
scheduling information about the flight control routine currently being 
executed. PRIOR is the routine's priority, UPDT the update frequency, and 
TSTOP the routine's stoptime. If IFLAG equals 0, it is set to a 1 by SANDF 
to indicate completion of the call and thus may be used as a first pass flag 



/ by the caller. If I FLAG Is not equal to 0, the location of the current job's * 

V bulldit data In the bulldft data table Is returned. * 

. . m 

. SCHED 

Subroutine SCHED (TSCHED , JOB, PRIOR, ISWITC) 

SCHED Is called by several Executive routines and may be called by flight 
control routines to perform one of three functions on the waitlist. The 
function to be performed is defined by ISWITC. If ISWITC * 0, the top flight 
control routine on the waitlist is removed. The preceding SCHED call is 
made by CONTRL during the simulation control loop If the current JOB does not 
reschedule Itself. If ISWITC = 1 , JOB (six-character BCD routine name) Is 
scheduled (Included on the waitlist) with priority PRIOR to be executed TSCHED 
seconds (flight time) in the future. If ISWITC a 2, JOB Is scheduled with 
priority PRIOR to be executed at TSCHED. The last two SCHED functions are 
intended for use by the flight control 

SELECT 

SELECT is a FORTRAN subroutine written by the Preprocessor for use in the 
subsequent Simulation. SELECT is called by CONTRL for the scheduled flight 
software execution in the simulator control loop. Entry points in SELECT 
represent the print routines and/or one of the plot routines PLOTR, PL0T1 , .... 
* PL0TR8* PL0TR9 requested by Preprocessor input. These entry points may be 

called by the flight control although the entry point SELECT must never 
\ be called except by CONTRL. Also, Entry point MAXMIN may be represented in 

r SELECT. The entry point MMPRNT for printing the maximum/minimum values is an 

entry point in select. 



SI DUMP 


Subroutine SIDUMP may be called by the flight control or scheduled to dump 
State Images onto tape or mass storage for later editing or restarts, SIDUMP 
is called by EXECT; RECYCL, and TERM IN. If scheduled, SIDUMP calls FREQEJ to 




r 


i 


i 
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SIPRIN 


Subroutine SIPRIN may be called by the flight control or scheduled to print 
State Images for diagnostics. SIPRIN is called by EXECT, RECYCL , and TERMIN 
and uses PNNOUT to accomplish the print. If scheduled, SIPRIN calls FREQEJ 
to be rescheduled and for termination. SIPRIN contains 48 entry points which 
may be scheduled or called for individual COMMON block prints. 

S0RT3V 

* 

Subroutine S0RT3V (ARRAY1 , ARRAY2, ARRAY3, NVALUS) 

S0RT3V is used by SCHED to order the waitlist but may be called by a flight 
control routine, S0RT3V sequences the first NVALUS locations of the three 
arrays ARRAY1, ARRAY2, and ARRAY3 in ascending order of the values In ARRAY1. 

SREAD 

Subroutine SREAD (LU) 

SREAD is used by the Executive to read in the State Image dumps from a logical 
unit (LU), which is normally a 7 for restart runs and an 8 for recycle runs. 

SUPTIM 

Subroutine SUPTIM (ICAU, I CCER, 10) 

SUPTIME computes the CAU, CC/ER, and I/O standard units of processing (SUP) 
charges accrued for the total runstream time to the time of the call. The 
routine is UNIVAC-dependent and is the basic SOAP timing algorithm. The basic 
time increment is 200 microseconds. 

SWRITE 

Subroutine SWRITE (LU) 

SWRITE is used by the executive to write out the State Image dumps to a an LU, 
normally logical unit 8. . -- ' 
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SYSPRT 


SYSPRT is a SOAP executive system print routine* which can be scheduled or 
called. SYSPRT will print the schedule data tables and the current waitlist. 
These two prints are also accessible via entry points PRTSDT and PRWAIT, 
respectively. See these entry points for description of contents of the 
schedule data table and the waitlist prints* 

TERMIN 

TERMIN is the simulator termination subroutine and may be scheduled or called 
for termination of the simulation, TERMIN should be called for all abnormal 
simulator terminations as it terminates plot and State Image tapes which may 
then be used for diagnostic purposes. TERMIN first calls MMPRNT to print the 
maximum/minimum summary if MAXMIN was scheduled. TERMTN calls HEDLNP, an 
entry point in HDLNS, for a reprint of any "headlines" printed during the simu- 
lation. TERMIN calls SIDUMP for a final optional State Image dump, SYSPRT for 
a final scheduling summary and SIPRIN for a final diagnostic State Image print. 
TERMIN then puts end-of-flles on the plot and State Image units and terminates 
the simulation. . 

TIMER 

Subroutine TIMER (CALLER, ISWTCH, I END) 

TIMER is called to obtain a measure of the real time consumed by a particular 
routine during a simulation run, As many as 30 routines may be timed on the 
same run. The calling routine must supply the three calling arguments. CALLER 
is the name of the calling routine (Hollerith). ISWTCH is set to integer 1 to 
turn TIMER on, and then back to 0 to turn TIMER off. 

In the most simple case, the calling routine would call TIMER at the beginning 
of its execution with ISWTCH = 1, and then at the end of its execution with 
ISWTCH = 0, However, if during its execution the routine to be timed passes 
control to some other routine, it may be desirable to turn TIMER off during 
this time and then back on when the routine being timed resumes execution. 
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3.3 SOAP SYSTEM COMMON BLOCKS 


The SOAP executive recognizes the following COMMON block names: 


ACCELC 

BENDC 

DPFSC1 

GRADC 

LNDGRG 

RCS2C 

ACSC 

CONST 

ENVCOM 

GRAVC 

MASPRC 

RENDEC 

ACTVEC 

DATBUS 

FSCOM 

GUSTC 

MASP1C 

RMC 

ACTV1C 

DAVEHC 

FSC1 

GYR01C 

MASP2C 

RNGNAC 

ACTV2C 

DAVE1C 

FSC2 

GYR02C 

NAVA I C 

SCRTCH 

AEROC 

DCONST 

FSC3 

HORIZC 

ORBALC 

SLOSHC 

AER1C 

DGRAVC 

FSC4 

H0RZ1C 

ORBPAC 

SNSCTS 

AER2C 

DIMUC 

FSC5 

H0RZ2C 

PASVEC 

STARTC 

AIRDAC 

DMASSC 

FS1C 

IMUC 

RADALC 

THRU3C 

ATMOS C 

DPDBC 

FS2C 

LECCOM 

RCSC 

TWCC 

BARALC 

DPFSC 

GGTORC 

LNDAIC 

RCS1C 

WINDC 


Most of these blocks are associated with Environment models and the definitions 
of the variables stored therein are to be found in the documentation for the 
Environment models. These definitions can change from one version of the model 
to another. The COMMON blocks whose names start with the letters FS are 
associated with flight control; the definitions of the locations in these 
blocks are documented in section 5. 

The SOAP system maintains the blocks CONST, DCONST, ENVCOM, LECCOM, and 
SCRTCH. The block SCRTCH is used for temporary storage during execution. No 
fixed locations are defined. Maps of the contents of the other four COMMON 
system blocks follow. 
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3.3.2 OCONST COMMON BLOCK LENGTH = *»8 PAGE 

LOC VAR NAM DIM UNITS DEFINITION AND INITIAL VALUE 


M 


t 


U 

U. 

0 

1 


o 



a 


4 


a 

a 


OJ 

a o 


2 


a 


H 

H 

X 

o 

M H 


O 


00 


o 

O 

H 

1 

O O 


o 


K> 


40 

4/) 

a 

Cj O 

V) 4/>3 


H « 

OOO 

K 


CL 

CL 

< 

aou 

CL an 


1* H 

1 Ut' 

Kl 


H 

H 

LJ 

hhq: 

M H 4 


4 oat 

• If 

O'* 


-J 

-J 


Oir>< 

H_Jfr*JQ 


040 

ooh 

O' 


-J 


Ul 

00V1MD 

4 JHJH 

*- 

HHOU 


O' 


U 

U 

X 

OCL HO 

OW ♦ Wv> 


c$cr(sj 

OOO 

3* 

H 



H 

4 MOO) 

Ul D H 


4 vo oiy; 

aoo 

Ot- 


CL 

a 


O *JU3 NUIXOUIH 

h 

0(0 0)0 

• • » 

a'S: 

< 

U) 

UlOU 

in-iroxoKiUioJUj^r 

1 1 Ul 

fTOCMiJ 

QHO 

Ch«* 

H 

X 

ICO 

NUlCM- 1 

0) X UO X H 

iT2 

MO^OOX 



V) 

o 

U 4 


o ojao 

•uituh 

xao 

fOHfTU 

DOOCH^) 

55 

4.0 

ViCMA 

HOTCO-lH 

OMineo 

H 1 U 

ano» 

CJOOH Z 

O 

H 

M K*D 

o uiauJcr Hmcmmco 

HO 

HKhO 

• • •4 110 

u 

U 

U* OH 

oxro to 

u. mu.fr 


M 


OCJ«*-«OOKH\I 

.jcoa: 4 uunmo 
wooaroeoojLj 

cOHOOJfnHD H 

or vd-xo^cmooo a 

*1 03H O CMSit^Oe h 

.jojMOinoinuj <1 

t5*~l<r-l,3- • X C 

2Tino:0^!niD I >* 

H 03 «H(\f QH 

Xfox<4 a jhh 

HO'h •!) Si^i «Z 

arsiaro # «x- u 

m hi: a 

u •uti a(\ia</. mm 


HpCQ UD CM aOO 

au rj^HU. ua< 
• OJ 4 «a fO ohq; 
HUKD20 


• OK>UI3*4n *-* 03 

(Sli/KSIXa'OU.auU 
HinninHoaoh 
iiu.ro o)o h uo 
kh u o *x wo on cna 


lon-j 

M U.OChMCM«W 03 HO) *4 JHu. • 0(SJ UM »A) M 3* 

oho^ho' xaxMHOior co xoxu 

a H CO <O«0t 4 <00<V 4 O X*£>_J«XiT *lcr 
IUH <<M~<7' Cf MOOHX M h-H«t • • 

UiOH ffMO'OJOJQI^HOJH HNLil CL 
2LOH fJ>h-CO)Oc:MO) H^OJCICM OH 
OCL> H<H-)H2HDHU.UMHrOO*3 2 
2 •Ih*CLCSi<tOM^O«30M ttfCZoaHa 

uickisj on XHXinuir* hcoH'O xujxuj 
uocofo o> iM.i'w non.C'Z^i • a: i a 
C< «2 »H • M *Z •UMCOUJ • HH<M< 
ttZ2HOtf lg3XO)<UlUjSOQaZDL o 
«toD o u u ui q, j»hu ouauo 

UlOVOil nuoll OHXll lOLu^UJII UJiOtr.oOLO 



csi 



cv» 

(Si 









* 



# 

# 








o 

# 



* 

4 








U* 

i/* 




4A 








uo 




V 









X 

ro 



to 

tY 






CJ 

CO 

o 

* 



# 

* 




in 



* 

< 




* 





V 




a 

X 

l • 1 • 

l 

1 1 x 

X 

X 

X 

X 

X 1 

• 

XX 

X 

r-4 

H 


o 

HH*H 

H 

H 

H 

H 

PlH 

H 

viM 

H 



COO 

o 

>UO 

X 




UIH 

(SI 




LJ 

nqu 

M 

too 

o 




H <1 

a 

ccrsi 

ro 

U 

O 

HMHO 

v 

MUlX 

X 

U 

U 

u 

Hj 

a 

SL 

u 

u 

Ul 

a CL a. a 

i 

ZQCU. 

u 

<: 

GO 

X 

JL 

Ul 

U< 

CD 

CL 

CL 

a cl a cl 

CL 

UwOOl 

a 

CL 

CL 

CL 

ua 

a. 

ua 

a 

o 

a 

oooo 

a 

OOO 

o 

O 

a 

o 

QD 

a 

oo 

a 

H 

csr 


4- 

»oh® 

a 

o 


ru 

K|j» 

in 

0)H 

€C 





. i 4» <« < 

H: 

ra 

(Si 

(Si 

coca 

N 

COM 

(SI 


3-27 


DCONST COMMON BLOCK CONTINUED PAGE 

LOC VAR NAM DIM UNITS DEFINITION AND INITIAL VALUE 
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LECCGM COMMON BLOCK-VERSION 800 CONTINUED 
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3.4 MATH. FLIGHT CONTROL. AND ENVIRONMENT UTILITY ROUTINES 

These routines supplement the FORTRAN mathematical library and serve to support 
Environment and/or Flight Control routines. Typically, the user will call these 
routines rather than schedule them. In places where both single- and double- 
precision versions of the routines are available, one of the writeups Is consid- 
erably more extensive, None of these routines has a dedicated COMMON block, 
although some use information stored in specified locations of other COMMON 
blocks, ' 

Routines described in this section include.* 


ABVAL 

dot 

MATYPR 

ACCOSH 

DOTDP 

MSU 

ARTAN 

DPDFQ 

MSUS 

ATOQX 

DPDFQS 

MTQ 

CROSS 

DPDTRM 

MXM 

CROSSDP 

DPEVAL 

MXV 

DABVAL 

DPOLY 

ORTHOG 

DCROSS 

DPOSSO 

PEVAL 

DDOT 

DPRK2 

POLY 

DETRM 

DPRK2S 

POSSUM 

DEULER 

DUN IT 

QCONJG 

DEtlLRD 

DVXM 

QNORM 

DEULRE 

DXPOSE 

QTM 

DIFEQ 

DXTRCT 

QTOAX 

DIFEQP 

DXYZMT 

QUTMLT 

DIFEQS 

DYZXMT 

QXFORM 

DIFQPS 

EULER 

RAP 

DMXM 

EULERD 

SGN 

DMXV 

EULERE 

SIG 

DNVRSE 

DORTHO 

EXTRCT 

INVERSE 

SIGNM 
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ABVAL 

Function ABVAL (VECTOR) 

ABVAL takes the absolute value of a single-precision three-element vector and 
returns a single-precision scalar, 

ACCOSH (X) 

Function ACCOSH (X) 

Function ACCOSH (X) returns the single-precision hyperbolic arc cosine of x 
in radians. 

ARSINN 

Function ARSINN (THETA) 

Single-precision real function ARSINN computes the angle whose sine is equal 
to THETA. ARSINN calls the system ARCSIN (ASIN) function but does some testing 
to ensure that ASIN will not abort the run. Angle output Is in radians In the 
interval from -j to +£. 

ARTAN 

Function ARTAN (TOP, BOTTOM) 

Single-precision real function ARTAN computes the angle whose sine and cosine 
are equal to TOP and BOTTOM, respectively, The system routine ATAN2 is called 
If the testing of the Input arguments are such that they will not abort the run. 
Angle output is in radians on the Interval from -tt to +tt. 

ATOQX 

Subroutine ATOQX (A, Q) 

ATOQX extracts a quaternion Q from an Euler sequence of rotations A2, A3, A1 
about the Y, Z, X axes, respectively, 
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CROSS 

Subroutine CROSS (VECTR1, VECTR2, VECTR3) 

Single-precision subroutine CROSS computes the vector cross product of vector 
VECTR1 with vector VECTR2, returning the resultant vector in VECTR3. 

CROSSDP 

Subroutine CROSSDP (VECTR1, VECTR2, VECTR3) 

CROSSDP computes VECTR3, the cross product of VECTR1 times VECTR2 (in that 
order) where all vectors are double precision. 

DABVAL 

Function DABVAL (DVECT) 

Double-precision function DABVAL returns the magnitude of the vector DVECT. 
DCROSS 

Subroutine DCROSS (DVECT 1, DVECT 2, DVEC73) 

Double-precision subroutine DCROSS computes the vector cross product of the 
vectors DVECT1 with DVECT2, returning the resultant vector in DVECT3. All 
quantities are double precision. 

DDOT 

Function DDOT (X, Y) 

Double-precision function DDOT returns the vector dot product of the double 
precision vectors X and Y. 

DETRM 

Subroutine DETRM (DETRMTy VALUE) 

DETRM evaluates the determinant DETRMT, returning the result in VALUE. DETRMT 
is a 3x3 matrix stored sequentially by columns. 
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DEULER 

subroutine DEULER (RPYH, RPYMAT) 

Double-precision version of EULER (see EULER). Computes a transformation matrix 
when given roll, pitch, yaw In that order. 

DEULRD 

Subroutine DEULRD (EULER, W, DW, DEUL, DDEUL) 

Double-precision version of EULERD. Calculates Euler rates and accelerations. 
DEULRE 

Subroutine EULRE (DT, DEUL, DDEUL, DE) 

Double-precision version of EULERE. Calculates double precision estimated delta 
Euler angle set from Euler rates and accelerations. 

DIFEQ 

Subroutine DIFEQ (NEQNS, X, DX, Y, I PHASE) 

DIFEQ solves a set of first-order differential equations using a four-pass, 
Runge-Kutta numerical integration technique. NEQNS specifies the number of 
first-order equations to be integrated. X is the independent variable, and 
DX is the integration step size. Y is an array, dimensioned 4NEQNS + 1, which 
contains the dependent variables in the first NEQNS locations, the corresponding 
derivatives in the second NEQNS locations, and a working area for DIFEQ in the 
remaining locations. I PHASE is a flag indicating which of the four passes has 
been executed and is set by DIFEQ. I PHASE is reset to 0 after the pass which 
completes the integration step. Before the first pass of an integration, 
the independent variable, dependent variables, derivatives, and the flag must 
be set. After each space, the derivatives must be reevaluated for the next 
step. 
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DIFEQP 

Subroutine DIFEQP (NEQNS, X, DX, Y, IPHASE) 

DIFEQP Is the same as DIFEQ in that it solves a set of first-order differential 
equations using a four-pass, Runge-Kutta numerical integration technique. All 
calling arguments have the same definition except for Y, which is an array 
dimensioned 5NEQNS + 1. DIFEQP differs in technique from DIFEQ in that 
during the last (fourth) pass of an update, DIFEQP calculates the roundoff due 
to computer word *ngth and saves this value, which is added back in during the 
next update, Use and Implementation of DIFEQP is the same as DIFEQ, but a 
more accurate result is obtained. , 

DIFEQS 

Subroutine DIFEQS (NEQNS, X, DX, Y, IPHASE) , 

DIFEQS is an entry point in DIFEQ and is like DIFEQ, except it is used for 

second-order equations. NEQNS is equal to twice the number of second-order 

equations for DIFEQS, Y is an array, dimensioned 4NEQNS + 1, containing the 

dependent variables in tho first NEQNS/ 2 locations, the first-order 

derivatives in the second NEQNS/2 locations and in the third NEQNS/2 

locations, the second-order derivatives in the fourth NEQNS/2 locations, 

*■ 

and a DIFEQS working area In the remainder of the array. The integration 
process for DIFEQS is the same as for DIFEQ, where the first NEQNS locations 
for DIFEQS are considered to be dependent variables and the second NEQNS 
locations are the derivatives. 

DIFQPS 

Subroutine DIFQPS (NEQNS, X, DX, Y, IPHASE) 

DIFQPS is an entry point In DIFEQP and is like DIFEQP, except it is used for 
integrating second-order equations. Definition of the calling arguments is 
the same as DIFEQP except for NEQNS, which is equal to twice the number of 
second-order equations. The Y array assignments are the same as outlined in 
DIFEQS. 
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DMXH 

Subroutine DMXM (DM ATI , DMAT2, DMAT3) 

Double-precision subroutine DMXM computes the product of the matrices DMAT1 
with DMAT2 returning the resultant matrix In DMAT3. 

DMXV 

Subroutine DMXV (DMAT1 , DVECTl, DVECT2) 

Double-precision subroutine DMXV computes the product of the matrix DMAT1 
with the vector DVECTl, returning the resultant vector in DVECT2. 

DNVRSE 

Subroutine DNVRSE (DMAT1 , DMAT2 ) 

Double-precision subroutine DNVRSE computes the Inverse of the matrix DMAT1 , 
returning the resultant matrix In DMAT2. 

DORTHO 

Subroutine DORTHO (AM) 

Double-precision version of ORTHOG. Orthogonal Izes a matrix, l.e. , creates a 
set of three orthogonal unit vectors from It. 

DOT 

Function DOT (VECTRl, VECTR2) 

DOT returns the vector dot product of VECTRl and VECTR2. 

DOTDP 

Function DOTDP (VECTRl , VECTR2 ) 

DOTDP returns the double-precision dot product of two double-precision vectors. 
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DPDFQ 

subroutine DPDFQ (NEQNS, X, DX, Y, IPHASE) 

DPDFQ is the double^precislon version of DIFEQ. All calling arguments are 
single precision except V* which is a 4NEQNS + 1 dimensioned double-precision 
array. Use and implementation of DPDFQ Is the same as DIFEQ. 

OPDFQS 

Subroutine DPDFQS (NEQNS, X, DX, Y, IPHASE) 

DPDFQS Is an entry point In DPDFQ and is like DPDFQ, except It is used for 
second-order equations. Definition of the calling arguments is the same as 
DPDFQ except for NEQNS, which is equal to twice the number of second-order 
equations. The double-precision Y array assignments are the same as outlined 
in DIFEQS. 

DPDTRM 

Subroutine DPDTRM (DETRMT, VALUE) 

Double-precision subroutine DPDTRM evaluates the double-precision determinant 
DETRMT, returning the result in double-precision VALUE. DETRMT is a three by 
three determinant stored sequentially by columns. 

DPEVAL 

Subroutine DPEVAL (A, N, X, PY) 

Double-precision version of PEVAL. DPEVAL evaluates a polynomial of N terms 
with coefficients A at X and stores the value In PX. 

DPOLY 

Subroutine DPOLY (X, FX, N, C) 

Double-precision version of POLY. DPOLY finds a polynomial of order N-l 
to relate N values of Independent variable X to N values of dependent variable 
FX. The coefficients are stored in array C. The number of points, N, must 
be 25 or smaller. 


DPOSSU 

Subroutine DPOSSU (TIte, RSMKS, RMHKS) 

Double-precision version of POSSUM. Computes position and apparent diameter 
of sun and moon using data from reference 5. 

DPRK2 

Subroutine DPRK2 (NEQNS, X, DX, Y, I PHASE) 

0PRK2 solves a set of first-order differential equations In double precision, 
using the second-order* two-pass Runge-Kutta numerical Integration technique. 

NEQNS and I PHASE are Integer variables, X, OX, and Y are double-precision 
variables. 

DPRK.2 differs from DIFEQ in that it is a second-order solution Instead of fourth- 
order, and the computations are in double precision instead of single precision. 
The definition of the calling arguments is the same as for DIFEQ, except for pre- 
cision, DPRK2 specifically uses the "Heun form" second-order Runge-Kutta solution 
During the first pass through DPRK2, time is advanced to X + |ox, and IPHASE is 
set to 1. During the second pass, time is advanced over the final one-third 
of the integration step, and IPHASE is set to 0. Prior to each pass, the calling 
routine is responsible for reevaluating the derivatives. 

DPRK2S 

Subroutine DPRK2S (NEQNS, X, DX, Y, IPHASE) 

DPRK2S is like DPRK2, except it is used for second-order equations. NEQNS and 
IPHASE are integer variables, X, DX, and Y are double-precision variables. 

Except for precision, the definition of the calling arguments is the same as 
for DIFEQS. DPRK2 and DPRK2S have the same relationship as DIFEQ and DIFEQS . 

DUNIT 

Subroutine DUNIT (DVECT1, DVECT2) 

Double-precision subroutine DUNIT unitizes DVECT1 , returning the resultant 
vector in DVECT2. • 


f 


* 


DVXM 

Subroutine DVXM (DVECT1, DMATl » DVECT2) 

Double-precision subroutine DVXM computes the product of DVECT1 with matrix 
DMATl, returning the resultant vector In DVECT2, 

DXPOSE 

Subroutine DXPOSE (DMATl, DMAT2) 

Double-precision subroutine DXPOSE transposes matrix DMATl, returning the 
resultant matrix In DMAT2, 

* 

DXTRCT 

Double-precision version of EXTRCT . Extracts an RPY Euler sequence from an 
attitude transformation matrix. 

DXYZMT 

Subroutine DXYZMT (EULERS, XMATRX) 

* 

Given the three double-precision Euler angles EULERS (In radians), double- 
precision subroutine DXYZMT defines the double-precision transformation matrix 
XMATRX, assuming an X (roll), Y (pitch), Z (yaw) order of rotation. EULERS is 
an array of roll, pitch, yaw angle sequence. 

DYZXMT 

Subroutine DYZXMT (EULERS, XMATRX) 

Given the three double-precision Euler angles EULERS (in radians), double- 
precision subroutine DYZXMT defines the double-precision transformation matrix 
XMATRX, assuming a Y (pitch), Z (yaw), X (roll) brder of rotation. EULERS is 
an array of pitch, yaw, and roll angle sequence. 
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EULER 




ROUTINE DESCRIPTION 

The EULER routine calculates a roll, pitch, yaw rotational 

* 

sequence transformation matrix. 

I 1 

MATH MODEL , ’ # 

The transformation matrix erected by EULER assumes a roll, 

f 

pitch, yaw rotation sequence. . The three angles are received 
via calling arguments and then substituted into the following 
equation. ’ i 


, RPYMAT (1) 

3 

COS (?) 

COS (Y) 



RPYMAT (2) 

- 

-COS(P) 

SIN (Y) 



RPYMAT (3) 

- 

SIN (P) 

. 



RPYMAT (4) 

XC 

SIN (R) 

SIN (P) 

COS (Y) + COS (R) 

SIN (Y) 

RPYMAT ( 5 ) 

m 

-SIN (R) 

SIN (P) 

SIN (Y) + COS (R) 

COS (Y ) 

RPYMAT (6) 

3 

-SIN (R) 

COS (P) 



RPYMAT (7) 

m 

-COS (R) 

SIN (P) 

COS (Y) 4- SIN (R). 

SIN ( Y) 

RPYMAT (3) 

at 

COS (R> 

SIN (P) 

SIN (Y) 4 SIN (R) 

COS (Y ) 

RPYMAT (9) 

3 

COS (R) 

COS ( P ) 



* *» 

where R, P, Y are the roll, 

pitch, 

yav; angles ., 



INPUT/OUTPUT 

All input/output of EULEP. is done via calling arguments. 

The call to EULER IS 

CALL EULER (RPYA, RPYMAT) • 

is the input array consisting of the roll, pitch, and 
yaw angles in that order. 

is the output array consisting of the nine elements of 
the transformation matrix. 


where 

RPYA 

RPYMAT 
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ROUTINE VERIFICATION 


This routine was verified by logic review and by comparing routine 
output to hand calculated values. 

# 

EULERD 


ROUTINE DESCRIPTION 

The EULERD routine calculates the Euler angle rates and 
accelerations* 

MATH MODEL * 

I 1 

Given the vehicle body rates and accelerations and the attitude, 
the Euler angle rates and accelerations can be calculated. 


where 


and 


4> 



p costii - q sln 'jj 

cose 

p sint}/ + q cost}/ 
r - <$ sine 


a a A • e • 

u « pi + qj+rk ■ ^ + F + t{7 


V = (p - ctl>) cos;!/ - (a » o-* 1 ) sin'!* 4- 1 9 sin8 

■ s COS& 

0 ■ (ptj/ + q) ccs^ + (p - q$) sint}/ 

w # '% 9 ■ " ' •• 

$ * r - $e cos9 - $ sin8 


The above equations become indeterminate when ‘the pitch angle, 
0, is + 90°. To Obtain* the Euler rates and accelerations at 
this point, a different set of equations is used. 
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9 

e 

9 


■ — =■ ^-p cosy + q sin y) + j 

■ p sin^ + q cos 9 

■ -4* (£> cos t|» - q sint|») + 

20 ' ' 2 


and 

to 

♦ ■ 
to 

6 ■ 

•« 

t - 


- 

(p9 + q) 
r - * 


2p) 


20 

sin^ 

0 

+ (P 


- + 2q) 

- q$) sinp 


cost!) 




Actually, the switch to the above equations for determining 
the Euler rates and acceleration is made as 0-*-+9O o . Even 
in the latter set of equations 0 can become zero in which 
case the equations below are used. 


and 


0 

to 

9 

m 

i 


6 « 

9 * - 

i * 


+ (P 2 + q 2 ) k 
r - (pq -qp)/(p 2 + q 2 )J 
r + (pc - qp)'/(p 2 + q 2 )] 


1 

r 

i 

T 


(p$ + q) cos# + (p - q#) sin# 

if: 


- [ #(q# - 2p|p — # (py + 2q)q /(p 2 + q 2 ) - &£| 


r - 


The rotational sequence • .assumed in deriving these equations 
is roll, pitch, yaw. . 



ss 


The definitions of the terms in the above equations are: 


p,q,£ - vehicle body angular rates 

p,q,r - vehicle body angular accelerations 

0 

♦ - Euler angles roll,' pitch, yaw 

i 

$ , 0 , iji - Euler angle rates . 

?,0,5 - Euler angle accelerations 

INPUT/OUTPUT 

All input and output of EULERD is via calling arguments. 

The input arguments are.: • 

$,0,if> * - Euler angles roll, pitch, yaw [PS: EULER] 

p,q,r -♦Vehicle body angular rates [PS: W] 

p,q-,r - Vehicle body angular accelerations [PS: DW] 

The output arguments are : 

• ■ • • 

- Euler angle rates [PS: DEUL] 

♦,5,$ - Euler angle accelerations [PS: DDEUL] 

where [PS : ] indicates the program symbol . 

« • 1 

ROUTINE VERIFICATION 

This routine was verified -by logic review and by comparing 
routine output to hand calculated values. 
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EULERE 


ROUTINE DESCRIPTION 

The EULERE routine calculates an estimated delta Euler 
angle set from Euler rates and accelerations. 

> 

MATH MODEL 

II <ru..niri n r ' .1 T I » 

The previously calculated Euler angle rates and accelerations 
are used to estimate the Euler angles change over some delta 
time . 

• l'« 2 

Ac a * c^AT + jc ^ A x 

where 

Ac^ ■ is the Euler angle change over the AT. [PS : DE] 

» is the Euler angle rates [PS: DEUL] 

* is the Euler angle accelerations [PS: DDEUL] 

AT » is the delta time period [PS: DT] 

INPUT/OUTPUT 

All input/output of EULERE is done via calling arguments. 

The call to EULERE is 

CALL EULERE ( DT, DEUL , DDEUL ,DE) 

where the arguments are defined above in [PS : J. 

ROUTINE VERIFICATION 

This routine was verified by logic review and by comparing 
routine output to hand calculated values. 
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EXTRCT 


ROUTINE VERIFICATION 

The EXTRCT routine extracts from an attitude matrix a set 

// 

of three Euler angles assuming a roll, pitch, yaw rotational 
sequence. 


MATH MODEL 

The EXTRCT routine assumes the attitude matrix was formed 


using a roll, pitch, yaw rotation sequence and to have the 
following form: 


COS (?) COS (Y) 


SIN (R) SIN (P) COS (Y) 
+COS (R)SIN(Y) 


COS (R) SIN (P) COS (Y) 
+SIN(R) SIN (Y) 


COS (P)SIN(Y) 


-SIN (R) SIN (P) SXN(Y) 
+COS (R) COS (Y) 


COS (R) SIN (P) SIN ( Y) 
+SIN (R) COS (Y) 


SIN (P) 


-SIN (R) COS (P) 


COS(R) COS(P) 


where R,P,Y are the roll, pitch, yaw Euler angles. 


To find the pitch angle, the arcsine is taken of matrix 
element 3,1. The yaw angle can be found by dividing the 
matrix elements 1,1 or. 2,1 by the cosine of the pitch angle 
and then taking the arccosine or arcsine respectively. The 
roll angle is found by the same procedure using matrix 
elements 3,2 and 3,3. ‘ . , 

As the pitch angle approaches -90'*, the matrix approaches 
a singularity. When this situation arises, the delta Euler 
angles calculated by EULERE are added to the old values to 
obtain an estimated value. 
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INPUT/OUTPUT 

I * 

All input/output of EXTRCT is done via calling arguments. 

The call to EXTRCT is 

.CALL EXTRCT (AM,DE,A) 

where inputs are: 

AM - is the attitude matrix 
DE - is the delta Euler angles 

output is: 

A - is the attitude array of .roll (R) , pitch (P) , 
and yaw(Y) angles. 

ROUTINE VERIFICATION 

This routine was verified by logic review and by comparing 
routine output to hand calculated values. 

INVRSE 

Subroutine INVRSE (XMATRX, YMATRX ) 

INVRSE computes the inverse of matrix XMATRX, returning the result in matrix 
YMATRX. 


ROUTINE DESCRIPTION 


The MATYPR routine extracts from an attitude matrix a set of three Euler 
angles 'assuming a yaw, pitch, roll rotational sequence. 


MATH MODEL 


The MATYPR routine assumes the attitude matrix was formed using a yaw, pitch 
roll rotation sequence and having the following form: 


cos. e cos p 


cos e sin p 


•cos p sin p cos $ cos p 

+ sin p sin 9 cos p + sin P sin 9 sin p 

sin p sin ^ -sin p cos p 

+ cos $ sin e cos tj; + cos $ sin 9 sin p 

where p, 9, p are the yaw, pitch, roll Euler angles. 


-sin 9 
sin $ cos 9 
cos p cos 9 


The pitch angle Is calculated as follows: 

P - ASIN t-sin 9) 

The yaw angle Is calculated as 

w _ cos 9 cos p 
*1 V cos P 

„ _ cos 0 sin d; 

or h COS P 

and finally Y » ACOS Y 1 

or Y • ASIN Y 2 

t * 

Likewise, the roll angle Is calculated as 

„ _ cos 9 sin 
R 1 * ■ <!osT 

B _ COS 9 COS ± 

or k 2 * cos P 
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where* 


R - ASIN R 1 

i 

R • ACOS R 2 

* 

Where Y, P, R are the yaw, pitch* roll Euler angles. 

As the pitch angle approaches ± 90°, the matrix approaches a singularity, 

When this situation arises, the delta Euler angles calculated by EULERE are 

# 

added to the old values to obtain an estimated value (ref. 4), 

INPUT/OUTPUT 

All input/output of MATYPR Is done via calling arguments. 

The call to MATYPR Is • 

i i 

CALL MATYPR (AM, CL, A) 

where Inputs are: 

AM Is the attitude matrix. 

0E Is the delta Euler angles. 

Output Is: 

A Is the attitude array of yaw (Y ) , pitch (P), and roll (R) angles, 

ROUTINE VERIFICATION 

This routine was verified by logic review and by comparing routine output to 
previously computed values. 

0 

SOURCE ■ . - ‘ 

' • 

MATYPR was derived primarily from reference 4. 




jtO0L ' .JLrh. . 


: \ 
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MSU 

Subroutine MSU (AN6LS1 , ANGLS2, ANGLS3) 

MSU performs the angular modular subtraction of the angle sets ANGLS1, ANGLS2, 
returning the results in ANGLS3. ANGLS1 , ANGLS2, and ANGLS3 are dimensioned 
.three and are in radians. ANGLS3 are angles in the range +tt with a counter- 
clockwise rotation positive. 

MSUS 

Subroutine MSUS (ANGLS1 , ANGLS2, ANGLS3) 

MSUS performs the angular modular subtraction of the single angles ANGLS1 and 
ANGLS2. The resulting angle is returned in ANGLS3. ANGLS1 , ANGLS2, and ANGLS3 
are single angles and are in radians. ANGLS3 is an angle in the range +ir with 
a counterclockwise rotation positive. 

MTQ 

Subroutine MTQ (M, Q) 

MTQ extracts a quaternion from a direction cosine matrix. The extraction process 
avoids singularities. A positive scalar part is arbitrarily chosen for the re- 
sulting quaternion. It is assumed that the input matrix Is unitary. 

MXM 

Subroutine MXM ( XMATRX , YMATRX, ZMATRX ) 

MXM computes the product of matrix XMATRX with YMATRX, returning the resultant 
matrix in ZMATRX. 

MXV 

\ 

Subroutine MXV (XMATRIX, VECTR1, VECTR2) 

MXV computes the product of matrix XMATRX with vector VECTR1 , returning the 
resultant vector in VECTR2. 
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ORTHOG 


ROUTINE DESCRIPTION 


The ORTHOG routine orthogonalizes an attitude matrix. 
MATH MODEL 

The attitude matrix [AM} is assumed to be defined by 




To orthogonalise [AM] it is first transposed. 


[AMT] 


[AM] T - [AM 1# AiM 4 , Afi ? ] 


Then, 


MS. * AMT 4 


The vector components AMT 7 and are unitized, and 


5M? 4 » AMT ? * A*'1T X 


The matrix [AMT] is now transposed back to [Af-3] . 
INPUT/OUTPUT 

The matrix [AM] is input, orthogonalized and then output 
through one calling argument to ORTHOG. 

ROUTINE VERIFICATION 

This routine was verified by comparing routine output to 
hand calculated values. 


1 * , 


PEVAL 

ROUTINE DESCRIPTION 


PEVAL evaluates a given polynomial at a given point. Used 
in conjunction with the curve-fitting utility routine POLY 
PEVAL provides a means of extracting current values from 
tabulated input. 

MATH MODEL 
The polynomial 

1 2 n—1 

P (x) ■ a Q + a^x + a 2 x + ... + a R _^ x 

*. 

is evaluated by Horner's method, defined recursively as: 

V*> * s n-l. 

V x > - p i-l (x) - x + 

with P(x) ■ P n _ 1 (x) 

INPUT/OUTPUT 

Subroutine PEVAL is referenced by 
CALL PEVAL (A,N,X,PX) 
with arguments defined as: 

A - array of polynomial coefficients: A(I) » a i-i 

N - dimension of A (N * degree (P) + 1) 

9 % 

X - independent variable for which polynomial is to 
be evaluated 

PX - resultant value of dependent variable 

• . * 

PEVAL uses no COMMON input, contains* no internal data, and 
makes no external references. 
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ROUTINE VERIFICATION 


Output from PEVAL was compared with hand-calculated values. 

POLY 

ROUTINE DESCRIPTION 

■ * 

POLY fits a polynominal curve to a set of tabulated data 

» 

points. Used in conjunction with the polynominal-evaluation 
utility routine PEVAL, POLY provides a means of extracting 
current values from tabulated input. 

MATH MODEL 

Given n data points (Xj # 1 (x^)) # POLY approximates the function 
f with the unique polynomial of degree n-1 or less defined by 

P(x) - 2 (x) tt) 

i-1 1 1 

where L.^ is the ith Lagrangian polynomial 

(x-x^)^ . . (x-x^) (x-x i+1 ) . . . (x-x n ) , 

lit (x) ■ : 

* • ‘(x^-Xj^) . . . (x i -x i-1 ) ( X i- X i +1 ) • • • 

' I 

This polynomial will duplicate the function f at the specified 
data points (except for round-off error) , but may vary from f 
between input points (as does linear interpolation) . Since 
intermediate values in the computation of the Lagrangian 
polynomials may become many orders of magnitude larger than 
the input values, double precision variables are used internally 
to avoid floating point overflow and underflow. However, 
double precision overflow is still a consideration. If POLY 
is called with too many data points, overflow and/or underflow 
may cause divergence from the input function f , in which case 
POLY should be called with fewer points specified. 
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POLY returns the polynomial P <x) as defined by equation 1 in 
the form 

P(x) ■ a Q + ajX + ... + (2) 

. ♦» 

INPUT/OUTPUT 

Subroutine POLY is referenced by 
CALL POLY (X,FX,N,C) 
with arguments defined as; 

X - an array of N distinct values of the independent 
variable; their order is irrelevant# but no two 
may be equal. * 

PX - an array of the N values of the dependent variable# 
corresponding by index to X - i.e., FX(X) * f (X(I) ) 

i 

N - the number of data points specified. 

C - the output array of coefficients of the returned 
polynomials C(I) ■ in' equation 2. 

, POLY uses no COMMON input and makes no external references. 

Internal arrays are dimensioned 25; this value must be changed 
if POLY is to be called with more than 25 data points# (This 
practice is not recommended due to the overflow considerations 
mentioned above). 

• * 

ROUTINE VERIFICATION 

■ ,v 

POLY was checked out both independently and in conjunction 
with utility routine PEVAL, with output coefficients and 
values compared with hand-calculated results. 
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P0SSIJM/80 


Versions: 


/80SP single precision 

/80DP double precision 


POSSUM is a routine which computes the earth-centered-inertial (ECI) Cartesian 
positions of the sun and moon, ft is meant to be called either from the gravity 
force model GRAV2 (with lunar and solar perturbation flags on) or from the Rela- 
tive Attitude Processor (RAP), a flight software routine. POSSUM is designed 
to use orbital data for the sun and moon as published in the American Ephemeris 
and Nautical Almanac (ref. 5) in years up to and including 1980. This data is 
in the form of orbital elements; semimajor axis, a, eccentricity, e, inclination, 
i , mean motion, M, longitude of ascending mode, ft, longitude of perigee, to, and 
longitude of sun or moon at epoch, t. 


Sample orbital element data is displayed in figure 3-6, The POSSUM requires that 
the epoch for the elements be input in the form of year since 1900, month (two 
digits), day (two digits). The date for the start of the simulation is needed 
in the same form. This routine converts these dates into day counts since 1900 

for both the epoch and the date of simulation. The difference between the epoch 

■ ■ ' ' . : : ' ' 1 . ■■ 

and the date of simulation is computed and a warning issued if it is inconsistent 

With the availability of new orbital data (elements are updated roughly monthly). 
The fraction of a day at the start of the simulation is then added to the number 
of days elapsed since epoch. Integer numbers of lunar and/or solar periods are 
removed from this elapsed time. 


On the first pass through POSSUM, the epoch is used as a reference time; on 
successive passes, the time at the last call is used for a reference. In either 
case, the true anomaly at epoch is computed using v Q -I - w. The eccentric 
anomaly at the corresponding time is computed via (ref. 6) 


c _ e + cos v. 

cos t Q = o 


1 + e cos v„ 


Kepler's equation is iterated using 

AE - AtM - (E t + e sin E f ) + (E + e sin El 

t X 0 0 


( 1 ) 


( 2 ) 
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where E t is a trial value of the eccentric anomaly. 

When AE falls below some tolerance, E f is adopted and a true anomaly is computed 
via 


cos v - e - cos E (3) 
e cos - 1 

At this point the distance to the sun or moon is computed assuming an elliptical 
orbit. The Cartesian geocentric equatorial position components are obtained by 
first finding ecliptic Cartesian coordinates (ECC) 


x EC c a ( cos cos * e ) R 

Y ECC s *e cos V R 

Z ECC = ^ R 

where 

is the ecliptic latitude 

is the ecliptic longitude 

R is the distance from the origin 

Figure 3-7 displays the geometry. From the law of sines 



sin A e - sin i sin (v + w) 

(5) 

while from the law of 

cosines 


.£ 

e = cos~^ £cos (v + w)/cos ,\j + f? 

(6) 

(see fig. 3-7) and 

0 _ >(1* e 2 ) 

(7) 


~ 1 + e cos \) 


The ECC coordinates are converted to ECI by rotating about the x-axis by e, 
the obliquity of the ecliptic. The rotation matrix is 



0 0 
cos e -sin e 
sine cose 
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( 8 ) 




and 


* X ECI 


’ X ECC* 

Y ECI 

= B 

Y ECC 

• Z ECI. 


. Z ECC. 


The solar coordinates obtained from equation (9) are relative to the barycenter 
of the earth-moon system. They are corrected back to the center of the earth 
by 

W 

h - v s - V“ do) 

z s = z s ' V" 

where 

jj is the earth- to-moon mass ratio 
s denotes solar 
M denotes lunar 

Unit vectors from the center of the earth to the sun and moon are then construc- 
ted. The apparent diameters of the sun and moon are also computed. 

This routine has been checked by comparing the resultant values with tabulated 
apparent semidiameters and directions (right ascension and declination). 
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INPUT 


Program 

symbol 

Common 

block 

name' 

Common 

block 

location 

Dimen. 

Description (units) 

TIME 

(Calling argument) 


Time start of simulation (sec) 

TLAST 

GRAVC 

11 


Time at which POSSUM was last 
called starting from the beginning 
of the simulation, computed 
internally after first pass (sec) 

I DATE 

GRAVC 

55 

56 

2 

Epoch for which the orbital elements 
of sun and moon are defined, then 
date of simulation. Both dates 
are in the format; year (two 
digits), month (two digits), day 
(two digits). (YYMMDD) 

FDATE 

GRAVC 

57 


Fraction of a day at start of 
simulation (days) 

A 

GRAVC 

58 

59 

2 

Semimajor axes of sun and then 
moon geocentric orbits (m) 

ECC 

GRAVC 

60 

61 

2 

Eccentricity of solar then lunar, 
geocentric orbit 

XINC 

GRAVC 

62 

63 

2 

Inclination of solar’ then lunar, 
orbit relative to the ecliptic (deg) 

XMMOT 

GRAVC 

64 

65 

2 

Mean motion of sun, then moon, in 
orbit (deg/ sec) 

PI AM 

GRAVC 

66 

67 

2 

Apparent diameter of sun, then moon, 
at a distance equal to orbital semi- 
major axis (deg) 

X0ME6 

GRAVC 

68 

69 

2 

Longitude of the ascending mode of 
orbit of sun, then moon, measured 
in the plane of the ecliptic (deg) 

WOMEG 

GRAVC 

70 

71 

2 

Longitudes of perigee for sun, then 
moon, measured from vernal equinox 
to node in the ecliptic plane, from 
node to perigee in orbital plane (deg) 

XLONGO 

GRAVC 

72 

73 

2 

Longitude at epoch of sun 5 then moon 
(longitude of the sun differs by 180° 
from longi tude of the earth) (deg) 
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Program 

symbol 

EL 

EU 

UV 

I FIRST 
IMNTH 

OBLQ 

EMMR 

OUTPUT 

Program 

s ymbol 

RSMKS 

RMMKS 

GPOS 

UV 

API) I AM 
R 


Common 

block 

name 

Common 

block 

location 

Dimen. 

• Description (units) 

GRAVC 

74 

75 

2 

Lower limit of eccentric anomaly for 
sun, then moon (rad) 

GRAVC 

76 

77 

2 

Upper limit of eccentric anomaly for 
sun, then moon (rad) 

GRAV 

78-86 

3,3 

Unit vectors from center of earth to 
sun (UV (1,1) - UV(3,1) ) and to moon 
(UV( 1,2) - UV(3,2) ) 

GRAVC 

87 


First pass flag 

GRAVC 

88-99 

12 

Number of the day of the year on 
which each month starts 

GRAVC 

100 


Obliquity of the ecliptic (deg) 

GRAVC 

101 


Earth-to-moon mass ratio 

Common 
bl ock 
name 

Common 
bl ock 
location 

Dimen. 

Description (units) 

(Calling argument) 

3 

Position of sun in EC I coordinates 
(m) 

(Calling 

argument) 

3 

Position of moon in EC I coordinates 
(m) 

GRAVC 

12 - 
17 

6 

GPOS(l) = RSMKS ( 1 ) 

GPOS (4) = RMMKS ( 1 ) 

These locations are used to store 
positions previously computed (m) 

GRAVC 

78- 

86 

3,3 

Unit vectors describing the position 
of the sun, then the moon 

GRAVC 

102- 

104 

3 

Apparent diameter of sun then moon, 
then earth (deg) 

GRAVC 

105 

2 

Distance of sun, then moon from 
earth (m) 


POSSUM/81 

% 

Version: /81SP single precision 

/81DP double precision 

POSSUM is a routine which computes the ECI Cartesian coordinates of the sun 
and the moon. It also computes unit vectors to the sun and moon and apparent 
diameters for each in degrees. This information is meant for use by the 
gravity force model GRAV2 and by the flight control utility routine RAP. 

The POSSUM routines use locations in the common block GRAV to store output. 

POSSUM stores an epoch for solar positioning (year* month, day) and a 
starting time for the simulation (year, month, day and fraction of a day 
in Universal time). On the first pass it computes the number of days elapsed 
between the epoch for solar constants and the time of the start of the 
simulation. At times subsequent to the start of the simulation, the time 
between the epoch and the Universal time at that point in the simulation 
is computed (in days). The low precision formulae for solar coordinates 
given in section D of The Astronomical Almanac (ref. 7) are used. These formu 
lae are expected to be accurate to 0.01° in right ascension and declination. 

Defining 

G = mean anomaly of sun 

d - time elapsed since epoch 

SC^j = constant from Astronomical Almanac (ref. 7) 

e = obliquity of the ecliptic 
X - ecliptic longitude of the sun 

The computation proceeds as 

L 

G 

X 

R 

x 

y 
2 


v SC (1) + SC (2) d 

= sc (3) + SC (4) d 

= L + SC^ 5 j sin (G) + SCg sin (2G) 

- 1 - SC| 7 x cos (G) 

- R • A cos (x) 

= R • A - cos (e) sin X 
= R > A -in (e) sin X 
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where A is the semimajor axis of the earth's orbit In meters, 

The apparent diameter, AD, is found from the diameter at distance A and 

AD ■ (A/R) * DIAHCl) 

where DIAM is the diameter of the sun when A * R. 

It should be noted that this formulation does npt include any variation of 
the angular velocity of the sun in its computation. 

The position and diameter of the moon are computed using the daily polynomial 
coefficients for the moon. These coefficients form fifth-degree expressions 
for the right ascension, RA, declination, 6, and horizontal parallax, ir, of 
the moon in terms of the difference in ephemeris time (ET) between the epoch 
(0 hr ET of the date in question) and the simulation time. This time 
difference is found by adding the fraction of a day at which the simulation 
starts to the time since the start of the simulation and also summing UTET, 
the correction from Universal time to Ephemeris time. 

Once values of RA, 6, and it are found, they are converted to EC I Cartesian 
coordinates relative to the true equinox of date by 

R = Re/ it 

X = R cos <5 COS (RA) 
y = R cos 6 sin (RA) 
z = R sin (fi) 

AD = (A(2)/R) DIAM(2) 

where 

A(2) is the semimajor axis of the moon's orbit, in meters 
D I AM ( 2 ) is the apparent diameter of the moon when A - R. 


I 


INPUT 


Prog. 

symbol 

COMMON COMMON 
block block 
name location Dimen. 

Description (unit) 

TIME 

(calling argument) 


Time since start of simulation 
(sec) 

SL 

GRAVC 55 

6,3 

Constants for computing lunar position 
and horizontal parallax. The constants 
for right ascension are In positions 
SL(J, 1), for declination In SL(J, 2), 
and for horizontal parallax In SL(J, 3). 
The position 5L(1, 1) holds the a R 
coefficient for right ascension. 

SL 0 * 3) ■ 0, always ( = a,- for horizontal 
parallax) 3 

A 

GRAVC 76 

2 

Semimajor axis of solar then lunar orbit 
(m) 

I DATE 

GRAVC 87 

2 

Epoch of solar constants then date of 
simulation in form of year, month, day 
for a total of six digits 

FDATE 

GRAVC 89 


Fraction of a day elapsed from 0 hr UT 
at the start of the simulation. 

DIAM 

GRAVC 90 

2 

Apparent diameter of sun then moon 
when seen from distance A (deg) 

SC 

GRAVC 93 

7 

Solar constants available from the 
Astronomical Almanac. See equations 
(2 through 8) for use. 

OBLQ 

UTET 

GRAVC 100 
GRAVC 101 


Obliquity of the ecliptic (deg) 

Correction to add to UT to obtain ET 
(sec) 
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OUTPUT 


Prog. 

symbol 

COMMON 

block 

name 

COMMON 

block 

location 

Dimen. 

Description (unit) 

RSMKS 

(calling argument) 


Solar position, ECI (m) 

RMMKS 

(calling argument) 


Lunar position, ECI (m) 

GPOS 

GRAVC 

12 

6 

RSMKS then RMMKS 

R 

GRAVC 

74 

2 

D1 stance of sun then moon from center of 
earth In meters. 

UV 

GRAVC 

78 

3,3 

Unit vectors from center of earth sun 
(UV(1,1) to UV (3,1) and to moon (UV(1,2) 
to UV (3,2)) 

APDIAM 

GRAVC 102 

* 

3 

Apparent diameter of sun then moon 
as seen from the center of the earth (deg) 


GRAVC COMMON BLOCK LOCATIONS USED 


Program 

Symbol 

COMMON 

block 

Location 

Dimen. 

Description (units) 

TLAST 

11 


Store simulation time of last call to 
POSSUM (sec) 

GPOS 

12-18 

6 

Stores EC I position of sun then moon (m) 

SL 

55 

6,3 

Constants for computing right ascension, 
declination, and horizontal parallax 
of moon 

I FIRST 

73 


First pass flag 

R 

74 

2 

Distance of sun then moon from center 
of earth (m) 

A 

76 

2 

Semirnajor axes of orbits of sun then 
moon (m) 

UV 

78 

3,3 

Unit vectors from center of earth to 
sun then to moon 

I DATE 

87 

2 

Epoch of solar constants then date of 
simulation. Each is in six-digit year, 
month, day form. 

FDATE 

89 


Fraction of a day elapsed since 0 hr 
UT at the start of the simulator 

DIAM(2) 

90 


Apparent diameter of sun then moon 
as seen from distance A (see Toe, 76) 
(deg) 

SC 

93 

7 

Constants for computing position of sun 
(see eqs. (2-8)) 

OBLQ 

100 


Obliquity of the ecliptic (deg) 

UTET 

101 


Correction to add to UT to obtain JET (sec) 

APDIAM 

102 

3 

Apparent diameter of sun then moon as 
seen from the center of the earth (deg) 

COBL 

105 


COS (OBLQ) 

SOBL 

106 


SIN (OBLQ) 

Note: This 

routine is 

designed to 

be used with the Relative Attitude Pro- 


cessor (RAP). For this reason some of the COMMON block locations are 
left blank. They will be used by RAP and should not be eliminated. 
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QCONJG 

Subroutine QCONJG (Q» QCONJ) 

QCONJG takes the conjugate QCONJ of the quaternion Q. 


QNORM 

Subroutine QNORM (Q) 

QNORM yields a properly normalized quaternion with a positive scalar part. If 
the incoming scalar is negative, QNORM complements the vector part. QNORM com- 
putes the scalar part based on the vector part. Numerical problems in the 
vicinity of ISO 0 are avoided, and the resulting algorithm corresponds to an 
exact 180° rotation. 

QTM 

Subroutine QTM (Q, AMAT) , 

QTM generates a direction cosine matrix AMAT from its associated quaternion Q. 

QTOAX . * 

Subroutine QTOAX (Q, PSI, THETA, PHI, SINTH, COSTH, COSPHI, SINPHI) 

QTOAX extracts yaw, pitch, roll sequence Euler angles (PSI, THETA, PHI) from 
quaternion Q. For THETA = +90°, PHI is arbitrarily chosen to be 0°, PSI is 
then extracted so that (PSI, THETA, 0) represents the same attitude as the 
quaternion Q. The sine and cosine of THETA and PHI are also returned, The 
angles are returned in radians. 

QUTMLT 

Subroutine QUTMLT (PP, NP, QQ, NQ, R) 

QUTMLT generates the product R of two quaternions, PP and QQ. The order of the 
input quaternions is important, since in general 'quaternion multiplication is 
not commutative. If NP is negative, then the conjugate of PP Is used in the 
product. If NQ is negative, the conjugate of QQ is used in the product. 
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QXFQRM 

Subroutine QXFORM (U, QQ, N, V) 

QXFORM transforms the vector U into the vector V via the transformation qua ter 
nion QQ, thus V«QQ*U*QSONJ . If N is negative, the conjugate of QQ is used. 
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RELATIVE ATTITUDE PROCESSOR (RAP) 


PROGRAM DESCRIPTION 

The RAP models a designated optical system and associated field of view (FOV) 
and determines whether a specified object is within the optical FOV. The 
routine operates on a pair of systems (vehicles) as defined by the user with 
each system having one or more subsystems (subvehicles). Each primary vehicle 
may also have any number of defined optical systems with each optical system 
having its own coordinate system and two-dimensional FOV. Each subvehicle 
also has a defined coordinate system and a three-dimensional shape* 

Within the optical FGV, the object is located by an angular position with 
projected dimensions* The position of the observed vehicle or subvehicle is 
also computed in the Local Vertical -Local Horizgntal (LVLH) system. 


MATH MODEL 


The RAP is a flight software utility routine which can be either called or 
scheduled in the Space Shuttle Functional Simulator (SSFS) program. The 
routine uses data which are generated by the active vehicle models GRAV2 and 
POSSUM. Currently, RAP is configured for three vehicles although it could 
be extended for N vehicles and M subvehicles. RAP requires the following 
user input: 


NPV 

NSV • 
NSUBP 


NSUBS 


Designated primary vehicle with values of 1 for the Space Shuttle 
Orbiter (SSO), 2 for the MTV, and 3 for another designated vehicle. 

Designated secondary vehicle with the same description as NPV. 

The designated primary optical system located on the primary 
vehicle. For the SSO, the values are 2 for the rear window, 3 
for the overhead window, and 4 for the Flight Support Station 
(FSS). For the MTV, the values are 2 for the front camera and 3 
for the rear camera. 

The designated viewed system on the secondary vehicle. For the 
SSO, the values are 1 for the SSO vehicle, 2 for the rear window, 

3 for the overhead window, and 4 for the FSS. For the MTV, the 

values are 1 for the MTV body, 2 for the front camera, 3 for the 
rear camera. 







3-77 



' i . •, .vdfc^Hiirrt ihnaftitir? '• 'itfrtr imiar*. <■ : ^ 



For example, if NPV = 1, NSV * 2, NSUBP - 2, and NSUBS = 1, then the Shuttle 
rear window is the designated optical system and the MTV body is the viewed 
object. The block data associated with the designated input parameters is 
presented in Appendix A. 

After the systems have been defined, the position vector from the primary 
to secondary system in Earth-Centered Inertial (EC I) coordinates is computed. 

RVTVl = TtTV - W (1) 

where 

RTV is the ECI position vector of primary vehicle 
W is the ECI position vector for the secondary vehicle 

The resultant vector, RVTVl, is transformed from ECI to the primary vehicle 
system. This transformation and other transformations can be easily under- 
stood by referring to figure 3-8, which relates primary and secondary coordinate 
transformations. The vector RVTV, which is the primary- to-secondary system 
vector expressed in the primary system, is computed by 

RVTSr - (TIP) (RVTVl) 

★ 

where TIP is the inertial -to-primary transformation matrix. 

This relative position vector is adjusted for the optical offset and trans- 
formed to* the optical system by 

RVTVP = (TP0)(RVTV - RVOFF) (2) 

where 

RVOFF is the offset vector of the optical system as defined in the primary 
system 

TPO is the primary to optical system transformation 

Any offsets from the secondary to viewed system are computed with the final 
relative position vector in the optical system given by 

RVWP = RVTVP + (TSO) (RTWFF) 
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( 3 ) 



where 

RTWFF is the viewed system offset defining the .location of the designated 
subsystem (referred to the secondary coordinate system) 

★ . 

TSO is the secondary to optical system transformation matrix 

The vector RVTV’P defines the Line-of-Sight (LOS) vector from the optical to 
viewed systems (see figs, 3-9 and 3-10). The magnitude of the LOS vector is 

VMAG = V(RVTVP(1) 2 + RVTVP(2) Z + RVTVP(3) 2 (4) 

The dimensions of the viewed object are transformed to the optical system with 
the apparent dimensions in the optical FOV Y and Z axes computed as 

V ext = l(W 2)(T0 '' (2a)) l + l ( W 2)(TOV(2 ’ 2)) l + l<W 2)(T0V(2 > 3)) l (5) 
z ext = l< X dim /2)(TV0(3 ’ 1)) l * KW 2)(TOV(3>2)) l * I (z d1m /2)(TW(3>3)) I (6) 

where 

★ 

TVO is the viewed-to-optical transformation matrix 

Xdim, Ydim, and Zdiiu are the viewed vehicle dimensions in the viewed system 
The viewed vehicle apparent angular size as observed in the optical FOV is 



where 

AXS, AYS, and AZS are the apparent angular dimensions of the viewed vehicle 
RVEC is the vector from the viewed to optical system 

-v : a 

v Y 7 

B’ B’ B are defined viewed system unit axes 
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VEHICIE(SSO) COORDINATES 


TV CAMERA 



*gui" e 3 10,._ Coordinates of subvehicle in optical FOV 



The position of the center of the viewed vehicle in the optical system defined 
as the horizontal (Y-axis) and vertical (Z-axis) •angles is 

HANG « tan" 1 / R VTVP (2 ) / RVT VP ( ) ) ( 8 j 

VANG = tan -1 (RVTVP (3 )/RVTVP( 1 ) ) 
where RVTVP is the LOS vector in the optical system. 

To determine whether the viewed object is within the FOV, the following compu- 
tation is made: 

HANGS * [HANG] - Y t /| RVTVP | 

VANG2 = | VANG | - Z ext /| RVTVP | ^ 

Computed angles HANG2 and VANG2 define the apparent angular distance from 
the viewed subvehicle to the closest point on the Y and Z optical axes. If 
both of these angles are less than the optical FOV dimensions, the viewed 
vehicle is within the FOV. 


The routine RAP also determines whether the sun, moon, or earth is within the 
optical FOV by receiving ECI unit vectors describing the position of each 
object from the routine POSSUM which computes the current lunar and solar 
positions relative to ECI coordinates. POSSUM is called by the gravity 
environment model and updated a*, required by the integrator. These unit 
vectors are transformed to the optical system with the horizontal and vertical 
angles determined as previously discussed, except that the LOS vectors are 
replaced by the earth, sun, and moon unit vectors. 

The Y gxt and Z gxt for the earth, sun, and moon are simply the apparent angular 
sizes of these objects. With this information, the FOV check is made and the 
occulation of the sun or moon by the earth is considered. 

The position of the viewed object in the LVLH system (fig. 3-11) is computed 
using the vector RLVLH, the vector from the Shuttle to the viewed object in 
the LVLH coordinate frame. The LVLH angles as defined are 
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0 ® cos” (kOFLH * ZLV) 

<f> a tan (RlvlhY/RlvlhH) 
where ZLV Is the LVLH unit Z-axis defined by 

* LV “ 4 sso 

A 

where R ss0 is the unit EC I position vector of the Shuttle. 


The relative position and velocity of the Shuttle and the viewed secondary 
object are given by 

RELP = TT_ - IT. . 


VREL = V c 


where 


R"obj is the viewed object ECI position vector 

V sso is the ECI velocity of the Shuttle 

Y*obj 1S velocity of the viewed object, respectively. 

RAP OUTPUT 

The following variables output from RAP are defined as follows: 

HANG2 Leading edge of object if within FOV; horizontal angle in degrees 

HANG Center of object in FOV; horizontal angle in degrees 

VANG2 Leading edge of object if within. FOV; vertical angle in degrees 

VANG Center of object in FOV*, vertical angle in degrees 

HANGSM(l) Leading edge of sun*, horizontal angle in degrees 

VANGSM(l) Leading edge of sun; vertical angle in degrees 

HANGSM ( 2 ) Leading edge of moon; horizontal angle in degrees 

VANGSM(2) Leading edge of moon; vertical angle in degrees 

HANGSM (3) Leading edge of earth; horizontal angle in degrees 

VANGSM(3) Leading edge of earth; vertical angle in degrees 

pun w ) 
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THETLV 


Local vertical position angles of viewed vehicle in degrees 
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AXS 

AYS 

AZS 

RELP 

VREL 

RVTVP 


Apparent angular size of viewed object in degrees 

* 

ECI relative position vector between viewed and primary vehicles 
ECI relative velocity vector between viewed and primary vehicles 
Primary system relative position vector between viewed and 
primary vehicle 
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BLOCK DATA 

C0MM0N/FS1C/FSIU 

DIMENSION OATl I I 7 * • J) 

EQUIVALENCE ( F S f ISO) , OATl (1,1,1 )) 

EQUIVALENCE (FSfinO) NSUSP) ,?FS 1101 ) ,NSU0$> 

EQUIVALENCE iFSIlZfi) r N>VI r t r S(127)»NSVi 

DATA N®V /.V 
DATA M5V / 1 ✓ 

DATA N5UBP /2/ 

DATA NSUCS t\f 

LOCATION OF 0»I3i*N OF 5S0 IN SS p FRAME 

DATA |3ATHI,i,l >,1=1,3) 70,, D. *3./ 

MATRIX S 50 BODY TO SS1ISS0) 

DATA OATl ( 4 i If 1 ) / 1* :r 7 

DATA 0 A T l ( H • 1 * 1 ) 71. fl/ 

data da 11(15,1,1) /i.n/ 

DIMENSIONS OF SSO JN SSKSSP) 

DATA (DATltl, l,U, 1 = 13, 15) 7144. 7, 97.06,25. 52/ 

LOCATION OF OP 15 IN OF SSO SUBSYSTEM 2 REL TO SSO. 

data oath 1,2,1) /ssn. / 

DATA OAT If 2, 2,1) 7-. 057 
DATA 0AT1I 3,2, 1) /-1 17./ 

TRANS SSO BODY T 0 SSO SUBSYSTEM 2 (REAR WINDOW) 

DATA OATH 4,2,1) 7-1 »D/ 

DATA D A T 1( 8,2,1) 71*0/ 

DATA OATH 12,2, 1) /-l. 3/ 

DIMENSIONS OF SUBSYSTEM 2 IN THAT FRAME 

0 A T A 0 A T 1( IT, 2,1 > ,0A T 1 ( 1 4 , 2 , 1 ) , DA THIS , 2 , 1 ) / 1 . , H. , 1 . 7 

FOV OF SUBSYSTEM ? SYSTEM f -REAR WINDOW) 

DATA OATH 16,2,1 ) »OA T H 17,2,1 •/ 31 .,33./ 

LOCATION OF C° 15 IN Of SSO (OVERHEAD WINDOW) SUBSYSTEM 3 REL TO SSO 

DATA OATH 1 ,3,1) 7576. 2 / 

DATA OATH 2,3,1) /-,?•*/ 

DATA OATH 3,3.1) 7-122.4/ 

TRANS SSO BODY TO SSO SUBSYSTEM 2 (OVERHEAD WINDOW) 


f ) 


A-l 

R-88 


- - . i.Vi / *«. 




DAT* 9ini6,3,l) /I.O/ 

DAT* 3*1118.3. II n,n/ 

DATA OATmr,Sn /* 1 • 3/ 

DIMENSIONS of SUBSYSTEM 3 IN THAT FRAME 

DATA DA T H 17,3,1) ,DA T1 1 1 4 , 3 « 1 ) , DA Tl 1X5*3* 1 ) /!• , n . , l. / 

FOV OP nUB SYSTEM % S YSTEM 1 0 VERHEA C WINDOW) 

DATA DA T H I S, 3,1 ) ,DA Til IT* 3 * 1 I / *♦ 2 • » ** ll 

LOCATION OF ORIGIN OF SSO SUPSVSTFM 4 REL TO SSO 

flight gu s por t station 

DATA OATH 1,4,1) /-3 » 4 / 

DATA 0 A T 1 ( 1 1 4 , I ) /O.n/ 

DATA OA TH 3,4,1) /-M 1.3/ 

TPANS SSO BODY TO SSO SUBSYSTEM 4 (FSS) 

OATA OATH 6,4*1) / 1 • H / 

DATA OAT 1(3 ,4,1) 71.0/ 

DATA OA TH 10,4, 1 ) /-l. 3/ 

DIMENSIONS CF SUBSYSTEM 4 IN THAT FRAME 

D A T A 04 ’ H I 3 , 4 , 1) , D A T H 1 4 , 4 , 1 ) , D A T 1( 1 5 , 4 , ! ) / 9 . , 1 5 . , 4 . / 

LOCATION OF ORIGIN OT 4 T V SUBSYSTEM 1 REL TO MTV 

DATA DA T H 1,1,2) , DAT 1 ( 2 , 1 ,2 ) , DAT! (3 , 1 , 2 ) / O • ,0. , D» / 

TRANS MTV B °D Y TQ MTV SUBSYSTEM 1 


DATA OATH 4, 1, 2) /i.O/ 
DATA DAT 1(6, i f 2) /l.n/ 
DATA DA U| 10,1,?) /I .fV 


DIMENSIONS OF MTV IN MTV SYSTEM 

DATA OAT 1! 11, 1,7) ,DATH 14, 1,2), DA T1(1S, 1,3) / 3. 5, 2 .53,2 .58/ 

LOCATION Of ORIGIN OF MTV SS2 IN SSI FR AMF ( FORW AT D CAMERA) 

DATA (OATH 1,2,7) ,1 = 1.3) /2.,0. ,3./ 

MATRIX SSi TO SS 2 FRAME 

DATA OATH 4,2,2) /l."/ 

DATA DAT 1( fa 2,2) /l.n/ 

DATA OATH 12,2.2 ) /l.n/ 

DIMENSIONS OF SUBSYSTEM 2 IN 2 FR AMF 

OATA (DATH 1,2,?) ,1=17, 15) /D.,0. ,0./ 


A- 2 
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FOV OF MTV FOR WARP CAMERA 

DATA PATH 16,2,?)fDATl( 1 7, 2 ,2 >/ 43 . , 4GU/ 

LOCATION OF ORIGIN OF T V SS3 IN SSI FRAME! REAR DAMIRA ) 
DATA OATH I,E,?> ,1 = 1, J) /-2.,0.,0./ 


MATRIX 


;i TO SS3 FRAvr 


DATA DA TJ ( 4*3 ,2) /-! .0/ 


DATA OAT 1(6, 3.2) Zl.iT/ 
DATA DA T 1 ( 1 < , 5 , 3 ) /-l. 3/ 


DIMENSIONS OF SUBSYSTEM 3 IN 3 FRAME 


DATA OATHIt3»?)»I=13» 15) /D.,0i,0./ 
FOV OF MTV REAR CAMERA 


DATA OATH 16,3,? I, OAT If 17,3.21/16., 16./ 

LOCATION OF ORIGIN OF MTV? SUBSYSTEM l REL TO M TV 2 


DATA DAT Hip 1,3 ) , DAT 1 (2 , 1 , 3 ) , PAT1 (3,1,3) / 0. , 0. , 0. / 


TRANS "1 T V2 ROD Y TO «TV2 SUBSYSTEM 1 


DATA DAT1(4,1,3) /l.P/ 
DATA DAT 1(6, 1.3) /l.P/ 
DATA OATH 12, 1, 3) /I .TV 


DIMENSIONS OF MTV2 IN M T V2 SYSTEM 

DATA DA T 1 ( 1 3, 1 , 3 ) ,DA T 1( 14 , 1 , 3 ) , DA T1 ( 15 , 1 , 3 ) / 3. 5 ,2 .58 ,2 .58/ 
LOCATION OF ORIGIN OF -I T V 2f SS2 IN SSI FRAME ( FDRW A RD CAMERA) 
DATA OATH T, 2, 3) ,1 = 1,3) /2.,D.,3./ 

MATRIX ssi TO S S 2 FRAME 


DATA DA T 1(4,2, 3) /l.P/ 

DATA DA T l ( 6 , 2,3) /l.P/ 

DATA DA T V(1?,2*3) /,1 .0/ 

DIMENSIONS OF SUBSYSTEM 2 IN 2 FRAME 


DATA ( D A T 1 ( 1 , 2 , 3 ) * I - 1 3 , 1 5 ) / 0 . , 0 • ( 0 . / 

FOV OF MTV2 FORWARD CAMERA 

DATA DAT IT 16, 2, 3), PATH 17, 2, 3)/ 43. ,40./ 

LOCATION OF ORIGIN OF MTV2 SS 3 IN SSI FRAME (REAR CAMERA) 
DATA ( D A T 1 ( 1 , 3 , 3 1,1 -l » 3 > /• 2 • ,0 • , 0. / 
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MATRIX SSI TO SS 3 FR AMi 


OATA 

DATA 


DAT1M,3, 3) 
OATH 8*3*3) 


/-i.n/ 

/ 1 . 0 / 

DATA DAT 1(12, 3, 7) /-JO/ 
DIMENSIONS OF SUBSYSTEM 


DATA I D A T 1 ( 1,3,3), T 


. «'» 


IS) 


IN 3 FRAME 

/ 0 • , 0 • , O 1 • 


FOV OF M TV 2 PEAR CAMERA 

DA TA DA TIC 16 , 1, 3 > ,DAT1 (17,3,3 1/ IS . , 1 
END 


* • 0 S SUM/ SP 




RAP 

APPENDIX B 

PROGRAM LISTING AND COMMON BLOCK LOCATIONS 
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SGN 

Function SGN (WORD) 

SGN returns a 1 if WORD is positive, a 0 if WORD is 0, and a -1 if WORD is 
negative, 

SIG 

Function SIG (WORDT, W0RD2, NSDIG) 

SIG returns a 0 if W0RD1 is equal to W0RD2 to NSDIG significant digits, a 1 if 
W0RD1 is greater than W0RD2 in NSDIG significant digits, and a -1 otherwise, 

SIGNM 

Function SIGNM (WORD) 

SIGNM returns a 1 if WORD is positive or 0, and a -1 if WORD is negative. 








3.5 PLOTTER ROUTINES 


These routines process data from a SOAP simulation and cause plots to be 
generated. In general, they read and rearrange the data stored on logical 
unit 20 (and possibly 29) and create the calls to plotting utility routines 
which will do the actual point-by»point plotting. The plots which can be 
generated from a given simulation can include only that data specified upon 
preprocessing. Other Information will not be saved in the appropriate form. 
Typically, the user will neither modify nor directly call these routines. It 
Is only necessary to supply the proper plotting instruction cards (see section 
2) and to execute the available code. 

Routines described in this section arei 


AXIS 

FSPLOT 

NULINE 

PLMAIN 

PRPLOT 

SCALIE 

SYMBOI 

ENDFLE 

LINCNT 

NUMBER 

PLOT 

QUIKML 

SKPFLE 

TREAD 

FACTOR 

M INMAX 

PACK 

PLOTS 

SCALE 

SPMAIN 

WHERE 
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Subroutine AXIS (XP, YP, IBCD, +NCHAR, AXLEN, ANG, FIRSTV, DELTAV) 

$ 

AXIS is a Cal Comp subroutine which draws the graph axes* draws tic marks each 
Inch* writes the value at each tic mark, and then labels the axes. AXIS calls 
NUMBER, PLOT, and SYMBOL. 

XP, YP are the coordinates, In Inches , of the axis starting point. 

IBCD Is the literal title. 

+NCHAR defines number of characters and position of IBCD. 

AXLEN is axis length in inches. 

ANG Is number of degrees from X axis. 

FIRSTV Is the starting value of first tic mark. 

DELTAV Is number of data units/inch of axis. 

ENDFLE 

Subroutine ENDFLE (NLU) 

ENDFLE writes a coded end-of-file on a specified logical unit, NLU, 

FACTOR 

Subroutine FACTOR (FACT) 

FACTOR is an entry point in the Cal Comp plot package which allows the user to 
scale up or down the entire plot. FACT is the scale factor applied to the 
plot. A value of 2.0 doubles the plot size. A value of 1.0 Is the default 
value. • ■ • 


FSPLOT 


Subroutine FSPLOT {NLU, NPLOTS, NWDR) 

FSPLOT contains most of the logic of the production plotter, FSPLOT Is called 
by PLMAIN where NLU is the Input logical unit and NWDR is the number of plot 
variables, NPLOTS is a dummy argument, required In the argument list but 
currently not used, FSPLOT first skips to the correct data file* FSPLOT 
then reads in plot request cards until either 10 cards or 10 different 
parameters are obtained. These variables to be plotted are then read into 
core. A large loop is then started to produce the requested plots. PRPLOT 
is used for printer plots while QUIKML is used for S0-4060 plots. For Cal Comp 
plots, PLOTS, PLOT, SCALE (or SCALIE), AXIS, SYMBOL, and LINE are used, After 
all plot requests are processed, control is returned to PLMAIN, 

LINCNT 

Subroutine LINCNT sets the top and bottom page margins to nonstandard UN1VAC 
values for accommodating SOAP page plots. 

MINMAX 

Subroutine MINMAX (A1 , K, AMIN, AMAX) 

MINMAX defines the minimum, AMIN, and maximum, AMAX, of an array of K numbers 
beginning with Al. 

NULINE 

Subroutine NULINE (XA, YA, NPTS, INC, LN, IQ,T) 

NULINE, a Cal Comp routine, is called to plot the actual data, one line or plot 
at a time, using subroutines WHERE, SYMBOL, and PLOT. 

XA and YA are data arrays, 

NPTS is number of points. 
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INC is frequency of plotted data. 

IN describes type of line. 

* 

IQ. Is the Cal Comp Integer equivalent of the desired symbol. 

T is the array of angles by which symbols IQ are to be rotated. 

I ! 

NUMBER 

Subroutine NUMBER (XP, YP, HT, FPN, ANG, NDEC) ! I 

«. . t 

•y. 

NUMBER Is a Cal Comp routine which translates numbers to BCD and writes them j j 

on a graph (using SYMBOL). ; | 

. . . s 

XP and YP are coordinates in inches from lower left corner. ' 

HT is height In inches of the plotted number. J 

FPN is the floating point number j 

ANG is the angle from the X axis*. 

NDEC is number of decimal digits for output. 

PACK 

Subroutine PACK (WORD, CHAR, IPOS) 

Pack is called by PRPLOT to pack the BCD word WORD, where CHAR indicates the j 

direction of shift of IPOS characters. j 
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PLMAIN 

PLMAIN is the main program of the absolute element PLOTER used In execution 
of the production plotter. PLMAIN is a small driver program which uses 
FSPLOT to produce the desired plots. PLMAIN first reads a plotter control 
card to determine where the data is. If it is oh tape, the data is trans- 
• ferred to a high speed drum unit to speed up the data reading. FSPLOT Is 
then called for the actual plots. When control returns to PLMAIN, another 
control card is read if specified and the process described above is repeated. 
After all control cards are read, PLOT Is called to terminate the CalComp 
plot tape. 

PLOT 

Subroutine PLOT (XP, YP, +IP) 

PLOT is a CalComp routine which generates and writes instructions on tape to 
the CalComp hardware to produce the plots. Instructions to PLOT include (1) 
move the pen (with it down if IP = 2 or up if IP = 3) from the current 
position to the position specified by XP and YP; (2) end the plot if 
IP = -2 or -3 and move the pen to a new origin specified by XP and YP; 

(3) end the plot tape if IP = 999. 

PLOTS 

Subroutine PLOTS (WORK, NLOC, LU) 

PLOTS is an entry point in the PLOT subroutine called for initialization of 
PLOT. WORK is the storage array name, NLOC is the size of the WORK array and 
LU is the logical unit on which the plot tape is to be mounted. 
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PRPLOT 


Subroutine PRPLOT (AX» AY, K, IH, IC, ITITLE, IOM) 

PRPLOT is called by SPMAIN and FSPLOT to. produce printer plots. 

AX is the array of X values. 

AY Is a two-dimensional array containing Y values for all. curves in the frame. 
K is the number of points in each curve. 

IH is an array of characters for each curve. 

. IC is the number of curves in the plot. 

ITITLE is a title array, 

I DM is the number of points in AX. 

QUIKML / 

Subroutine QUIKML (+L* XMIN, XMAX, YMIN, YMAX , ISYM, BCDX, BCOY, NP, X, Y). 
QUIKML is a UNIVAC 1110 routine used to produce SD-4060 microfilm plots. 

L is the number of grids to be drawn on one frame. If L is negative, the 
frame will be advanced, 

XMIN, XMAX are minimum and maximum values of X array. 

YMIN* YMAX art minimum and maximum values of Y array. 

ISYM is the symbol to be used in plotting, 

BCDX, BCDY are alphanumeric titles for the X and Y axes, respectively.. 
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NP Is the number of plot points. 


X, Y are Initial locations of X and Y arrays. 

SCALE 

Subroutine SCALE (X, S, N, K) 

SCALE is a CalComp routine which defines a scale for a CalComp plot. SCALE 
defines the scale based on the desired plot size and the range of data such 
that is easy to interpolate between divisions on the resultant scale. 

X is the data array to be scanned. An adjusted minimum value will be stored 
in X (N*K + 1). An adjusted DX (maximum - minimum) will be stored In X 
(N* K + K + 1) 

S is the length over which data is to be plotted, 

N is the number of data points in the X array. 

« t . . : 

K is the repeat cycle of a mixed array. 

SCALIE 

Subroutine SCALIE (X, S, N, K) 

SCALIE is an alternate scale routine used with CalComp plots which produces 
larger plots than SCALE but less easy interpolation between divisions. 

X ( 1 ) is minimum value to be plotted. 

X(2) is maximum value to be plotted. 

X(3) is the output adjusted minimum value. 

X (4) is the adjusted delta value per inch. 

S is the maximum length of plot in inches. 


N and K are dummy arguments. 


* , * 


SKPFLE 

Subroutine SKPFLE (NLU) 

SKPFLE is used to skip a data file on a specified logical unit* NLU. 

SPMAIN 

SPMAIN is the main program of the absolute element SPDPLT used In execution 
of the high speed plotter. To reduce execution time, SPMAIN plots data from . 
core instead of continually reading a data tape. To enable this, SPMAIN 
restricts data tapes to 30 variables with 1350 points per variable per file.’ 
SPMAIN ignores additional data or parameters if either limitation is exceeded. 
SPMAIN first reads a control card to find out where the data is and then brings 
the data into core. If CalComp plots have been specified, the maximum and 
minimum for each variable is determined* assuming that all will be plotted. A 
large loop is then entered to produce the requested plots. For CalComp plots, 
PLOTS, PLOT, SCALE (or SCALIE), AXIS, SYMBOL, and LINE are used. For printer 
plots, the data to be plotted is reduced to 100 points or less corresponding 
to the available page field width. PRPLOT is used for printer plots, if 
requested. 

¥ » 

SYMBOL 

Subroutine SYMBOL (X, Y, HT, IBCD, ANG, N) 

SYMBOL is a CalComp routine which will plot N alphanumeric characters, IBCD, 
at a height of HI inches and angle of ANG degrees beginning at coordinates 
specified by X and Y, 
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TREAD 

# « 

Subroutine TREAD (NLU, NW, MW, IP, IK, BUFF) 

TREAD is used to read 1 a data record, BUFF, from the input logical unit, NLU. 
NW and MW are the first and last indices of the BUFF array, 

IP is a parity indicator, 

IK is a read operation indicator. 

WHERE 

Subroutine WHERE (X, Y, F) 

WHERE is an entry point in the Cal Comp routine PLOT, WHERE is called by LINE 
to define the current CalComp pen position where X and Y are plot coordinates 
and Fisa scale factor. 



4. ENVIRONMENT MODELS 


The Environment models serve to calculate the state vector (position, 
attitude, and rate of change of each) of each vehicle. They support the 
computation of forces and torques, the updating of mass properties, and 
the integration of the equation of motion. Environment models have a 
special position within the SOAP system since they allow input and output 
through individual COMMON blocks with names which are known to the SOAP 
system, since they perform their computations in a monolithic block so 
that all parameters are updated at once, and since model specifications 
(typically what number of vehicles) are made at the time of preprocessing 
and at the time of mapping. Typically, Environment models are not scheduled. 

It should be noted that default values, which occur in lieu of specification 
of the Environment model at the preprocessing and mapping stages, will 
result in the use of a dummy model and the function of that model will be 
omitted. 


4.1 ACTIVE VEHICLE 


PROGRAM NAME 

Single precision ACTVEH 
Double precision AVEH/AVEH50 

PROGRAM DESCRIPTION 

ACTVEH is a driver which calls active vehicle models for each free flier. 

These models update the position and attitude of each vehicle. ACTVEH resets 
the flag IPDYN in ENVCOM prior to calling each active vehicle. This flag 
informs the vehicle models whether to update the dynamics derivatives (linear 
and angular accelerations), whether to update the active vehicle state, and in 
Which order to perform these actions. ACTVEH updates the active vehicle 
internal times (AVTIM) to the update time required. 

ACTVEH has no calling arguments and accepts no commands. 


PR OGRAM NAME 

AVES24 
AVES19 

AVED19 
AVED24 
AVED242 

The single- and double-precision versions of this model differ in the names of 
some of the subroutines called but are otherwise the same in operation and 
formulation. 

PROGRAM DESCRIPTION 

The active vehicle programs provide six-degree-of freedom simulations of a 
rigid body vehicle in three-dimensional inertial space. The vehicle dynamics 
use the forces torques caused by firing jets of an Attitude Control 


Single precision 
Double precision 
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I 


* 

* * . 

Propulsion System (ACPS), which is modeled by the Reaction Control System 
(RCS). The gravitational acceleration of earth (GRAV model), the torque due 
to the gravity gradient (GGT model) and the aerodynamic forces and torques 

( (AERO model) are used in the vehicle dynamics. Computation of the vehicle's 

* 

state vectors 'is performed relative to a nonrotating orthogonal coordinate 
system with the origin at the center of the earth (earth-centered inertial 
(EC I ) ) . The rotational sequence used is roll, pitch, yaw. Each model is 
supported by a COMMON block (e.g., ACTVEC, ACTV1C, ACTV2C) and a block data 
routine which initializes the COMMON block. 


MATH MODEL 

The AVES24 program maintains the vehicle state vectors in ECI coordinates. 
The contact forces and torques are in the vehicle body frame. By summing 
these forces and torques, the contact acceleration of the vehicle is 
calculated. 


4v = + ^AERO 

(1) 

^AV = ^KCS + ^AERO 

(2) 

^C = ~av /m av 

(3) 


where 

F - Forces 
T = Torque 
M = Mass 

A V = Active vehicle 
C = Contact 
A = Acceleration 



The torque due to the rotational rate of the vehicle, T , is calculated and 

— UlV 

is combined with the torques due to the contact forces and with the torques 
produced by the gravity gradient, T^q, calculate the vehicle angular 
acceleration. 


T 

— <i) v 


U> 

-v 


(I “v) 
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(4) 


where 


Wy is the vehicle angular velocity 
I is the active I vehicle inertia tensor 

The total acceleration (Ay) of the vehicle due to the contact forces and 
torques and the vehicle inertia is calculated. 

Aj = 4 + Ry X >< “y + % >< (iy X ^ My) (6) 

The total acceleration of the body is transformed from the vehicle body frame to 
the inertial frame and then summed with the gravitational accelerations, Au, 

v3 

iiiv = \ — t + As 

where Iv denotes inertial vehicle. The matrix J M y transforms from vehicle to 
inertial reference. 

A fourth-order Runge Kutta integrator (DIFEQS for single precision or DPDFQS 
for double precision) is used to integrate to obtain £ and £ in ECI 
coordinates. The time step used is either the active vehicle time step AVSTEP 
(from ENVCOM), the difference in time required to bring the active vehicle 
state up to the time required by the rest of the simulation, or TSTEP (passed 
in ENVCOM), whichever is smaller. 

The updated vehicle body rates , are computed by integrating <Ly using a 
fourth-order Runge-Kutta integrator (DIFEQ for single precision or OPDFQ for 
double precision). These updated body rates are than used to find the 
derivative (DAMVI) of the vehicle-to-inertial attitude matrix (AMVI) via 



0 

" u v,3 

“v,2 

DAMVI = 

w v,3 

0 

■ u v,l 


--w v ,2 

“v,l 

0 . 


DAMVI 


AMVI 


( 7 ) 


The updated veh1cle-tO"1nertial-attitude matrix Is computed by again 
Integrating, using DIFEQ or DPOFQ. 

Once the updated attitude matrix has been computed, the Euler angles are found 
using the Euler angle derivatives (obtained by a call to EllLERD or OEULRD) to 
compute a change In Euler angles (by a call to EULERE or DEULRE), The values 
of the Euler angles are then updated by adding the changes computed. This 
procedure requires a call to EXTRCT or OXTRCT. 

The integration is repeated until the time associated with the vehicle state 
equals the time requested. 

The active vehicle models have two initialization passes during which the 
inertial -to-vehicle attitude is initialized, its transpose is computed, and 
initial values of the environmental forces and torques are computed. 

The active vehicle models have no calling arguments and accept no commands. 
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4.2 AERO 
AER0/AER050 

AERO is the driver routine which calls models of the aerodynamic force and 
torque for each vehicle in the simulation sequentially. It has no calling 
arguments and no COMMON block requirements. It is specified on the environ- 
ment input cards (see section 2). 


AER0/AER037 
PROGRAM DESCRIPTION 

AER037 models the aerodynamic forces and torques on an on-orbit vehicle. The 
vehicle is approximated by three flat plates. In the case of the Space Shuttle 
Orblter (SSO), the aerodynamic characteristics modeled are for the 147/140B 
approximation to the 140 A/B Space Shuttle Vehicle Orblter (ref, 8), This 
model includes a calculation to approximate the atmospheric density for the 
1962 U.S. Standard Atmosphere (ref. 9). 


MATH MODEL 

AER037 approximates an on-orbit aero model for the orbiter using the Newtonian 
Impact Theory for a flat plate. For the orbit altitude range being considered, 
the drag coefficient, C D , is assumed constant. Further, at this altitude, the 
lift coefficient, C^, is zero. Therefore, torque on a vehicle can be approxi- 
mated by modeling the vehicle as a combination of three flat plates and de- 
termined from the drag forces acting on each plate (ref. 10). 


Calculation of the relative atmospheric density is dependent upon the geometric 
altitude above mean sea level. For geometric altitudes greater than 500,000 
feet, the equation is 


P f = P 0 * exp j - [24.2 + 0.002889Z - 2605. /Z] j 
or for geometric altitudes from 280,000 to 500,000 feet. 


II 


Pf * P 0 * exp 


- [39.57 - 0.01 044 Z - 6956. /Z] 


where 

is the relative atmospheric density in pounds per cubic foot 

P 0 Is the atmospheric density at the earth's surface In pounds per cubic 
foot 

Z Is the geometric altitude In kllofeet 

Then to convert the density to metric units, p * 16.01843 p^ kilograms per 
cubic meter. 

The dynamic pressure, q, Is defined by q « 0.5 pV* newtons per square meter 
where Is the magnitude of the inertial velocity In meters per second. 

The aerodynamic forces, Oj, are calculated by 

5|' * C D d Al i k * T 0 i? 0 newtons 

, 0 2 - C D q A s |j • T 0 |T 0 newtons 

h ’ c d q ’ V l 0 newtons 

where 

A| Is a vector normal to the true planform area having a magnitude (Aj) 
equal to the true planform area 

Cq Is- the drag coefficient 

I Q Is the unit vector of the inertial velocity vector 

The vehicle axes with the origin at the vehicle center of mass are defined 
as ■ 

Xy - Out the nose 

t 

Yy - Out the right wing 

Zy - Down 
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The vector Aj Is normal to an X-Y plane plate and in the direction of Z. The 
A g is normal to the X-Z plane plate and in the direction of Y. The A 3 Is 
normal to a Y-Z plane plate and In the direction of X. 

The origins are the centers of pressure of A y A 2 , and A 3 * respectively. 

The torques, Hi* are calculated 

* * 

Mj ■ r^ x Bj newtons -meters 
H 2 * F 2 x D 2 newtons-meters 
M 3 * r 3 x D 3 newtons -meters 

Where rj is a position vector from the center of mass to the center of 
pressure of A i in meters. 

The total aero forces, Dj, and torques, My, (ignoring shading of one part of 
the vehicle by another) is 

Dy » E D.J » Z C Q q | (3V| • T 0 )|T 0 newtons 
My * Z r^ x newton-meters 


INPUT/OUTPUT 

This model has no calling arguments. 

The following data input is required: 

i 

Locations of the center of pressure in the Fabrication frame. (Rp^) 
Drag coefficient (Cq) 

Magnitudes of the true planform areas (A^ , A 2 , A 3 ) 

All the above data can be changed using the SOAP data input scheme. 
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Input from the active vehtcle and mass properties common blocks Is required, 
AER037 requires on each update from the active vehicle the following: 


Inertial position vector (R JV ) 
Inertial velocity vector (R JV ) 


Vehicle to Inertial attitude transformation matrix (1 m ) 


Only on Initialization of AER037 Is Information required from the mass prop- 
erties common block. That Is 


Fabrication to vehicle position vector (Rpy ) 

Fabrication to vehicle attitude transformation matrix (^p) 


The following list gives the program mnemonic, location, and stored value of 
these variables. 



Variable 

Symbol 

Location 

Description 

Value 

Units 

r 

r fpc 1 

AERSTA(l) 

7 

Station locations of the 

1111.5 

inches 

AERSTA 2 

8 

center of pressure of A1 

0.0 


i i . 


AERSTAC3) 

9 

363.0 


\ 

r fpc 2 

AERSTA (4) 

10 

Station locations of the 

1064.4 

inches 

\ 

1 

AERSTA 5 

11 

center of pressure of A2 

0.0 



AERSTA (6) 

12 

414.4 


* 

Rpn/' 

AERSTA (7) 

13 

Station locations of the 

1300.0 

inches 


fpc 3 

AERSTA (8) 

14 

center of pressure of A3 

,1 

0.0 




AERSTA (9) 

15 

406.3 


I : 

CD 

CD 

16 

Coefficient of drag 

^.0 


*■ 

A1 

A1 

17 

Magnitudes of the true 

397.76 

meters 

\ 

f 

A2 

A2 

18 

planform areas 

280.95 


i • ; 

A3 

A3 

19 


99.81 
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This Information Is used to compute ty r 2 , and r 3 positions of the center of 
pressure. 





- Rpy) meters 





Rp V ) meters 


r 3 * V Mp ^ R FCP 3 " meters 

Output from AER037 Is 

Sum of the aero forces (Fy) In newtons 

Sum of the aero torques (My) in* newtons -meters 


Output of these forces and torques Is optional and Is controlled by an aero 
print; control switch, AEROPC. When AEROPC is 0.0, AER037 will not print. 

When AEROPC is 1.0, the output parameters are printed at the specified delta 
time contained in AERPDT. When AEROPC Is 5.0, the output parameters are out- 
put every time AER037 is updated. Starting the print at some time later In 
the run can be accomplished by setting the desired start time In AEROPT. All 
print Is nominally turned off. 

STORED VALUES 

AER037 defines Initial values for several parameters located In AEROC common, 
which may be altered via input cards. 

SOURCE 

AER037 was developed from AER018. .Improvements in AER037 Include calculations 
for three plates Instead of two as described In reference 9, and calculation j 
of the atmospheric density as defined In reference 8. 1 
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VERSION AVAILABILITY 


Both single and double-precision versions of AER037 are available, Block Data 
routines (BDAERO) are also available to support them, 


PROGRAM VERIFICATION 

The AER037 model subroutine was verified In closed-loop SSFS simulation runs. 
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AERl 

PROGRAM DESCRIPTION 

AERl computes forces and torques on a vehicle on orbit, using the same math 
model employed by AER037; the only differences are In the Initialization pass 
and the Input data, AERl bypasses use of a fabrication frame, therefore plate 
positions, normal vectors, and areas are Input directly In vehicle body coordl 
nates. All Input and output ‘Is Identical with that for AER037 with the excep- 
tion that Rpy and are not used. 

VERSIONS 

Both single and double-precision versions are available. 
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4.3 GRAVITY 
GRAV2 


PROGRAM DESCRIPTION 

GRAV2 provides a gravity model for earth orbital operations. 
Various degrees of fidelity in gravity calculation can be 
selected. For the lowest degree of fidelity, GRAV2 returns 
the spherical earth gravity; whereas for the highest, GRAV2 
adds perturbations caused by the sun, the moon, and the 
nonspherical nature of the earth. The user can select any 
combination of the perturbations caused by the sun, the moon, 
and the nonspherical nature of the earth to be added to the 
spherical earth gravity (for example: nonspherical only). 

MATH MODEL 

Spherical gravity is calculated using 

... i i • • • .... 

-u K 

S - — — - 

; ** m 3 

where 

E is the spheri.cal earth gravitational acceleration. 

®P ... 

is the gravitational constant of the earth. 

IT is the position vector of the point at which 
the gravitational acceleration is desired, in 
earth- centered- inertial (EC!) coordinates. 
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The nonspherical earth perturbation includes up to the . 
fourth harmonic of the earth. ‘The equation used is: 

, n * HS « ZS(7X • 3) ♦ FK(6x - t - 9x*j 

w (L 1 w . 

♦ T-2FJ • ZS * HS (0 .6 - 3k) » fK * ' b -- j ^~| 

where 

^N-s * s t ^ ie n o ns P^erical gravity perturbation 

is the gravitational constant of the earth 

K is the vehicle position vector in earth-centered 
inertial (ECI) coordinates 

FJ is a constant relating to the nonsphericity of 
the earth 

» i 

HS is a constant relating to the nonsphericity of 
the earth 

FK is a constant relating to the nonsphericity of 
the earth 

ZS is the z-component of the position vector E 

a is' a constant equal to 0.4235714286 

b is a constant equal to 1 . 714285714 

x is ZS 2 /|E| 2 

I is a unit vector in the Z (R(3) ) direction 


The sun perturbation acceleration is calculated using 


777? • ? .v> 

I r v* I 


where 


v* 


es 


is the sun's gravitational perturbation 
acceleration 

is the gravitational constant of the sun 
is the vector from the vehicle to the sun 
is the vector from the earth to the sun 
is the vector from the earth to the vehicl< 


•v 

and q is defined: 


X 


11 i m i. X>1 


[1 ♦ X][(l * X) z + /TTT] 


where 







‘ * ‘ 

The moon perturbation is the same as above with 
scripts changed to m and the position Vectors 
the moon's position instead of the sun's. 


The units of the symbols in the preceding equat 
Shown below. 


s sub- 
reflect.ing 

cns are 







F 

m 

a 

m 


m/sec 2 

b * 

m 

•p 




M-S 

m/sec 2 

X 

m 


m/sec 2 

r 

m 


m/sec 2 

7 v. . 

m 


m 3 /sec 2 


m 

M. 

m 3 /sec 2 

F .v 

TO 

^n» 

m 3 /sec 2 

*vm 

m 

FJ 

m 2 

r «m 

m 

HS 

m 3 

q 

- 

FK 

m 4 

X 

» 

ZS 

m 




A complete discussion of the mathematical model appears in 
reference 11 , 

All measurements are made in the earth-centered* inertial 
(ECI) coordinate system, 

INPUT/OUTPUT 

GRAV2 is called with three calling arguments: 

GTIME, FOS, tjRVTY 

where 

GTIME is the current time, 

FOS is the position at which the gravitational 
acceleration is to bo calculated, 

SRVTV is the gravitational acceleration vector. 
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GTIME and 70S are input variables and (TRVTY- is output. 

JPERT , a flag located in cell number one of GRAVC COMMON, 
is used to determine which perturbations (if any) are to be 
added to the spherical earth gravity. JPERT is an octal 
integer and can take on any value from 0 to 7, Each of the 
three least significant bits of JPERT is a flag indicating 
whether a particular perturbation is to be included in the 
gravity computations as follows: 

1st (least significant) bit Nonspherical perturbation flag 
2nd bit Sun perturbation flag 

3rd bit Moon perturbation flag 

The quantities reflected in the calculation of SRVTY , 


depending upon 

the value 

of JPERT, 

are shown below: 

JPERT 

3rd 

Bit 

(MOON) 

2nd 
* Bit 
(SUN) 

1st 

Bit 

(NONSPH) 

&RVTY Reflects 

0 

0 

’ 0 

0 

SPHER only 

1 

0 

0 

'1 

SPHER + NONSPH 

* 

2 

0 

1 ‘ 

0 

SPHER + SUN 

3 

0 

1 

1 

SPHER + NONSPH + SUN , 

4 

1 

0 

0 

SPHER + MOON 

5 

1 

0 

1 * 

SPHER + NONSPH + MOON 

6 

1 

1 

0 

SPHER + SUN + MOON 

7 

1 

1 

* 

X 

SPHER + NONSPH + SUN + MOON 


where SPHER stands for "spherical” and NONSPH stands for 
"nonspherical." 


GRAV2 initially defines the value of JPERT equal to zero. 

* * 

This value may be changed by input if desired. On the 
initialisation pass GRAV2 extracts the perturbation flags 
from JPERT and stores them in GftAVC COMMON as follows: 

GRAVC(2) Nonspherical perturbation flag 
GRAVC(3) Sun perturbation flag 
GRAVC(4) * Moon perturbation flag 

GRAV2 does not alter the contents of cells GRAVC(2) - (4) after 
initialization. Thus, if for any reason, the user desires to 
alter JPERT after initialization, he must also alter locations 
two through four of GRAVC to reflect the new value of JPERT. 

GRAV2 defines initial values (via DATA statements) for the 
following common block variables: 

JPERT FJ 

HS 

U, FK 

>*» • • 

The initial values will appear in the COMMON Block Map 
section of the SSFS User's Guide. They may be changed by 
input if desired. 

<► ’ 

■ * ■ 

GRAV2 calls one additional subroutine ,' POSSUM, to calculate 

’ t. 

the position of the sun and moon when perturbations due to 
those bodies are requested. Several versions of POSSUM exist. All of them 
use tabulated data. Typically, they use positions 55 and following in GRAV 
COMMON. 


STORED VALUES 

GRAV2 defines initial values for two local variables, a 
and b , used in the calculation of the gravity perturba- 
tion dw to the nonsphcrical nature of the earth. 

a - 0.4285714286 

b - 1.714285714 

COORDINATE SYSTEMS 

GRAV2 uses the earth-centered- inertial (ECI) coordinate 
system. 

VERSIONS AVAILABLE 

Single and double-precision versions of 6RAV2 are available. 

PROGRAM VERIFICATION 

~ ' ITm ' r ~ ' » » 

Open-loop checkout of GRAV2 consisted of comparison runs 
between GRAV2 and the gravity portion of the UNIVRS sub- 
routine used on the AGC bit-by-bit simulator. Test cases 
were run using the following inputs: 


Test 

Case 

PERT 

(for 

UNIVRS) 

JPERT 

(fOT 

GRAV2) 

Time 

(Sec) 

Position 

Vector (meters) 

1 

0 . 

0 

SO. 

6.4E+06, 

7.SE-06, -7.0E+06 

2 

1. 

7 

SO . 

6.4E-06, 

7.SE-06, -7.0E-06 

3 

0 . 

0 

ISO. 

-7.0E+06, 

5 . 0E-07 , 6.0E-06 

4 

,i. 

7 

ISO. 

-7.0E-06, 

S.0E+07, 6.0E-06 

S 

0 . 

0 

1000. 

8.0E+06, 

0 .0 V 0.0 

6 

1. 

7 

1000. 

8.0E-06 , 

0.0, 0.0 


The results of the comparison are listed below. 


Test 

Case 

* * 

Subroutine 

* 

Gravitational Acceleration C meters /sec 3 ) 

1 

GRAV2 

UNIVRS 

-1 .4429S56E+00 , *1 .6909636E+00 , 1 ,S782327E*00 

•1 ,44295365*00 , -1,39096365+00, l.S7S2326E*00 

2 

CRAV2 

UNIVRS 

-i.4423145E*00, -l.6904477E*00, 1.37917915+00 

-1.442514SE+00, -1.69044775*00, l. 379179 IE+C0 

3 

GRAV2 

UNIVRS 

2 , 1 2 2 9 79 7E * 0 2 , *1.316 4141E -01, -1,31969 69 E-02 
2.1229796E-02,’ -1.S164140E-01, *1,319696*5-02 

4 

GRAV2 . 
UNIVRS 

2.1231229E-02, -1 . 3164S34E-01 , -1.3196672E-02 
2.123122S5-02, -1 . S164533E-01 , -1.3196671E-02 

5 

GRAV2 

UNIVRS 

•6.2281439E+00, 0.0, 0.0 

•6 , 2281438E+00 , 0.0, 0.0 

6 

GRAV2 

UNIVRS 

*6 , 2 3 4 S 7 3 5 E+ 0 0 , 1.2922539E-07, *1. 06 9509 75 * 03 
-6.23437335*00 , 1 . 2922641E-0? , -1.0S9S097E-03 


The al»<?ve comparison verifies the correctness of the GRAV2 
gravity model, 

SOURCE * 

The GRAV2 math model has been extracted from the UNIVRS 
subroutine of the AGC bit-by-bit simulator. 
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For locations 55 and subsequent, see POSSUM 




4.4 GRAVITY GRADIENT 
GGT1 

PROGRAM DESCRIPTION 

•» 

GGT1 models the generalized gravity gradient torque, which 
is produced about the vehicle's center of mass by the earth 
gravitational acceleration acting on heterogeneous and 
asymmetrical vehicles. A spherical earth is used since the 
torque produced due to the earth's oblateness is negligible 


MATH MODEL . 

The gravity gradient torque components (T ,T ,T ) about the 

i w .x y Z» 

vehicle reference axes are 
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where u * earth's gravitational constant [PS:EGC] 

E » vehicle inertial position vector expressed 
in the vehicle body frame [PS:R] 

I * elements of vehicle inertia tensor in the 

vehicle center of gravity originated system. 
[PS: VI] 

D m direction cosines between the radius vector 
R and the vehicle body axes. [PS;D] 


The equations above use the negative products of inertia. 
All calculations, are made in the vehicle body frame. In' 
the above [PS: ] indicates program symbol. 
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INPUT/OUTPUT 

GGT1 has two calling arguments. * , , 

K ■ inertial state vector expressed in vehicle body 
coordinates, (meters) 

I - vehicle inertia tensor (kg-m ) 

» 

One additional parameter is required for input and is loca- 
ted in locations 7 through 12 of the GGT1. 

JOT ■ gravity gradient torque axial multiplies (one 
for each axis) , which are used to scale or * 
turn off the torque about a particular axis* 

T 

(nominally set to 1) 

*’ * 

Output from GGT1 is 

t 

• T * gravity gradient torque [PS:GGT] (M-NT) about 
1 the vehicle center of gravity (M-Nt). [PS:GGT] 

VERSIONS AVAILABLE 

Both single and double-precision versions are available. 

PROGRAM VERIFICATION 


The verification of GGT1 consisted of running several test 
cases and making hand calculations. One test case was run 
for a complete circular orbit to illustrate the values 
obtained from GGT1. In all cases, the same inertia tensor 
was used. 
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2725000.0 

•20000.0 

-333000.0 

I • 




m 

•20000.0 

19883000.0 

» 

SOOO.O 


K.. 


I *». 


•333000.0 

• 

5000.0 

21495000.0 


A spherical gravity was used and the initial inertial state 
was defined using 


• R JV * 6600000.0, 0.0, 0.0 meters 

V jv ■ 0.0, 7770.0, 0.0 meters/second 

♦ it 

In this case, the GGT1 routine calculated the values in 
table 4-1 for the circular orbit of 5400 seconds. The GG 
Torque plot samples the Torque vector every 20 seconds. 


l 


# 
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* 
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TABLE 4 

* 

-I.- GRAVITY 

GRADIENT TORQUE 

VECTOR 

i *■ 

Run 

GCT * 

GGTy . 

GGt 

Time (Sec) 

(NT-M) 

(NT-M) 

(NT-M) 

• 

0 

- ,06998861 

44.89743710 

- .04040712 

600 

3.11959702 

29 • 31005049 

-39.59090602 

1200 

-1.31008284 

-45.24714470 

-24,70406365 

l 800 

.17021475 

-22.90914392 

2S. 18761992 

2400 

-1.40059164 

5. 56837219 

47.49012947 

,?000 

1.96025780 

-27.64657974 

» si * * 

39.47186518 

4 - « 

. 3600 

,05934277 

- 3.77952236 

- 3.82524449 

4200 

( - .84780970 

49.77582121 

5.13785291 

4800 

-3.68393409 

22.61737800 

35.36430979 

S400 

.43922012 

-29,66196251 

27.02822924 
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4.6 HORIZON SENSOR 
HORIZ 


> I f 


Subroutine name: H0RIZ1 , H0RIZ2 
PROGRAM DESCRIPTION 

These routines model the performance of a Barnes Model 13-166-9. This apparatus 
Is pictured In the figure below. It senses pitch and roll at a frequency of 
•12,5 Hz. If the vehicle Is either pitched or rolled more than 5° from local 
vertical, she sensor output Is 5°. No other inaccuracies are modeled. Yaw 
Is not sensed. 



Figure 4-2.- Horizon scan configuration, 


MATH MODEL 

The horizon sensor model obtains a unit vector antiparallel to the local 
vertical, ZLVM, by transforming RIV, the position vector of the vehicle in 
inertial coordinates, to body coordinates via the transformation AMIV 


ZLVM = AMIV •Unit (RIV) 
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The dot product of the body X vector with the local vertical vector gives 
the cosine of the angle between the two. If pitch equals zero, this angle 
should be 90°. For negative pitch, the angle Is greater than zero, so the 
pitch In degrees (PHS) Is found to be 

PHS « 90°- ^ . cos* 1 (ZLVM(l)) (2) 

Similarly, for roll (RHS) 

RHS * 90°- yjjjy • cos * (ZLM(2)) (3) 

If either pitch or roll exceeds 5° In magnitude, It Is limited to 5° and the 
sign Is retained. 

This model has no calling arguments. 

INPUT/OUTPUT 

The symbol # denotes a vehicle number. 


Prog . 
symbol 

I/O 

COMMON 

block 

name 

COMMON 

block 

location 

D / i men. 

Description 

(units) 

RIV 

I 

ACTViPC 

1 

3 

Inertial to vehicle 
position vector (m) 

AMIV 

I 

ACTVIC 

100 

9 

Inertial to vehicle 
attitude transformation 
matrix 

LCBHS# 

I* 

LECCOM 
# = 1 
# = 2 

1348 

1406 

1 

Horizon sensor COMMON 
block length 

IFGHS# 

l* 

me 

70 

1 

Horizon sensor on flag 
1 = on 


2 = off 

♦Input required to use this model. 
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COMMON 

block 

name 

COMMON 

block 

location 

Dimen. 

Description 

(units) 

HSIC 

1 

1 

Roll angle (deg) 

HSIC 

2 

1 

Pitch angle (deg) 

HSIC 

3 

3 

Unit vector from origin 
of EC I coordinates to 
vehicle In the vehicle 
body coordinate system. 
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MASPRO 

MASPRO Is the mass properties driver, It calls mass properties routines MASPRO 
MASPR1 ... , etc, It maintains no dock and has no COMMON block. There are 
no calling arguments, 

MASPRO 

PROGRAM DESCRIPTION 

MASPRO defines the mass properties of the Orbiter with an aft payload for 
Operation Mission 2 (sortie). Values are chosen at pre-OMS deorbit burn as 
in reference 1. The center of mass and vehicle inertia are held constant 
while the mass can be decremented to model fuel consumption by the RCS 
jets and main propulsion engines* The nominal mass properties represent the 
Orbiter with a 32,000 pound payload. An optional mass properties data set 
may be selected to represent the Orbiter with no payload, 

MATH MODEL 

Except for the change in mass, MASPRO performs all its calculations during 
the initialization pass only. 

* 

The structural characteristics of the Orbiter vehicle are defined in the 
fabrication (F) coordinate frame. The location of the vehicle's center of 
mass (X) frame and the body (V) frame are expressed In the fabrication frame, 
This mass properties model calculates the position of the vehicle's center 
of mass In the body frame. 

Initial parameters are input to MASPRO via data statements in English units 
and converted internally to the MKS system. 

The attitude of the vehicle (V) body frame in the fabrication (F) frame, 

Apy ■ (0.0, 180.0, 0.0) degrees (1) 

is converted to radians to define the fabrication to vehicle body frame 

transformation matrix, V u . 

M F 

. * ■ 
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The vector ^ Is the position vector, expressed In the vehicle system, 
from the origin of the vehicle coordinate system to the actual center of 
gravity. The equation used to define W vx is 

j. 

V 

• tv) (2) 

where 

T?P X Is the position " the vehicle's center of mass (X) In the fabrication 
(F) coordinates. 

and 

Rpy Is the position of the vehicle's (V) body frame in the fabrication (F) 
coordinates. 

The inertia input is In fabrication coordinates and the product elements are 
defined by positive integrals. However, the irfertia tensor as used by the 
vehicle routine is referenced to vehicle coordinates and the product ele- 
ments are defined as negative Integrals. Therefore, MASPR0 performs the 
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following two operations on the Input inertia elements to get the output 
inertia tensor; 

1. Negates the product of Inertia elements 

2. Transforms from fabrication to vehicle coordinates* 


'v • [%] 'F [\f (3) 

The Inverse of the inertia tensor Is calculated, completing the initial call. ' 

On subsequent calls to MASPR0, the current Orb Iter mass is computed as 
follows: 

m « m 0 - &" RCS - ^RUST ^ 



where 


m 

m 0 

*"RCS 

^THRUST 


Is the current vehicle mass (kg) 

Is the initial vehicle mass (kg) 

is the total mass of RCS fuel consumed (kg) 

is the total mass of main propulsion engine fuel consumed (kg) 


INPUT/OUTPUT V 

This mass properties model requires the following input parameters: 



AVMASS - Active vehicle mass; input units are pounds, which are converted to 
kilograms for simulator calculations and output 

RFX - Position of the center of mass (X) in the fabrication frame; input 
units are inches, which are converted to meters for simulator calcu 
lations and output 
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RFV - Position of the vehicle (.V) body frame In the .fabrication (F) frame; 
input units are inches, which are converted to meters for simulator 
calculations apd output 

AFV - Attitude of the vehicle (V) body frame In the fabrication ..(F) frame;* 
input units are degrees which are converted to radians for simulator 
calculations and output 

AVINRT - Active vehicle inertia tensor. Input as positive integrals in slugs* 
ft 2 in the fabrication system. Converted to negative integrals in 
kg-m 2 in vehicle system for SOAP calculation and output 

MPMAS - Mass properties model activity switch (nominally set to -1 for a 
constant mass model ) 

NOLOAD - PAYLOAD selection switch, nominally set to 0 for 32,000 pound payload, 
or input as nonzero for no payload 

* 

The above input parameters are initially input in English units and have been 

defined using data statements. However, the internal system of units for 

MASPR0 is MKS. Payload independent values are as follows: 


0.0 

rad 

(0.0°) 


3.1415927 

rad 

(180.0°) 

(5) 

0.0 

rad 

(0.0°) 



MPMASS — 1 

Values for both nominal (with payload) and off-nominal (no payload) configu- 
rations are given in the attached tables. 

I? a decremented vehicle mass is desired, i.e., to account for fuel burned by 
the RCS and main propulsion systems, the model activity switch (ENVCOM loca- 
tion 389) must be changed from -1 to fl. This is done by resetting MPMAS via 
data input. 

The above data values were taken from reference- 12. These values can be 
changed by input If desired. Common block locations may be obtained from the 
attached listing of the MASPR0 common block. 



' 

f 

j. 

I 


) 




I 
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The output from MASPR0 is; 


Parameter 

Description 

Units 

AVMASS 

Active vehicle mass 

kg 

RVX 

CM location in vehicle body frame 

m 

AVINRT 

Active vehicle inertia tensor 

kg-m 2 

AVINV 

Inverse of active vehicle inertia tensor 

T/kg-i 

RFV 

Position of V in F 

m 

AMFV 

F to V attitude matrix 

m 

AFV 

F to V attitude 

rad 

RFX 

Position of CM in F 

m 

AVMASI 

Initial active vehicle mass 

kg 


MASPR0 accepts no commands and has no calling arguments. 
PROGRAM VERIFICATION 


MASPRjtf has been verified by open loop checkout, runs. 


NOMINAL VALUES 32,000 lb PAYLOAD 


f 




z 

01 

51 

c « 

01 

I 0* XI 3 
0» <D <D O 



V> 

2 £■ 

OJ 4- 
•M 

e c 

•p o 


<M 


01 <u 
> Cd>~ 

V'P U4J 
I 3 

OVf- U Q 
3 l/l -Q C 

Q (V 

V) S.*- *— 


o 

CO 

in 

• 

• 

o 

wo 


CO 

o> 

o 

CM 

co 

00 

r* 

• 

* 

* 

CM 


p H 

CM 

1 

CO 

CM 


• 

VO 




CM 


co 

wo 

00 


O ro 

* * 


01 


o 

00 

m 




CM F- 
CO I 



o> 

UP> 

• 

00 

wo 


wo 

LO 

• 

00 

CM 


CM 

o 


o> 

00 

CM 

os 

VO 


O 

VO 

CM 

VO 

P* 

CM 

CO 

in 

o 

CO 

• 

Os 

• 

00 

CM 

• 

• 

OS 


o 

00 

cr> 

vn 

VO 

wo 


o 

«T 

• 

• 

• 

<o 

r— 

00 

os 

WO 

CM 

CO 

VO 

in 

• 


* 


o 

P- 

CO 

r* 

o 


i 

CO 

• 


:l 


01 CM O 


m 

00 

VO 

m 

V o 
m 
00 


co wo 
wo as 
co <o 
• * 
r** cm 

I CM 
CM 



II 


II 



4-40 




i U .ta£^ 4 il.uf. ■ ■ .i-A- — . 



^ f**, in in a\ in in ov 

< oo in o ^ in o 



o 



■fK 

tn 

o 

r*. 

CO 

CM 

*r 

ocv 


i n 

O 

« 

* :■ 

■ •■ 

• 

• 

• 

K 

r> 

.. .*** 

00 

*r 

CO 


CM 

00 

in 

iO 

00 

O 

CM 

, CM 

CO 

r* 

MO 

• 

•; 

•» 


» 

m 

o* 

CM 

CM 

rs, 


r— 

i 

o 

1 

i 

as 

» 


OV 

* 


00 

*■ 

in 




MASPRS 


PROGRAM DESCRIPTON 

» 

MASPRS defines the mass properties of the maneuverable television (MTV) or any 
gas-fueled vehicle. The center of mass (CM), the vehicle inertia matrix, and 
the mass are updated to model fuel consumption by the RCS jets. 

MATH MOOEL 

MASPRS calculates a center of gravity (CG), an active vehicle inertia tensor, 
and an active vehicle mass each time it is called. These calculations are 
performed after setting their initial values to zero each time. This mass 
properties model calculates the vehicle's CG position in the body frame. 

Initial parameters are input to MASPRS via data statements In MKS units. 

VEHICLE BOD'Y 



Z 


RVX is the position vector of the CG expressed in the vehicle body coordinate 
system, from the origin of the vehicle coordinate system to the actual CG. 

The inertia input is in body frame coordinates and the off-diagonal (product) 
elements are defined by negative integrals. The vehicle is broken into 
several subvehicles or pods, each with its own inertia tensor, The mass of 
each pod may change, but the relative magnitudes of the elements of the 
inertia tensor of each pod remain the same. The initial inertia tensor 
associated with each pod is also input and defined by the same integrals as in 
the active vehicle inertia tensor. As the mass of each pod changes so does 
its inertia tensor, thereby changing the entire active vehicle inertia 
tensor. Therefore MASPRS' performs two coerations to obtain the output inertia 
tensor: (1) Calculates current mass of each pod, and (2) Calculates current 

inertia tensor of each pod. 
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A total vehicle Inertia tensor is then found and the inverse of the inertia 
tensor is calculated by a call to INVRSE'. 

Therefore, on each call to MASPRS the current inertia tensor and its inverse 
are calculated. 

The current vehicle mass is also computed as follows: 

M » £ M p (1) 

where M p is the mass associated with each pod and 


M » M nn - AMn r c 
p po RCS 


( 2 ) 


where Mp 0 is the Initial mass associated with each pod and aM rc ^ is the 
mass of RCS fuel consumed by each pod. 

The position of the vehicle CM is recomputed at every pass using 


RVX = l MpRVPOO (3) 

where the vectors RVPOD are positions of the pod CM in the vehicle body 
frame. These positions are assumed to be unchanging although the mass 
associated with each pod can change. 

By convention, the pod with index 1 corresponds to those parts of the vehicle 
which do not move, while pods 2 through 5 represent RCS fuel tanks. 
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INPUT/OUTPUT 


Prog. 

Svmbol 

COMMON 

block 

name 

COMMON 

block 

location 

Dimen. 

Description (unit;) 

AVMASS 

MASPIC 

1 


Active vehicle mass (kg) 

RVX 

MASPiKC 

2 

3 

Position of CM In vehicle system coor- 
dinates (m) 

AVINRT 

MASP#C 

5 

9 

Active vehicle inertial tensor, input 
as negative integral In vehicle system 

o 

coordinates (kg-m ) 

XMASO 

MASPIC 

23 

5 

Initial mass of each pod (kg) 

RVPOD 

MASPiPC 


3,5 

CM of each pod, in vehicle body frame 
coordinates 

FUELUS 

RCS#C 

8 

50 

Total fuel used by each jet (kg) 

JETPOD 

RCS#C 

1133 

50 

An array which identifies pod membership 
of each RCS jet 

INPUT 





Prog 

Svmbol 

COMMON 

block 

name 

COMMON 

block 

location 

Dimen. 

Description (unit) 

AINRT 

MASPIC 

43 

9,5 

Inertia tensor associated with each 
2 

pod (kg-m ) 

OUTPUT 





Prog 

Svmbol 

COMMON 

block 

name 

COMMON 

block 

location 

Dimen. 

Description (unit) 

AVINV 

MASP1C 

14 


Inverse of active vehicle Inertia 


2 

tensor (1/kg-m ) 


MASPRS accepts no commands and has no calling arguments. 

. ... . * . ■ 
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The Input/output for MASPRS/MEC is identical to that for MASPRS with the 
following exceptions; 


INPUT 


Prog. 

symbol 

COMMON 

block 

name 

COMMON 

block 

location 

Dimen. 

Description (units) 

♦RADGYR 

MASPIC 

43 

i 

9,5 

Radius of gyration squared associated 
with each pod (m 2 ) 

OBJVOL 

MASPIC 

95 

1 

Volume of object (MEC without antenna) 
(m 3 ) 

SPHVOL 

MASPIC 

96 

1 

Volume of J.^herical gas bottles (m 2 ) 


♦RADGYR is defined as follows; 


iiri 


1j 

* 

where I is the moment of inertia and m is the mass of each pod 


MASPRS/MEC 


MASPRS/MEC is a mass property version similar in description to MASPRS except 

for the following changes and assumptions: 

a. Two spherical gas bottles, volumes stored in SPNVOL , and treated as POO 12 
and POD #3 

b. One antenna whose radius of gyration and mass is treated as POD #4, 
(RADGYR(J ,4) , XMASS(4)) 

c. Rest of MEC is an object of volume OBJ VOL, which includes the gas bottle 
volume but not the antenna. It is treated as POD #1, (RADGYR ( J , 1 ) » 
XMASS(l)), with the calculation for its radius of gyration using a 
homogeneous body formula. 


VERSIONS AVAILABLE 


Both single- and double-precision versions of MASPRS are available. The 
supporting routine BDMSPIC has both single- and double- precision versions, as 
well as versions for several vehicles. 
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HASP 4IC COMMON BLOCK-VERSION FOR MTV LENGTH = 90 PAGE 1 OF 

LGC VARNAM DIM UNITS DEFINITION AND INITIAL VALUE 
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4,8 REACTION CONTROL SYSTEM 
RCS 

This routine calls subroutines which model the forces and torques due to each 
vehicle In the simulation, l.e., RCSO, RCS1, RSC2 , It has no calling argu- 
ments, Input or output variables, or COMMON blocks, 

RCSO 

PROGRAM DESCRIPTION 

The RCSO subroutine provides a dynamic model of 38 Attitude Control Propul- 
sion System (ACPS) bl propellant primary thrusters <870-1 b nominal vacuum 
thrust) and six vernier thrusters (25-lb nominal vacuum thrust), These 
thrusters are arranged as specified In references 13 and 14 for the Rockwell 
International orbiter vehicle. Each of the 44 thrusters Is Individually con- 
trolled and commanded. Each type of thruster, either primary or vernier, has 
the same on-delays, off-delays, minimum Impulse, specific Impulse, total 
Impulse and fuel -loss-per-fl ring. The forces, torques, and fuel are calcu- 
lated for each engine. The resultant forces and torques, Including vacuum 
impingement Increments,* are calculated and filed for Incorporation In the * 
vehicle's dynamics. The total fuel consumed Is calculated and filed for 
updating the vehicle's mass properties, 

MATH MODEL 

* 4 

RSCG models 38 ACPS engines of 870-lb thrust each and Six vernier engines 
of 25-lb thrust each, according to the June 1976 data given In reference 15. 

With appropriate input, RSCO can model up to 50 ACPS engines which have 
constant thrust magnitude. Also, up to three different-size engines can be 
modeled. 

The force direction (see fig. and table 4-1 1) depends on the location of 
the engine on the vehicle. All engine locations are initially defined in the 
fabrication (F) frame, from which the locations of the ACPS engines in the 
vehicle's (V) body frame are calculated during Initialization. Knowing the 
vehicle’s initial center of mass (X), the torques (n^) of the engines can be 
calculated. 


4-47 


A i’ rnr II I • 



If the mass properties of the vehicle are constant, then the position l'«vx) 
of the vehicle's center of mass will be constant, and the torque ( N.) 

calculated during initialization will be constant. If the mass properties 
are not constant, then the ACPS engine torques and torque impingement incre- 
ments are recalculated each update. 

* Bfv) 

-i “ (*VJ. ’ *VX * 3-i) 

The symbols in the above equations are defined below. The program symbol is 
given in the brackets [ ]. 

Rp V - Position of the V body frame In the F frame [RFVl 

Ry X - Position of X in the V body frame [RVXl 
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TABLE 4 -II,- RCSO INTERNAL VALUES FOR JET SIZES, STATIONS, AND FORCES 
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Mote: 1. RCS Jet Stations (R e , j are given in the fabrication (F) frame, 

in inches. '~ FJ i / 

2. RCS Jet Forces (oj are given in the vehicle (V) frame. The forces 
shown do not include impingement increments. 

3. RCS Jet Stations for all forward jets (JETPOD ■ 7) include the 
distance % [DISH] from the thruster attach point to the point of 
effective thrust due to nozzle scarfing effects. 
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ORIGINAL PAGE IS 
OP POOR QUALITY 


Note: 


Transformation matrix from p frame to V frame 

■* * 

Position (station) of the ith ACPS engine (J) In the F 
frame refiled from RCSSTI [RCSSTA] 

* 

Position of J In the V body frame IftCSPOS 

Force vector ’of the ith ACPS engine refiled in newtons from RCSFI 
[RCSFOR] 

Torque vector of the Ith ACPS engine [RCSTOR] 

For the 14 main and two vernier jets In the forward jetpod, a distance 
i [OISTL(I)] (see table 4-III) Is added to ^ . This is the distance 

from the jet attach point to the effective thrust point as a result 
of nozzle scarfing effects. 

< » 

TABLE 4-III.- OISTL 


OISTL contains the distance (in Inches) from the thruster attach point (RCSSTI) 
to the effective thrust point for all forward RCS engines, the last six loca- 
tions are for the verniers, Jets 509, 510, 513, 514, 541, 542 Include cant 
angles (37°). 


JET 



ID 

X 

Y Z 

501 

- 22.88 

, 0,0 , 0.0 

502 

- 22.88 

, 0.0 , 0.0 

503 

- 22.88 

, 0.0 , 0.0 

504 

0.0 

,-15.28 , 0.0 

505 

0.0 

,-15.28 , 0.0 

506 

0.0 

, 0.0 , 16.22 

507 

0.0 

. 0.0 , 15.80 

508 

0.0 

, 0.0 , 16.22 

509 

0.0 

,-12.385,-16.436 

510 

0.0 

,-12.385,-16.436 

511 

0.0 

, 15.28 , 0.0 

512 

0.0 

, 15.28 , 0.0 

513 

0.0 

, 12.385,-16.436 

514 

0.0 

, 12.385,-16.436 

541 

0.0 

,-11.976,-15.893 

mm 

0.0 

. 11.976,-15.893 



The fuel mass consumed by the firing of an engine Is computed using a 
variable specific Impulse. This variable specif 1c. Impulse is a function of 
the electrical pulse width (EPW). EPW Is defined as 

EPW » COMMAND qff - COMMAND^ 

where 

CQMMANDqji - Time at which the engine is commanded on (does not Include an 
on-delay) 

COMMANDqpp - Time at which the engine Is commanded off (does not Include an 
off-delay) 

Since all commands received by RCSO include on and off delays, EPW is cal- 
culated by RCSO as 

EPW = TIME off - TIME qn + DONOFF 

where 

TIME qn * COMMAND on + DELAY qn 
T1ME off * COMMANDqpp + DELAY 0FF 
DONOFF » DELAY qn - DELAYq FF 

The symbols In the above equations are defined below with program symbols in 
brackets, < 

EPW - Electrical pulse width [EPW] 

TIMEqn - Time at which PCS engine actually fires [ONTIME(I)] 

TIMEq FF - Time at which RCS engine actually ceases firing [OFFTIME(I)] 

DELAYqn - On-delay for type of engine (main or vernier) [RCSONO(J)] 

DELAYq FF - Off -del ay for type of engine [RCSOFD(J)] 

DONOFF - Difference between on and off delays [OONOFD(J)] 
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If an update Is requested at some point before an engine is turned off, then 

the mass consumed up to that point Is calculated with EPW defined as 

EPW • TIMEyp - TIME on + DELAY on 

where TIMEyp * Time at which update Is requested [AVTIMEj. 

The fuel mass consumed by an engine firing is calculated as 

imp t 

aM fuel * TmpJ + m loss , 

where 

IMP t * (EPW - 0|LAY on )THR 

IMP S * f(EPW) * 

The symbols in the above equations are defined below with program symbols in 

brackets . 

AM FUEL - Delta fuel mass consumed by an engine, this firing as of this 
update [aFUEL] 

IMPj - Total impulse as a function of EPW [TIMP] 

IMP<j ■* Variable specific impulse [SPl] as a function of EPW found by 
evaluating a curve-fitted polynomial f(EPW) [SPIT8L(K,J)](see 
table 4-IV). For st'eady state firings, i.e. , EPW 1.0 second, the 
variable specific impulse is set to its steady state value taken from 
the input table SIMP. IMP S « SSIMP S [SSMSPI(K)] 

M L0SS ” l° ss P er ^iTing (dribble penalty, etc.) for each type of engine 
(main or vernier) refiled’ in kg from FUELI [FUELOS(J)] 

THR - Total centerline vacuum thrust for each type of engine (main or 
vernier) , in kilograms [THRK(J)1 


TABLE 4- IV,- SPIMP 


Variable specific impulse as a function of electrical pulse width for main and 


vernier engines. 








EPW (seconds) 

*• 

0.04 

,0.16 0.32 

0.52 

0.64 

0.80 

0.88 

1.00 

Variable) Main 
specific | 

185,0 

266.0 279.0 

285.0 

286.5 

288.0 

288.5 

289.0 

impulse ) Verniers 

240.0 

256.0 262.0 

266.7 

268.7 

271.0 

271.8 

272.0 


The total fuel mass consumed by all engines is calculated as 


M FUEL " M FUEL “ M FLAST i + AM FUEL i 


where 

M FUEL 

m flast 1 

aM FUEL. 


- Total fuel mass consumed by all firings of all engines as of this 
update [SFUEL] 

- Mass: of fuel consumed by jet i during this firing prior to this 
update [FUELST(I)] 

- Mass of fuel consumed by jet i for this firing as of this update 

[dfuel] 


The total fuel mass consumed by each engine over all its firings is main- 
tained [FUEL(I)]. Total on-times for all engines [SONT] and for each engine 
[RCSQNT(I)] are also maintained. 

In computing 'forces and moments for the aft main RCS engines, vacuum impinge- 
ment effects are included (ref. 16). The aft RCS vacuum impingement force . 
increments (AN, aY, aA) (see table 4-V) are computed using coefficients 
( N ft , Ny i N n ) (see table 4-VI)initially defined in the F frame. 

During initialization, these coefficients are converted to the V frame 
and the force increments are calculated as 

AN * N m • THR 
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TABLE 4-V.- RCS FORCE IMPINGEMENT INCREMENT 


JETPOD 

X-axis 

Y-axis 

Z-axis 

1 

131.57839 

321 .20608 

.00000 

2 

131.57839 

-321 .20608 

.00000 

3 

638.54221 

247.67698 

1180. 33559 

4 

638.54221 

-247.67698 

1180.33559 

5 

7.73991 

.00000 

61.91924 

6 

7.73991 

.00000 

61,91924 

7 

.OOOOO 

.00000 

.00000 


TABLE 4-VI RCS FORCE PLUME IMPINGEMENT EFFECT COEFFICIENTS 


i 

JETPOD 

X 

Y 

Z (Fabrication fram«) 

1 

-0.034 

0.083 

0.0 

2 

-0.034 

-0.083 

0.0 

3 

-0.165 

0.064 

-0.305 

4 

-0.165 

-0.064 

-0.305 

5 

-0.002 

0.0 

-0.016 

V.: 6 

-0.002 

0.0 

-0.016 

7 

0.0 

0.0 

0.0 





AY « N y • THR 


AA « N A • THR 


The symbols used In the above equations are defined below with program symbols 
in brackets. 


aN - Normal component (along z-axis) of force increment [RCSFM(d,3)] 
AY - Side component (along y-axis) of force increment [RCSFII(J,2)] 
aA - Axial component (along x-axis) of force Increment [RCSFII(J,1 )] 
N N - Normal component of impingement force coefficient (RCSFI I ( J ,3)] 


Ny - Side component of impingement force coefficient [RCSFII(J,2)j 
N a - Axial component of impingement force coefficient [RCSFII(J,1 )] 

THR - Centerline vacuum thrust of aft RCS engine in newtons [THRN(l)] 

Data for the aft RCS vacuum impingement moment coefficients (fL, N^, N n ) 

(see table 4-VII) were derived separately and were Initially defined in the 
V frame, thus no conversion Is necessary. During initialization, the 
moment increments (aN 4 , AN m> AN p | (see table 4-VIII) are calculated as 


AN fl » N. • A, ur • THR 
l i AVG^ 


4 W A AVG m • ™ 

(T1 


aN „ ■ N n • a /wg„ • THR 


where 


AN i - Roll component of. impingement moment increment [RCSTII(1 ,1 )] 
aN^ - Pitch component of impingement moment increment [RCSTII(I ,2)] 
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TABLE 4— VII*— = RCS TORQUE IMPINGEMENT INCREMENT 


[Newton-meters] 


JETPOD 

X-axis 

Y-axis 

Z-axIs 

1 

,2393.58603 

-271.24019 

-3308.55228 

2 ‘ 

-2393.58603 

-271 .24019 

3308.55228 

3 

-3447.64221 

12399.55164 

-1654.27614 

4 

3447.64221 

12399.55164 

1654.27614 

5 

614.86613 

503.73178 

-157.55011 

6 

-614.86613 

* 

503.73178 

157.55011 

7 

.00000 

.00000 

, 00000 


Note: These forces and torques due to plume Impingement effects are added 
directly to the forces and torques of each Individual jet In the 
appropriate jet pod. 


TABLE 4-VIII.- RCS TORQUE PLUME IMPINGEMENT EFFECT COEFFICIENTS 





dN n - Yaw component of impingement moment increment [RCSTII(I,3)] 

N^ - Roll component of impingement moment coefficient [fsCSIIT(1 ,1 )] 


N m - Pitch component of impingement moment coefficient (j&S IIT (1, 2)] 
- Yaw component of impingement moment coefficient [RCSIIT(I ,3)] 



Roll component of average moment arm for aft RCS engines [ARM(l)] 


A avg - Pitch component of average moment arm for aft RCS engines [ARM (2)] 
m 

A AVG ’ ^ aw con, P° nent °Y average moment arm for aft RCS engines [ARM(3)] 
n 

Individual RCS nominal forces [RCSFI(J)] (without impingement effects) are 
input via DATA cards and converted to newtons. Individual RCS nominal 
torques [RCSTOR(O)] are calculated at initialization. Total forces [RRCSFJ(J)] 
and torques [RRCSTJ(J)] for each jet are then computed as the sum of the jet's 
nominal force or torque and the impingement force or torque increment. It • 
should be noted that the values for the impingement -increments apply to one 
jet, i.e., the increments are those created by the firing of one RCS aft main 
engine., but they will always be equal for all the engines in a particular 
group. Thus, these increments are calculated for groups but added to the 
forces and moments of individual jets. 


Total forces [SRCSF(I)] and total torques [SRCST ( I )] generated by firing RCS 
engines are the sums of individual forces and torques over all engines which 
are on. If, during a simulation run, the mass properties of the vehicle 
change, the moment arms, impingement moment Increments, nominal moments and 
the individual total moments are recalculated. 


It is possible to exclude the impingement effects from all force and moment 
calculations by changing the input value of a flag called IMPFLG (Location 
1225 in the RCS COMMON block). Normally, IMPFLG is defined as one and 
impingement effects are Included. By changing its value to zero, no impinge- 
ment effects will be incorporated. 
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INPUT/OUTPUT 


This model has no calling arguments. However, values for three parameters 
are obtained from the mass properties routine for its initialization. They 
are: 

Ry X - Position vector of the vehicle's center of mass (X) in the vehicle (V) 
body frame (meters). [RVX] 

Rpy - Position vector of the vehicle (V) body frame in the fabrication (F) 
frame (meters), [RFV] 

V Mp - Fabrication (F) to vehicle (V) body frame attitude matrtx. [AMFV] 

All other data values necessary for this model's computations are contained 
in internal data statements (see ref. 13) which can be changed using the SOAP 
data input scheme. The input parameters to this model are: 

fcj. - Positions (stations) of the ACPS engines (J ) in the fabrication (F) 

1 frame; input units are inches, which are" converted to metars for 
simulator calculations. [RCSSTl] (See table 4- 1 1 . ) 

* Force vectors of the ACPS engines in the vehicle body frame; input 
units are pounds, which are converted to newtons for simulator 
calculations. [RCSFl] (See table 4-II.) 

AXIS - Hollerith array (used for output print only), which contains the 
force direction for each engine, e.g. , 4H - X for engine 1, 4H + Y 
for engine 5, etc. 

PCSONO - jth size ACPS engine on-delay (0.016,0.010,0.0 sec) 

% 

RCSOFD - jth size ACPS engine off -delay (0.0,0. 0,0.0 sec) 

TMINPL - jth size ACPS engine minimum impulse time (0.04,0.04,0.0 sec) 

... %' : 

FUELI - Mass loss per jth size ACPS engine firing; input units are pounds* 
which are converted to kilograms for simulator calculations and 
output. (0.0,0. 0,0.0 lb) 


JSCFLG * ACPS engine size change flag* This array identifies all ACPS engines 
with a particular size engine group (see table 4-11)* 

NJ - Number of ACPS engines that are to be used; no units* (NJ * 46} 

All of the above parameters may be changed by input. 

During initialization, six parameters are calculated: 

Ryj^ - Positions of J In the V body frame. [RCSPOSj (m) 

Nj - Torque vectors of the ACPS engines [RCSTOR], This will remain 

constant if the center of mass of the vehicle is not changing, (N*m) 

RCSFII - Force Impingement increment for each jet group » 

RCSTII - Torque Impingement Increment for each jet group 

RRCSFJ - Resultant or total force for each jet (sum of RCSFOR for jet j and 

RCSFII for the jet group of jet j) 

* 

RRCSTJ - Resultant or total torque for each jet (sum of RCSTOR for jet j and 
RCSTII for the jet group of jet j) 

* i 

The values input and calculated above are printed out during initialization. 
Table 4-21 shows the values of JSCFLG, JETPOD, R^ , and Also printed at 

initialization are the steady state and minimum impulse values for the 
variable specific impulse. 

When an ACPS engine(s) is turned on (1) or off (0), optional print is 
available. This print is controlled by a print flag (RCSPC). When RCSPC 
is 0.0, no print is received when commands are issued. When RCSPC is 
1.0, the following information is printed: 

1. Time of this ACPS action 

2. The number of ACPS engines on and which ones they are 
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3. The command being issued (501 through 540) 

. 4. The axis along which the force is directed 

5. The state (on or off) of the command 

6. The dumber of times this ACPS engine has been fired 

7. The total on-time for this ACPS engine 

8. The fuel used by this ACPS engine 

9. The total fuel used by all ACPS engines 

10. The force vector for this ACPS engine 

11. The torque vector for this ACPS engine 

12. The resultant or total force vector for all ACPS engines which are on 

13. The resultant or total torque vector for all ACPS engines which are on 

When a jet command is incorporated, the vehicle program updates its dynamics. 

A vehicle dynamics print flag, PAVRDC, can be set to foree the vehicle to 
print its response to the jet command(s), PAVRDC is defined initially by the 
vehicle model to be -1.0 via a data statement and is in the vehicle COMMON 
block (location 131 ). PAVRDC has three settings: -1.0, 0.0, and 1.0 If 

equal to -1.0, no print will be obtained and the RCS model cannot reset it, 
The value- of PAVRDC can only be changed from -1.0 by data Input. When PAVRDC 
is equal to 0.0, this RCS model will set it to 1.0 Indicating to the vehicle 
model to print its dynamics response to a command. The vehicle model resets 
PAVRDC to 0.0 for future commands. 

I ' > , 

••••••••• * _ 

RSCO uses two external subroutines, POLY and PEVAL . At initialization, POLY 
is called to construct a polynomial curve to a set of data points for variabl 
specific impulse. PEVAL is an associated routine which evaluates the poly- 
nomial for a given electrical pulse width. 


SOURCE 

The RCSO subroutine Is essentially identical to RCS13 with the following 
exceptions. 

1. The incorporation of RCS vacuum impingement force and torque increments 

2. The incorporation of a variable specific Impulse for fuel calculations 

3. New data on jet positions, forces, on and off delays, minimum Impulse 
times, and mass loss per firing 

Data for ACPS engine locations and thrust direction were taken from table 
2.4.1 .3 of reference 15. This data Is for the RI Orbiter Vehicle 5 as shown 
In Rockwell drawing number VC 70-00Q002A in reference 13. The jet numbering 
system is from reference 14. The specific impulse data is from reference 15. 
The Impingement effect data is from reference 16. The subroutines POLY and 
PEVAL are described in references 17 and 18, respectively. 

PROGRAM VERIFICATION 

The RCSO subroutine has been verified with closed-loop SOAP runs and the 
printed output agrees with hand calculations. 
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RCS 

Subroutine names RCS1, RCS2. 

PROGRAM DESCRIPTION 

These RCS models compute the force, torque, and fuel use which result from 
commands Issued to the Reaction Control System by the flight control routines. 
The RCS routines call the appropriate mass properties model to Initiate 
update of the inertia matrix, The various Subroutine names supplied allow 
the SOAP to distinguish among the models which correspond to various vehicles, 
but denote no change in algorithm. 

The RCS models maintain a COMMON block with location, thrust level, thrust 
direction, fuel container, etc., for each of up to 50 Jets, Commands to 
individual jets are issued by the flight control system along with time tags 
for the commands. Each RCS model responds only to those commands designated 
for it and ignores commands for other vehicles. 

Upon receipt of a command to either turn, a jet on or off, the RCS model 

a. Delays the action by an on or off-delay if desired. 

b. Delays off-commands until a jet has been on for some minimum time. 

c. Computes the thrust due to each jet. If desired, the jet response can 
be modeled as having a specific impulse which is a function of time. The 
routine POLY (or DPOLY) is used to generate the polynomial describing the 
variation of specific impulse and PEVAL (or DPEVAL) evaluates the impulse. 

d. Computes the fuel use associated with this jet firing. This computation 
may include the effects of the varying specific impulse noted above (so 
fuel use is not simply proportional to time on) and may include a loss of 
fuel by the jet. This loss is modeled as a fixed amount per firing. 

e. Causes the mass properties of the vehicle to be updated. 

f. Computes the torque due to each jet and the total torque about the 
current center of mass of the vehicle. 

g. Axiliary, nondynamic quantities such as which jets are on, how long each 
is on, how many times each has been used, etc., are computed. 

As currently configured, RCS1 accepts commands numbered 551 * 551 + n and RCS2 
accepts commands 601 •+ 651 + m where n and m are the number of jets specified 
to be in use. 
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INPUT/OUTPUT 

RCS has no call log arguments. 


Prog. 

Symbol 

i/o 

COMMON 

block 

name 

COMMON 

block 

location 

Dimen. 

* 

Description (unit) 

PAVRDC 

0 

ACTVIC 

131 

1 

Print flag 

AVTIME 

I 

ACTVfC 

138 

1 

Active vehicle time (sec) 

IN IT 

I 

ENVCOM 

322 

1 

Initialization flag 

NCCMDS 

I 

ENVCOM 

323 

1 

Number of commands on push list 

TPUSH 

I 

ENVCOM 

431 

50 

Time associated with command in 
CMMD (sec) 

CMMD 

I 

ENVCOM 

481 

50 

Command Identifier 

DAT I 

I 

ENVCOM 

531 

50 

* 

Command, 

1 * turn jet on, 
0 * turn jet off 

LCBRCS 

I* 

LECCOM 

1713 

1 

Length Of RCS COMMON block 

RVX 

I 

MASP#C 

2 

3 

Center of mass location in body 
fixed coordinates (m) 

SRCSF 

0 

RCSIC 

1 

3 

Sum of RCS jet forces (N) 

SRC ST 

0 

RCSIC 

4 

* 3 

Sum of RCS jet torques about 
vehicle center of mass (N*m) 

SFUEL 

0 

RCSIC 

7 

1 

Sum of RCS jet fuel used (kg) 

FUEL 

0 

RCSIC 

8 

50 

Fuel used by each jet (kg) 

NJ 

I* 

RCSIC 

58 

1 

Number of jets present 

RCSDT 

0 

RCSIC 

59 

1 

Time to which RCS is updated 

RCSSTA 

I* 

RCSIC 

61 

3xNJ 

(But 

<150) 

RCS jet locations in body coordinates 
The order is x, y, z for jet 1, x, 
y, z for jet 2, etc. , (in) 

RCSPOS 

0 

RCSIC 

211 

150 

RCS jet locations (m) 

RRCSFJ 

0 

RCSIC 

361 

150 

Forces due to RCS jets individually 
(N) 


♦Values which require initialization, either by the user or by default 
^ A symbol which depends on the vehicle involved . 
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Prog 

Symbol 

m 

COMMON 

block 

name 

COMMON 

block 

location 

Dimen. 

» 

Description (unit) 

RRCSTO 

0 

RCSIC 

511 

150 

Torques due to RCS jets Individually 
(N*m) 

CROSS 

0 

RCSIC 

661 

50 

Commanded RCS jet states, 

0 ■ off, 

1 » on 

PRCSS 

0 

RCSIC 

711 

50 

Present RCS jet states, 

0 * off, 

1 * on 

RCSPF 

0 

RCSIC 

761 

50 

RCS print flag 

QNTIME 

0 

RCSIC 

811 

50 

Time RCS jet is commanded on (sec) 

RCSONT 

0 

RCSIC 

861 

50 

Time for which RCS has been on 

FIREN 

0 

RCSIC 

911 

50 

Number of times a jet has been fired 

TMINPL 

I* 

RCSIC 

961 

3 

Minimum time for Jet to be on, 
by categories (sec) 

RCSFMG 

I* 

RCSIC 

964 

1 

RCS jet force magnitude 

SONT 

0 

RCSIC 

965 

1 

Sum of the RCS jet firing times (sec) 

RCSFLG 

0 

RCS#C 

966 

1 

RCS jet state change flag 

RCSPC 

0 

RCSIC 

967 

1 

RCS print control 

IRCSS 

0 

RCSIC 

968 

1 

Number of jets on 

IRCS 

0 

RCSIC 

969 

1 

Jet on indicator 

AXIS 

I* 

RCSIC 

1019 

50 

Direction of RCS jet force, literals 
used for prints only 

FUEL I 

I* 

RCSIC 

1071 

3 

Fuel lost per jet firing* three 
categories (kg) 

RCSOND 

I* 

RCSIC 

1074 

3 

Delay between jet-on command and 
jet firing, three categories (sec) 

RCSOFD 

I* 

RCSIC 

1077 

3 

Delay between jet- off command and 


actual jet turn off, three cate- 
gories (sec) 


COMMON COMMON 
Prog • block block 

Symbol I/O name location Pi men. Description (unit) 


DONDFD 

I* 

RCSrflC 

1080 

3 

Difference between on and off 
delays (sec) 

JSCFLG 

I* 

RCS#C 

1083 

50 

Jet size change flag 

JETPOD 

tit. 

* 

RCS#C 

1133 

50 

Jet pod membership numbers. 
Identify fuel source for each 
jet 

RCSFI 

I*/0 

RCS#C 

1226 

150 

RCS jet force vectors in the 
order x, y, z force due to 
jet 1; x, y, z force due to jrt 
2, etc. (Input in lb, changed 
to N Internally) 

RCSTOR 

0 

RCS#C 

1376 

150 

RCS jet torques (N*m) 

SPIMP 

I* 

RCS#C 

1526 

16 

Table of variable specific 
impulses (sec) 

CrFTIME 

0 

RCS#C 

1542 

50 

Time at which jets are commanded 
off (sec) 

SPI 

0 

RCS#C 

1592 

1 

V 

Specific impulse (sec) 

EPWTBL 

I* 

RCS#C 

1593 

8 

Electrical pulse width table (sec 

DFUEL 

0 

RCS#C 

1601 

1 

Change in fuel used by a jet 
this firing 

EPW 

0 

RCSfC 

1602 

1 

Electrical pulse width 

TIMP 

0 

RCS#C 

1603 

1 

Total impulse (kg/s) 

THRNI 

I* 

RCS#C 

1604 

2 

Centerline vacuum thrust of jets, 
two types (lb) 

FUELST 

0 

RCS#C 

1606 

50 

Fuel used by each jet in this 
firing prior to this update 

THRK 

0 

RCS#C 

1659 

2 

Centerline vacuum thrust (N) 

SSMSPI 

0 

RCS#C 

1661 

4 

Steady state and minimum impulse 
values for specific impulse 
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4.9 SENSOR CONTROL 
SNSC3 


ft * 


PROGRAM DESCRIPTION 

SNSC3 determines the times at which the sensor models are called and calls 
each active model at the appropriate time. Optional printout of these times 
may be selected if desired. 


MATH MODEL 

At initialization, all sensor model update times are set equal to the time to 
which the active vehicle is to update. The update times for those models 
which are not active are set to a very large number on the second pass initi- 
alization. During a normal update, the value stored as the current sensor 
control time, t $ , is the next most current update time of either the vehicle 
or a sensor model. t $ is then compared with the update times for each of the 
active sensor models to see if the time has been reached for that particular 
model to be called, i.e., 

t as tn_ 

• s m 


where tn m is the sensor model update time. The comparisons are repeated every 

update until the time which has been set for a model to be called or started 
is reached. When this time has been reached, the model is called and the time 
for the next call is determined as follows: 




+ At 


m 


where 

t m Is the current time at which the model is being called 

tn m is the model update time 
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X 


* 

« t . * 

I 

tn m Is the next time at which the model Is to be called, or the next 
update time 

At Is the model's maximum update Interval 

The two times calculated above may be printed after each update. 

INPUT/ OUTPUT 

SNSC3 uses three print control flags. If SNSPC is 1.0, a test Is made to see 
if it is time to print and if the print time, SNSPT, should be updated. If 
it is time to print and the print update interval, SNSPOT, is greater than 
0.0, SNSPT is updated before the next print. 

SNSC3 requires as input from CONST 

ROCDEL Small number used to check for roundoff 

FINITY Large number used if sensor model activity switch is off 

The model activity switches and model update intervals for all of the sensor 
models called come from ENVCOM. These models are RADALT, BARALT, LNDAID, 
RNGNAV, and AIRDAT. 

PROGRAM VERIFICATION 

SNSC3 was checked out in open-loop simulations. 
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5. FLIGHT CONTROL SYSTEMS 


The flight control system (also termed flight software) models the activities 
of each vehicle and pilot in response to both the environmental conditions 
and to some overall mission plan. Functions which are elsewhere termed 
guidance, navigation, and control are all included in this category. 

The flight control routines are arranged into blocks which model a specific 
vehicle. Typically, these blocks interact mainly within themselves. Infor- 
mation about the outside world can be obtained from the environment models, 
the pilot model, or from sensor models. Communication within a flight 
control system is typically handled by a COMMON block whose name begins with 
the letters FS. The structure of these blocks is determined by the vehicle 
involved. 

Flight control routines are either called (by other flight control routines 
within the same vehicle) or are scheduled. In the latter case, the SOAP 
executive calls each routine at the times specified. Scheduling can be 
performed at the outset of the simulation by input or can be performed by 
other routines during the course of the simulation. 

The flight control routine documentation is arranged by vehicle. It includes 
a general description of the operation of the routines which are specific 
to a vehicle. Pilot models and utility routines are excluded. 

The flight control system models available include the Maneuverable 
Television (MTV), the Detached Equipment Carrier (DEC), and the Space 
Shuttle Orbiter (SSO). 

5.1 MTV AND DEC (MMU TYPE) 

These vehicles are based on the control system of the Manned Maneuvering 
Unit described in reference 19. A block diagram of the operation of this 
control system is shown in figure 5-1. 

The pilot (a user-supplied routine) issues commands to the hand controllers 
and/or the attitude hold button. These commands are carried in COMMON 
block locations to MTVFSW, the main control system routine. MTVFSW in 
turn calls the phase plane routine PPMTV, which may issue rotational 
commands to either halt the vehicle or to move it back to zero attitude. 
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HORIZON 

SENSOR 

(HORIZ1) 





The rotational and translational commands are then combined by MTVFSW and sent 
to the routine JETSL where the appropriate choice of jets to fire is made.- 
These commands are then filed by MTVFSW with the SOAP command executive. 

The operation of the jet selection logic is simply an implementation of the 
selections shown in table 5-1, The jet locations and numbers are shown in 
figure 5-H. The properties of the jets, their locations, thrust levels, etc,, 

are not known to the flight control system. They are specified as desired 

* 

by the user. The jet numbers and locations are indicative of the situation 
for which the control system was designed, Table 5-1 displays the jets 
which are chosen in response to commands or combinations of commands. The 
switch between prime mode and backup modes Is initiated by setting jet 
string failure flags. 

Figure 5-3 displays the phase plane logic. In the region indicated by 
"Thrusters pulsing" the flight control commands a 0.008-second rotational 
command and then waits 0.32 second before performing any further action. 

There is a separate phase plane for each axis. 

The rates and attitudes which are used by the phase plane are body rates 
and attitudes. The attitudes are obtained by integrating the body fixed 
rates (from environment) and then filtering these rates. Two filters are 
applied, The limb motion filter (modeled in LIMBF) was designed to prevent 
the vehicle from responding to bodily contortions of the astronaut using the 
manned maneuvering unit. This is a first-order attitude filter with a time 
constant of 1 second. There is also a second-order f i 1 ter with u> = 13. Hz 
and s " 0.707 applied to the sensed rates. This is termed the ripple filter 
and is modeled in FILT2. The attitude associated with any body axis is 
initialized after attitude hold is initiated when the attitude rate drops 
below 0.2 deg/sec. The attitude hold on an axis is released either when 
a rotation about that axis is commanded or when attitude hold is turned off. 

The Detached Equipment Carrier control system differs in that the pitch 
and roll attitudes received by the phase plane are generated by a horizon 
sensor. The pitch rate which the phase plane receives is biased by a 
constant equal to the negative value of the orbital angular rate (i.e., 2-tr 
radians per period). Yaw is controlled by a rate gyro. The deadband limits 
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TABLE 5-1,- JET SELECTIONS 


t 


(a) X, theta (pitch), psl (yaw) logic, prime 


No. 

Command 

Prime mode 

■ 

X 0 <|> 

12 7 10 19 1 16 13 4 

1 

+ 

1 1 . 1 l 

i 


1 1 1 1 

3 

+ 

1 l 

4 

- 

1 1 

* 

5 

* 

i -i • * , ».* . 

6 

m 

i . i 

7 

+ + 

i i 

8 

+ - 

i i 

B 

- + 

i i 

□ 

- - 

i i 

D 

+ + 

l l ^ : '' 

□ 

+ - 

l l 

13 

- + 

l l 

14 

- 

l l 

B 

+ + 

1 l 

1 

+ - 

l l 

ra 

- + 

l l 

m 

■v mm 

1 1 

19 

+ + .+ 

l ... . . - 

20 

+ + - 

l 

21 

+ - + 

l 

22 

+ - - 

1 - : 

23 

- + + 

1 

24 

- + - 

: ./ ■ /;> ‘ X . .... : _ 

25 

• • + 

1 

26 

— • - 
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(b) Y, phi 


No. 


Command 


Y <f> 


8 


17 


27 

28 


29 

30 

31 

32 


+ » 


33 

34 


+ + 

+ - 


35 

36 


- + 


it-- • 


37 

38 


+ + 
+ , - 


39 

40 


41 

42 


+ + 
+ - 


43 

44 


* + 


45 

46 


+ + + 

+ + *■ 


47 

48 


+ - + 
4- - - 


49 

50 

51 

52 


- + 4- 

«• + - 
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oil), psi (yaw) logic, prime 
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are changed to 1° and 5°. Pilot commands can initiate or terminate the use 
of the horizon sensor, inertial attitude hold, and free attitude mode. 

A diagram of the control system implementation is shown in figure 5-1 , 
Subroutine names are in parentheses, The routines MTVFSW and RCALSW are 
transparent to the user but also exist as part of the flight control system. 
The map Instructions and control statements required to obtain a simulation 
using these flight control systems are available in system files as outlined 
In section 2. 

INPUT/OUTPUT 

All input/output is via the COMMON blocks FS1C and FS2C. Interfaces with 
the pilot are defined as follows: 


Prog. 

Symbol 

iZi 

COMMON 

block 

name 

COMMON 

block 

location 

Dimen. 

Description (units) 

WV 

I 

ACTVIC 

26 

3 

Body attitude rates in rad/ sec 

HDCTRL 

I 

FSIC 

5 

6 

Hand controller commands for roll 
pitch, yaw, x, y, ft where 
+1 ” positive 
-1 * negative 
0 = no command 

ISFL 

I 

FSIC 

23 

1 

Jet string fail flag 

0 * no failures 

1 = string A failed 

2 * string B failed 

ISWON 

I 

FSIC 

12 

1 

Attitude hold on flag, 
1 « on 

ISWOF 

I 

FSIC 

13 

1 

Attitude hold off flag 
1 = off 

TH 

I/O 

FSIC 

14 

3 

Array of body angles used by 
phase plane, deg 

THD 

I/O 

FSIC 

17 

3 

Array of body rates used by- 
phase plane, deg/s 

1 Denotes 1 or 

2 
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Prog. 

Symbol 

IATHDL 

IFGHSI 

NPPCMD 


Routines 

Routine 

MTVFSW 

FILT2 

LIMBF 

PPMTV 

JETSL1 

RCALSW 


HI 

COMMON 

block 

name 

COMMON 

block 

location 

Dimen. 

* 

Description (units) 

0 

FS#C 

20 

1 

Last pass attitude hold flag, 

■ 1 for first pass through phase 
plane after initiating attitude 
hold 

i 

FS#C 

70 

1 

Horizon sensor input flag. If 
1 use horizon sensor, otherwise 
use gyro 

0 

FS#C 

90 

3 

Number of phase plane commands 
issued 


In MTV flight control systems: 

Scheduled Interval 
0.16s, typical 
0.01s 
0.01s 

Called by MTVFSW 
Called by MTVFSW 
Scheduled by MTVFSW 


Purpose 

Control system main driver 

Ripple filter 

Limb motion filter 

Phase plane logic 

Jet selection logic 

Produces a minimum 
duration jet firing 
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The onorbit flight software is' configured to the Space Transportation System- 1 
(STS-1) Level C as described in reference 20. The overview of the onorbit 
Digital Autopilot (DAP) is presented in figure 5-4. Also included in the 
software package are the Attitude processor which receives Inertial 
Measurement Unit ( IMU) attitude information for further processing and the 
Universal Pointing processor which generates guidance command information to 
be used by the autopilot. For detailed information concerning these 
processors and the onorbit DAP, consult reference 20. Table 5-II is a list of 
the major flight software routines that can be correlated to figure 5-4. 

Diagrams of the operation of the DAP are given in figure 5-4 and the routines 
which model each block are named. Detailed descriptions of the algorithms 
used are available in reference 20. 

In addition to the DAP functions, the Shuttle flight control system accepts 
input and output via the routines which model the Universal Pointing Processor 
(UNPTPR) and the DAP Load Select (DAPSLT). 

The inputs to these routines correspond to inputs via the CRT display 2011 
(figures 5-5 and 5-6) and via a control panel (C3, fig. 5-7). Reference 21 
contains detailed descriptions of the switch and display settings required to 
cause the Shuttle to perform the desired maneuvers. On figures 5-5 and 5-6, 
the COMMON block locations (and in the case of fig. 5-6, the required values) 
which correspond to pilot inputs are noted. On figure 5-5, the functions 
"Land Site Update" and Vent Door Control" are displayed. These functions are 
not modeled by the SOAP. 

In addition to the switch panel and the two CRT displays, the Shuttle flight 
control also accepts input from rotational and translational hand controllers, 
Setting these inputs is discussed in the section on user input requirements. 

Communication within the Shuttle flight control system is through the COMMON 
block FSCOM. Definitions of the variables within this block correspond 



Figure 5-4.- Overview, onorbit digital autopilot. 






TABLE 5-1!*- FLIGHT SOFTWARE FUNCTIONS 


Routine 

Associated with 

Function 

Called by 

RECON 

DAP 

DAP main driver 

OFC 

AMTRAK 

DAP 

Auto maneuver track 

RECON (as required) 

TPULSE 

DAP 

Translation pulse 

RECON (as required) 

TRAACC 

DAP 

Translation acceleration 

RECON (as required) 

ROD I SC 

DAP 

Rotation discrete 

RECON (as required) 

RPULSE 

DAP 

Rotation pulse 

RECON (as required) 

ROTACC 

DAP 

Rotational acceleration 

RECON (as required) 

RCSERR 

DAP 

RCS errors 

RECON (as required) 

PPLANE 

DAP 

Phase plane 

RECON (as required) 

P1FLTR 

DAP 

State estimator 

RECON (as required) 

P2FLTR 

DAP 

State estimator 

RECON (as required) 

PVJSL 

DAP 

Jet select 

RECON (as required) 

DAPSLT 

I/O 

DAP load select 

RECON (as required) 

ATTPRC 

I/O 

Attitude processor 

Scheduled at 0.16 
sec intervals 

UNPTPR 

I/O 

Universal pointing processor 

Scheduled at 1.92 
sec intervals 

UNPTSP 

I/O 

Universal pointing specialist 

Scheduled at 1.92 
sec intervals 

OFC 


Scheduled FSW driver 

Scheduled at 0.08 
sec intervals 
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Figure 5-5.- Universal Pointing Processor display. 

(Bold numbers denote location in FSCOM associated with the CRT input position) 
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Figure 5-6.- DAP configuration spec display. 

(Bold numbers denote location in FSCOM associated with the CRT input position) 







closely to the definitions in reference 20. An extensive list of these 
definitions is given in the appendix. 


USER INPUT REQUIREMENTS 

The mission-dependent data associated with various flight control modules 
(Hoad data) is consistent with STS-1, If the user requires changes to any of 
this data, the following FSCOM locations will need to be modified: 


Item 

FSCOM location 

Jet group 

Translation pulse, DAP A 

139 

Primary 

Translation pulse, DAP B 

140 

Primary 

Rotation discrete rate, DAP A 

112 

Primary 

Rotation discrete rate, DAP A 

111 

Vernier 

Rotation discrete rate, DAP B 

113 

Vernier 

Rotation discrete rate, DAP B 

114 

Primary 

Pulse size, DAP A 

136 

Primary 

Pulse size, DAP A 

135 

Vernier 

Pulse size, DAP B 

137 

Vernier 

Pulse size, DAP B 

138 

Primary 

Attitude deadband roll, DAP A 

101 

Primary, Vernier 

Attitude deadband roll, DAP B 

104 

Primary, Vernier 

Attitude deadband pitch, DAP A 

102 

Primary, Vernier 

Attitude deadband pitch, DAP B 

105 

Primary, Vernier 

Attitude deadband yaw, DAP A 

103 

Primary, Vernier 

Attitude deadband yaw, DAP B 

106 

Primary, Vernier 

Deadband rate limit, DAP A 

119 

Primary 

Deadband rate limit, DAP B 

122 

Primary 

Deadband rate limit, DAP A 

128 

Vernier 

Deadband rate limit, DAP B 

125 

Vernier 

Jet options for pitch, DAP A 

115 

Primary 

Jet options for pitch, DAP B 

116 

Primary 

Jet options for yaw, DAP A 

117 

Primary 

Jet options for yaw, DAP B 

118 

Primary 


This information corresponds to the DAP configuration specification display as 
shown in figure 5-6 and further explained in table 5- I I I . When changing the 


TABLE 5-1 II.- DAP CONFIGURATION SPECIALIST USAGE 
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HOAD data, the user must also set the following discretes if the data changes 
are made during (rather than at the start of) the simulation. 

FSC0M(773) * 1 

This is the DAP load change discrete flag which will cause the new ILOAD data 
to be incorporated into the appropriate flight control modules. 

The following data are initialization default values for the flight control 
software: 

Item FSCOM location Default value 


DAP load A 

775 

1 

Roll discrete submode 

905 

3 

Pitch discrete submode 

906 

3 

Yaw discrete submode 

907 

O 

o 

Primary RCS 

890 

1 

X- pulse submode 

917 

2 

Y-pulse submode 

918 

2 

Z- pulse submode 

919 

1 

Manual mode 

757 

0 


This data corresponse to the panel "switches" shown in figure 5-7, To change 
any of these "switches" the following locations must be set: 


Item 

FSCOM location 

Value 

Translation high Z 

856 

1 

X normal submode 

917 

1 

Y normal submode 

918 

1 

Z normal submode 

919 

1 

Roll acceleration submode 

905 

1 

Pitch acceleration submode 

906 

1 

Yaw acceleration submode 

907 

1 

Roll pulse submode 

905 

2 
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Item 

FSCOM location 

Value 

Pitch pulse submode 

906 

2 

Yaw pulse submode 

907 

2 

Auto mode 

757 

1 

Vernier RCS 

890 

0 

DAP load B 

775 

2 


When switching from either primary to vernier or vernier to primary, the 
following discrete must also be set: 

FSC0M(774) - 1 

This will cause the I LOAD data associated with that jet group to be used in 
the software. 

Figure 5-8 depicts the QPS-2 orbit display associated with the Universal 
Pointing Processor. Referring to figure 5-8, the following user inputs are 
required in order to execute guidance maneuvers: 


Item FSCOM location Value 


Required inertial roll attitude 

1891 

deg 

Required inertial pitch attitude 

1861 

deg 

Required inertial yaw attitude 

1900 

deg 

Desired pitch body vector 

1847 

deg 

Desired yaw body vector 

1849 

deg 

Omi cron 

1859 

deg 

Maneuver task (MNVR) ) 

1857 

1 

Maneuver task ) 

1858 

1 

LVLH task (LVLH) ) 

1857 

2 

LVLH task 

1853 

1 

LVLH task ' 

\ 

1855 

1 

Rotation task (ROT) ( 

1857 

3 

Rotation task ' 

1853 

1 

Stop/attitude hold task 

1857 

3 
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Item 


FSCOM location 


Value 


•* 1 t 


Total attitude error 

1854 

1 

Total DAP error 

1854 

0 

Current roll attitude 

1892 

deg 

Current pitch attitude 

1862 

deg 

Current yaw attitude 

1901 

deg 

Required roll attitude 

1890 

deg 

Required pitch attitude 

1860 

deg 

Required yaw attitude 

1899 

deg 

Total error roll 

1850 

deg 

Total error pitch 

1851 

deg 

Total error yaw 

1852 

deg 

Roll rate 

390 

deg 

Pitch rate 

391 

deg 

Yaw rate 

392 

deg 


The rotational hand controller (RHC) and translational hand controller (THC) 
inputs are required as follows: 


Item 

FSCOM location 



Value 


Roll RHC 

394 

±1 

for 

On, 

0 

for 

off 

Pitch RHC 

395 

±1 

for 

on, 

0 

for 

off 

Yaw RHC 

396 

±1 

for 

on, 

0 

for 

off 

X THC 

413 

±1 

for 

on, 

0 

for 

off 

Y THC 

414 

±1 

for 

on, 

0 

for 

off 

Z THC 

415 

±1 

for 

on, 

0 

for 

off 


i 


The RHC defaults to the discrete rate mode and the translational THC default 
to the pulse mode in all areas. The area must select the primary or vernier 
jets (default is primary)} however, the verniers are incompatable with the 
translation manual modes. RCS rotational mode is selected on an axis-by-axis 
basis (see fig. 5-7) with mixing of modes allowed. If the RHC exceeds 
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"soft stop" and the acceleration mode was not selected, the mode selected 
becomes active (on axis Involved) when RHC returned to detent. To exceed the 
soft stop, the following location must be set: 


Item 

Location 

Value 

Roll RHC soft stop 

406 

1 to exceed, 0 not 

Pitch RHC soft stop 

ml 

1 to exceed, 0 not 

Yaw RHC soft stop 

408 

1 to exceed, 0 not 


Free drift Is in effect after the "soft stop" when returning to detent with 
attitude hold in detent. If the acceleration mode was selected, acceleration 
command is effected when the RHC is out of detent with free drift when in 
detent. If in the automatic mode and a RHC command is effected, then the 
manual mode will be automatically selected. 
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APPENDIX 
FSCOM LOCATIONS 


Program 

symbol 

Location 

Type Dimen. 

Definition 

AUTO 

1 

L 

Auto manual switch, true a auto 

AXIS 

2 

I 

Counter, set internally on 1, 2, 3 

DEGTOR 

3 

R 

Degrees to radians 

DYINIT 

4 

I 


FTOM 

5 

R 

Feet to meters 

GTMPSS 

6 

R 

Gravity, meters per sec 2 

HOLO 

7 

I 

Attitude hold flag 

I DM AT 

8 

R 9 

Identity matrix. 

INTOM 

17 

R 

Inches to meters 

LBFTON 

18 

R 

Foot pounds to newtons 

LBTOKG 

19 

R 

Pounds to kilograms 

MANUAL 

20 

L 


MVRTRK 

21 

I 


NMTOFP 

22 

R 

Newton meters to ft lb 

NUL 

23 

I 


OFF 

24 

L 

- FALSE, a mnemonic 

ON 

25 

L 

* TRUE, a mnemonic 

PI 

27 

R 

3.1415 

PR I MR Y 

30 

L 


RTODEG 

34 

R 

Radians to degrees 

SLTOLB 

35 

R 

Slugs to pounds 

TDAP 

36 

R 

Period of flight control minor cycle 

TVC 

38 

I 

Thrust vector control 

TWOPI 

39 

R 

6 .28 ♦ »« 


Program 

symbol 

Location 

Type 

Dimen. 

Definition 

ARCOMP 

89 

R 

12 

Array compensation thresholds 

ARDB 

101 

R 

6 

Attitude deadband limits RPY for 
dapload A then B 

ARDISR 

107 

R 

4 

Attitude discrete mode rates, 
nominal then vernier, A then B 

ARMVR 

111 

R 

4 

Maneuver rate 

ARNPL 

115 

I 

2 

Nominal nose/tail jets, pitch A 
then B 

ARNYWL 

117 

I 

2 

Nominal nose/tail jets, yaw A then B 

ARRLIM 

119 

R 

12 

Rotational rate limits, RPY, nominal 
then vernier, A then B 

ARROHL 

131 

R 

4 

Rotational rate limits, RPY, nominal 
then vernier, A then B 

ARROPL 

135 

R 

4 

Rotational pulse level, nominal then 
vernier, A then B 

ARTPLS 

139 

R 

2 

Translational pulse level x, y, z 
for dapload A then B 

ATERPL 

141 

R 

3 


ATI6HZ 

144 

R 

3 

Attitude increment, low roll rate, 
deg 

ATTER 

147 

R 

3 

Attitude error PYR 

DVBSEC 

150 

R 

3 


ETSEP 

153 

L 


External tank separation, not used 

HFAIL 

154 

L 

44 

JET FAILURE, TRUE - failed 

HFLCHG 

198 

L 


Change in jet failure status has 
occurred this pass if true 

IMUFLG 

199 

L 


IMU ATT GOOD FLAG 1 = FALSE 

INERTR 

200 

R 

3 

Diagonal components of vehicle 
inertial ratio xx, yy, zz 

MCACC 

243 

R 

3 

Magnitude of control accelerations 
along x, y, z 

MTRKOP 

246 

L 


Maneuver track option enable, 
1 * enabled 

NEWJON 

247 

L 


Turn-on jet command present for ith 


jet 


Program 

symbol 

Location 

Type 

Dimen. 

Definition 

PLEXT 

371 

I 


Payload extended flag 

PPACEL 

372 

R 

3 

Phase plane accelerator levels, 
depends on primary vernier, or the 
P/L extended flag 

qBMSO 

375 

R 

4 

M50 to body quaternion measured 

QBM50C 

379 

R 

4 

M50 to body quaternion required 

QBM50P 

383 

R 

4 

M50 to body quaternion display 

RCMDPL 

387 

R 

3 


REST 

390 

R 

3 

IMU derived body rate estimates about 
x, y, z 

RESTRT 

393 

L 



RHC 

394 

R 

3 

Rotational hand controller state 
-1, +1, roll, pitch, yaw 

RHCST 

397 

I 

3 

Rotational hand controller state 

ROTJC 

400 

R 

. 3 

RCS rotation command 

SFTSTP 

406 

L 

3 

Soft stop, status of RHS, each axis 

TCMD 

412 

R 


GMT of required body quaternion 

THC 

413 

I 

3 

Translational hand controller 
commands x, y, z, -1, +1 

TRCMD 

418 

R 

3 

x, y, z components of required body 
attitude rate 

VMASS 

422 

R 


Vehicle mass 

AGPRI 

425 

R 


Acceleration gain, primary 

AGVER 

427 

R 


Acceleration gain, vernier 

AG1PRI 

428 

R 


Attitude gain 1, primary 

AG1VER 

430 

R 


Attitude gain 1, vernier 

AG2PRI 

431 

P 


Attitude gain 2, primary 

AG2VER 

433 

R 


Attitude gain 2, vernier 

AMM 

434 

R 


Approximate/exact al gorithm threshol d 

RG1PRI 

468 

R 


Rate gain 1 primary 

RG1VER 

470 

R 


Rate gain 1 vernier 

RG2PRI 

471 

R 


Rate gain 2 primary 

RG2VER 

473 

R 


Rate gain 2 vernier 


Program 

symbol 

Location 

Type 

Dimen. 

Definition 

RUM 

475 

R 


Phase phase rate limit 

TRNINC 

480 

R 

42 

Translational Increments 14 jets, 
x y z for each 

T1 

522 

R 


Threshold for use of second vernier 
jet 

T2 

523 

R 


Threshold for use of third vernier 
jet 

WFRATE 

524 

R 


Scale factor for off-axis vernier 
preference 

ANGINC 

530 

R 

150 

Angle Increments 26 jets, x, y, z 
for each 0.08 sec firing 

JETMAP 

680 

R 

38 

Jet numbers 

MCAALT 

718 

R 

15 

Magnitude of control acceleration 
vernier alternate 

PR I MM I 

733 

R 


Rotational rate change due to 
minimum impluse burn with primary 
jet 

PTRM1E 

734 

R 



TNBBOD 

735 

R 

9 

Transformation to body coordinates 

TRANMI 

744 

R 



VERNMI 

745 

R 


Largest rotational rate change due 
to minimum impulse burn with vernier 
jet 

AEST1 

748 

R 

3 

Attitude estimate 1 

AEST2 

751 

R 

3 

Attitude estimate 2 

AGAIN 

754 

R 


Acceleration gain chosen 

AGAIN! 

755 

R 


Attitude gain 1 chosen 

AGAIN2 

756 

R 


Attitude gain 2 chosen 

ATOMAN 

757 

L 


Automatic or manual, AUTO = TRUE 

ATT 

758 

R 

3 

Attitude degrees 

ATTHLD 

761 

L 


Att hold flag 

BYPASS 

766 

L 

13 

Bypass' this axis 

CLCNTR 

769 

I 


Counter in RECON 
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Program 

symbol 

Location 

Type 

Dimen, 

Definition 

DAPCHG 

773 

L 

2 

Change In dapload present, first 
element true for change, second 
element true for primary/vernier 

DAPLD 

775 

I 


Dapload index - 1 * dapload A, 
2 = dap load B 

DATINC 

776 

R 

3 

Desired attitude increment 

DATT 

779 

R 

3 

Desired attitude 

DB 

782 

R 

3 

Deadband to phase plane 

DBREQ 

785 

R 

3 

Deadband requested 

OBRT 

788 

R 

3 

Desired body rate 

DISRT 

791 

R 



DVMIMP 

795 

R 

3 

Change in body angle rate due to 
minimum impulse 

DWRCS 

798 

R 

3 

Rate change due to jet firing 

FLG6HZ 

849 

L 


* 

FORFIR 

850 

L 

3 

Rate damping flags 

HILOZ 

856 

L 


Manual mode translation, z high 
regardless of T or F 

IAMTRK 

857 

I 


Initiate auto maneuver track 

INITJS 

860 

L 


Initiate jet select flag 

IP1FTR 

861 

L 


Initiate part 1 filter flag 

IRCSER 

862 

L 

3 

Initiate RCS errors flag 

IRDISC 

865 

L 


Initiate rotation discrete flag 

IRPLSE 

866 

L 


Initiate rotation pulse 

ITPLSE 

867 

L 


Initiate translation pulse 

LOPTN 

871 

L 


Tail of nose jet selection for low 
pitch rotation level, 1 * tail 
0 = nose 

LOYWTN 

872 

L 


Tail or nose jet selection for low 
yaw rotation 

LROTOP 

873 

I 

3 

Last rotation option 

MAGMVR 

876 

R 


Magnitude of auto maneuver rate 
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Program 

symbol 

Location 

Type 

Dimen. 

Definition 

PVINDX 

889 

I 


Primary or vernier jet index, 
1 3 P, 2 = V, depends on PVSW 

PVSW 

890 

L 


Primary or vernier jet switch, 
true 3 normal 

REST2 

895 

R 

3 

Rate estimate 2 

RGAIN1 

898 

R 


Rate gain 1 

RGAIN2 

899 

R 


Rate gain 2 

RLAMP 

900 

I 

3 

Rotation underway lamp 

ROHILO 

903 

L 

2 

Pitch then yaw acceleration level 
selection 1 3 nominal 

ROTOPT 

905 

I 

3 

Manual mode rotation, R, P, Y 
1 3 accel , 2 = pulse, 3 * discrete 

RPLSES 

908 

R 



RTER 

909 

R 

3 

Attitude rate error 

RTLIM 

912 

R 

3 

Attitude rate limit 

TQPT 

917 

I 

3 

Manual translation mode for x, y, z 
1 3 normal, 2 = pulsed 

TPLSES 

920 

R 


Translational pulse size, depends on 
deadband 

TRANJC 

921 

I 

3 

Translational jet command 

UNDSAC 

930 

R 

3 

Undesired acceleration, each axis 

AAOPT 

933 

I 


Auto maneuver option, 1 3 hold, 
2 = maneuver 

FCCNTR 

940 

I 


Counter in RECON 

OLAAOP 

946 

I 


Last pass auto maneuver option 

OLBYPS 

947 

L 

3 

Last pass bypass phase plane flag 

OLDAM 

950 

L 


Last pass auto maneuver change flag 

OLDTQP 

955 

L 

3 

Last pass translation HC option 

OLMODE 

958 

I 


Last pass auto or manual mode 

OLROPT 

960 

I 

3 

Last pass rotation HC option 

OLSSTP 

963 

L 

3 

Last pass soft stop flag 

AM 

976 

R 


Auto maneuver eigen angle 

BIAS 

977 

R 


Auto maneuver bias in phase plane 
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ewmhni m Location Type Dimen. Definition 


BIAS V 

978 

R 

3 

DELTA 

981 

R 

3 

MTRACK 

986 

L 


OLDMT 

987 

L 


OMGCMD 

988 

R 

3 

GCB 

991 

R 

4 

QGMD 

995 

R 

4 

VR 

999 

R 

3 

WEST 

1002 

R 


AZANG 

1006 

R 


COP 

1007 

R 


COSP 

1008 

R 


COSR 

1009 

R 


COSY 

1010 

R 


CP 

1011 

R 


CRMN 

1012 

R 


CRPL 

1013 

R 


CTR 

1014 

R 


CY 

1015 

R 


DELGA 

1016 

R 

4 

DP 

1020 

R 


GA 

1021 

R 

4 

GAPREV 

1025 

R 

4 

GAREAD 

1029 

R 

4 

GATOBD 

1033 

R 

4 

IATPRC 

1045 

L 


IRANG 

1047 

R 


MROLLB 

1061 

R 

9 


Auto maneuver bias vector 
Not used 

Auto maneuver flag, T * auto, 

F * hold 

Last pass auto maneuver flag 
Commanded auto maneuver rate 
Commanded quaternion 
Commanded quaternion 
Auto maneuver rotation vector 
Auto maneuver parameter 
Compensated azimuth gimbal angle 
Cosine 1/2 DP 
Cosine of pitch 

Cosine of roll ' 

Cosine of yaw 
Cosine of gimbal angle 3 
Cosine of gimbal angle 2 
Cosine of gimbal angle 1 
Inner loop, outer loop flag 
Cosine of gimbal angle 4 

Change in gimbal angles 
IMU roll /pitch nonorthogonality 

Gimbal angles 
Previous gimbal angles 
Gimbal angles relative to ADI 
Gimbal angles to body 
Initiate attitude processor flag 

Compensated inner roll gimbal angle 
IMU roll to body system transformation 

Holder for old gimbal angle sines 


OLDSIN 1070 R 
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Program 

symbol 

Location 

Type 

Dimen. 

Definition 

ORANG 

1071 

R 


Compensated outer roll glmbal angle 

PANG 

1072 

R 


Compensated pitch glmbal angle 

QBPREV 

1073 

R 

4 

Previous quaternion body to M50 

QBROLL 

1077 

R 

4 

Body roll vyrt glmbal angle quaternion 

QROLSM 

1081 

R 

4 

Roll wrt stable member quaternion 

QSMM50 

1085 

R 

4 

Stable member to M50 quaternion 

QUPD 

1089 

R 

4 

Update quaternion 

SOP 

1093 

R 


Sine 1/2 DP 

SINP 

1094 

R 


Sine of pitch 

SINR 

1095 

R 


Sine of roll 

SINY 

1096 

R 


Sine of yaw 

SP 

1097 

R 


Sine of gimbal angle 3 

SRMN 

1098 

R 


Sine of gimbal angle 2 

SRPL 

1099 

R 


Sine of gimbal angle 1 

SY 

1100 

R 


Sine of gimbal angle 4 

TCLM50 

1101 

R 

9 

Transform cluster to M50 

TNBROL 

1110 

R 

9 

Transformation, roll to body 

TNOW 

1116 

R 


Time now, sec 

ARE SI 

1119 

R 

3 

Rate estimate 1 

ARES2 

1122 

R 

3 

Rate estimate 2 

ATTINC 

1125 

R 

3 

Attitude Increment 

DTOMGA 

1128 

R 

3 

Rate estimate variable (local) 

MSDATT 

1131 

R 

3 

Rate estimate variable (local) 
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Program 

symbol 

Location 

Type 

Dimen. 

Definition 

TMEAS 

1134 

R 


Time 

TMSSQ 

1135 

R 


DAP time squared 

C 

1136 

R 


DELTHA 

1137 

R 

3 


DRSTRT 

1138 

L 



JETCMD 

1141 

R 



OLOOR 

1142 

L 

3 

> Local variables In phase plane 

OLDFF 

1145 

L 

3 


SUD 

1148 

R 



Sll 

1149 

R 



TEMP 

1150 

R 



XXIS 

1151 

I 



XI 

1152 

R 



X2 

1153 

R 



CHECK1 

1174 

L 



CHECK2 

1175 

L 


Two jet failures In nose group 
if true 

CHECK3 

1176 

L 



CHECK4 

1177 

L 



CHECK5 

1178 

L 

4 


CHECKS 

1182 

L 

4 


CHECK7 

1186 

L 



CHECK8 

1187 

L 



COMCMD 

1188 

R 

3 
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Program 

symbol 

Location 

Type 

01 men. 

Definition 

OVRCS 

1191 

R 

3 


DYFAIL 

1194 

L 

2 

Denotes failure of both y jets 
for a forward cluster If true 

DZFAIL 

1196 

L 



DZFL2 

1197 

L 


Denotes two failed -z jets on a 
side If true 

HFLJSL 

1204 

L 

38 


HIPIT 

1242 

L 


High pitch, true * high rate 

HIROLL 

1243 

L 


High roll, true * high rate 

HIX 

1244 

L 


High x rate, true = high rate 

HI YAW 

1245 

L 


High yaw, true ■ high rate 

HIZ 

1246 

L 


High z translation command present 

I 

1247 

R 



J 

1248 

R 



JETTAB 

1249 

I 

38 

Divides jets Into 14 bunches 

JONJSl 

1287 

L 

38 


K 

1325 

I 

) 


L 

NODUMP 

1326 

I 


, Local variables jet select logic 
(JSL ) 

1327 

L 

) 

NOM 

1328 

L 


Nominal 

NOSECK 

1329 

L 

6 

Nose check z, true * low nose 
rotation, false see TAILCK 

OLDJON 

1330 

L 

6 


OLDRJC 

1336 

I 

3 


RACCUM 

1358 

R 

3 \ 

i 

RCSDMP 

1361 

L 

; 


ROTMI 

1362 

R 

3 | 

f 

ROTVAL 

1363 

R 


} Local variables In JSL 

SUMDV 

1366 

R 

3 I 

( 

SZFAIL 

1369 

L 

. 3 . 


SZFL 

1372 

L 

V 

1 



1 



Program 

symbol 

Location 

Type 

Dimen. 

Definition 

TACCUM 

1373 

R 

6 


TAILCK 

1379 

L 


Tall check z, true - low rotation 
tall, else, see FSSR 

THRENZ 

1380 

L 



TWOZFL 

1381 

L 



VFAIL 

1382 

L 


Denotes fallure(s) exlst(s) In 
vernier jets 

VRNCMO 

1383 

R 

3 


XFEED 

1386 

L 


Crossfeed Interconnect flag, 
true * crossfeed open 

XFMINV 

1387 

R 

3 


ZFAIL 

1390 

L 

5 


OLDTHC 

1398 

I 

3 

Old translation HC command 

ORHCST 

1404 

I 

3 

Old rotational HC command 

FL 

1603 

L 

46 s 


JETFLG 

1649 

I 

2 1 
i 

Variables that Interface FSW 
> with RCS model 

JPFLG 

1651 

I 



P5501 

1652 

R 

46 ' 


RFLG 

1698 

I 

3 

Region flag In phase plane 

SSFSNO 

1701 

I 

44 

RCS jet numbers for RCS model 

ANG 

1836 

R 

3 

Current inertial angles 

ARBBCS 

1839 

R 



ATTERR 

1840 

R 

3 

Attitude error 

ATTRAT 

1843 

R 

3 

Rate error 

BODVP 

1846 

R 



BODVPD 

1847 

R 


Pitch component of body, desired 

BODVY 

1848 

R 



BODVYD 

1849 

R 


Yaw component of body* desired, 
deg 

ERR 

1850 

R 

3 

Total error in degrees, R, P, Y 

IBDVD 

1853 

R 


Total attitude error 

IERSEL 

1854 

R 


Error select flag 


A- 11 


; n . ji ■ i> fc^»y . ■'-&'<■ • • y • • r ' ‘ : i>7¥V» •>...'••• ' * •***.*-*...#> -- 


Program 

symbol 

Location 

Type 

Dimen,, 

Definition 

IOMID 

1855 

R 


Task flag 1 * LVLH 

I0PACT 

1856 

R 


1 roll or yaw, auto maneuver, 

2 * WLH alignment 

IOPSEL 

1857 

R 


Task flag, 1 * MNVR, 2 « LVLH, 3 » rf 

MNVD 

1858 

R 


Maneuver option discrete 

OMICRN 

1859 

R 


Oml cron 

PITCHC 

1860 

R 


Required pitch inertial 

PITCHI 

1861 

R 


Initial input for pitch (ADI) 

PITCHN 

1862 

R 


Current pitch 

PLOS 

1863 

R 

3 

Desired rotation vector 

QADIM 

1866 

R 

4 

Quaternion ADI to Inertial 

QBAOI 

1870 

R 

4 

Quaternion body to ADI 

QBADIC 

1.874 

R 

4 

Quaternion body to ADI commanded 

QBBC 

1878 

R 

4 

Quaternion body to body commanded 

QMAOI 

1882 

R 

4 

Inertial to ADI quaternion 

QMNVR 

1886 

R 

4 

Commanded quaternion 

ROLLC 

1890 

R 


Required (roll, Inertial) 

ROLL I 

1891 

R 


Initial input for roll (ADI) 

ROLLN 

1892 

R 


Current roll 

RVMAG 

1893 

R 

3 

Attitude error magnitude 

TOTERR 

1894 

R 

3 

Total error 

YAWC 

1899 

R 


Required yaw, inertial, deg 

YAWI 

1900 

R 


Yaw required. Inertial, deg 
(ADI) 

YAWN 

1901 

R 


Current yaw 

EULER 

2055 

R 

3 

Euler angle array 

CNTACA 

2073 

I 


Control acceleration index 

CNTACB 

2074 

I 


Control acceleration index 




Program 

symbol 

Location 

Type 

Olmen. 

Definition 

MCACCP 

2075 

R 

3 

Mean acceleration level, primary 

MCACCV 

2078 

R 

3 

Mean acceleration levels, vernier 

MCACPL 

2081 

R 


Mean apceleratlon levels with 
payload extended 

CNTLACC 

2096 

I 


Control acceleration index 
* CNTACA for dapload A i 
■ CNTACB for dapload B 

OLIM 

2172 

R 


1° tolerance factor 

DOTT 

2173 

R 


Local variable (dot product) 

EA 

2174 

R 


Angle of rotation 

INI OPT 

2176 

I 


Task flag 1 = LVLH, 2 * rotation 

INPASS 

2177 

I 


Initial pass flag 

IOPPRC 

2178 

I 


Task flag 1 = LVLH, 2 = rotation 

MTP 

2179 

R 

9 

Transformation matrix M50 to 
required body 

MTT 

2188 

R 

9 

Local variable matrix 

ONINTY 

2197 

R 


90° 

PLOSA 

2198 

R 

3 

Desired pointing vector 

PTM 

2201 

R 

9 

Required body to M50 
transformation 

QATM50 

2210 

R 

4 

Commanded quaternion 

QBAM50 

2214 

R 

4 

Commanded quaternion 

QMNVRA 

2218 

R 

4 

Commanded quaternion 

QRBB 

2222 

R 

4 

Local variable 

RAVGAB 

2226 

R 


Magnitude of position vector 

ROLL 

2227 

R 


90° 

RRABOD 

2228 

R 

3 

Unit vector on y-axis 

RRATA 

2231 

R 


Desired maneuver rate 

RRBOO 

2232 

R 

3 

Local variable 

RRM50 

2235 

R 

3 

Local Variable 

RXVAVG 

2238 

R 

3 

Local variable 

SEA 

2241 

R 


Local variable 


c~ y 




Program 

symbol 

Location 

Type 

Dimen. 

Definition 

TLOS 

2142 

R 

3 

Unit position vector 

TLXVAV 

2245 

R 

3 

Local variable 

UVMXRM 

2248 

R 

3 

Loc.il variable 

VBXRRB 

2251 

R 

3 

Ltcal variable 

VBXYN 

2254 

R 

3 

Local variable 

VECllOO 

2257 

R 

3 

Desired pointing vector 

VECM50 

2260 

R 

3 

Unit position vector 

VMXRRM 

2263 

R 

3 

Local variable 

VHXYT 

2266 

R 

3 

Local variable 

VVMXRM 

2269 

R 

3 

Local variable 

MB I 

2272 

R 

3 

Inertial rotation rate 

YN 

2275 

R 

3 

Local variable 

YT 

2278 

R 

3 

Local variable 

YVEC 

2281 

R 

3 

Unit vector along Y{“ 0, 1, 0) 

ZVEC 

2284 

R 

3 

Unit vector along Z(® 0, 0, 1) 

SPAV6G 

2287 

R 

3 

Shuttle position M50, inertial 
single precision 

SVAVGG 

2290 

R 

3 

Shuttle velocity, M50, inertial 
single precision 

UVVXRM 

2293 

R 

3 

Local variable 

CROLL 

2296 

ft 


Cosine of roll 

SROLL 

2297 

R 

3 

Sine of roll 

SVAVGN 

2298 

R 

3 

Unit vector of velocity 


